Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,680 G→T 24.7% L1253M (CTG→ATG)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,6800G→T24.7% 238.2 ‑0.5 ‑0.6316L1253M (CTG→ATG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (41/37):  ref base (115/123):  total (156/160)
Fisher's exact test strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.61e-01

Only 100 of 331 total aligned reads displayed.

GTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGT  >  NC_007779/3451580‑3451780
                                                                                                    |                                                                                                    
 tgtATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAgg                                                                                                     <  1:8945261/101‑1 (MQ=255)
     tgctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGt                                                                                                 <  1:8457017/101‑1 (MQ=255)
      gctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTa                                                                                                <  1:10341831/101‑1 (MQ=255)
      gctgcGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTa                                                                                                <  1:2139730/101‑1 (MQ=255)
         gcgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGa                                                                                             >  1:11889021/1‑101 (MQ=255)
         gcgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGa                                                                                             >  1:2683553/1‑101 (MQ=255)
          cgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGaa                                                                                            <  1:1069739/101‑1 (MQ=255)
          cgcCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGaa                                                                                            >  1:3714459/1‑101 (MQ=255)
            ccGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAAcc                                                                                          <  1:5899666/101‑1 (MQ=255)
            ccGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAAcc                                                                                          <  1:2292947/101‑1 (MQ=255)
            ccGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAAcc                                                                                          >  1:3206019/1‑101 (MQ=255)
              gTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCgg                                                                                        >  1:176965/1‑101 (MQ=255)
              gTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCgg                                                                                        >  1:11268666/1‑101 (MQ=255)
                aTGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTg                                                                                      <  1:6209952/101‑1 (MQ=255)
                   cTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGaa                                                                                   <  1:4069545/101‑1 (MQ=255)
                   cTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGaa                                                                                   <  1:4491305/101‑1 (MQ=255)
                   cTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGaa                                                                                   >  1:8998444/1‑101 (MQ=255)
                       cAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAAcgcg                                                                               >  1:3508308/1‑101 (MQ=255)
                         ttgcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAAcgagcg                                                                             <  1:9646473/99‑1 (MQ=255)
                          cgcCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGt                                                                            <  1:14073327/101‑1 (MQ=255)
                           gcCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTg                                                                           >  1:7441142/1‑101 (MQ=255)
                            cccACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGc                                                                          <  1:11920433/101‑1 (MQ=255)
                             ccACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTgaa                                                                         >  1:2661526/1‑99 (MQ=255)
                              cacaCTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCAt                                                                        >  1:6399639/1‑101 (MQ=255)
                               acacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATc                                                                       >  1:13906044/1‑101 (MQ=255)
                               acacCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATc                                                                       >  1:13340773/1‑101 (MQ=255)
                                 acctCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGGGAATCtt                                                                     >  1:1470422/4‑101 (MQ=255)
                                 acTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCtt                                                                     >  1:12609409/1‑101 (MQ=255)
                                 acTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCtt                                                                     <  1:6790628/101‑1 (MQ=255)
                                   ctCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCttct                                                                   <  1:11111831/100‑3 (MQ=255)
                                    tCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtc                                                                  <  1:11065400/101‑1 (MQ=255)
                                      cATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcgt                                                                >  1:2730966/1‑101 (MQ=255)
                                       aTCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTgtcgtc                                                               >  1:1223414/1‑101 (MQ=255)
                                        tctcACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCtcgtcg                                                              <  1:13629124/101‑1 (MQ=255)
                                        cccccccAAACCGTTGACCCCCTATCTGTGCCTTACCACCCAGCGGCTGCTGAGTACCCAGGCTGTAAAAACCGGTGGAACGCGCGTGCATCTTGtcgtcg                                                              <  1:9391561/98‑1 (MQ=255)
                                           ccccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGAcc                                                           <  1:12328671/101‑1 (MQ=255)
                                            cccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCa                                                          >  1:5829370/1‑101 (MQ=255)
                                            accGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCa                                                          <  1:331541/100‑1 (MQ=255)
                                             ccGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAg                                                         <  1:4328130/101‑1 (MQ=255)
                                              cGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAgg                                                        <  1:10925293/101‑1 (MQ=255)
                                                 aaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                     <  1:5324682/101‑1 (MQ=255)
                                                 aaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                     >  1:9515228/1‑101 (MQ=255)
                                                 aaCGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAggtgg                                                     <  1:2346132/101‑1 (MQ=255)
                                                   cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTTGTCGACCAGGTGGtt                                                   <  1:13822522/101‑1 (MQ=255)
                                                   cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGtt                                                   <  1:12538511/101‑1 (MQ=255)
                                                   cGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGtt                                                   >  1:4128009/1‑101 (MQ=255)
                                                    gCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTc                                                  >  1:2404763/1‑101 (MQ=255)
                                                     cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                 >  1:12474386/1‑101 (MQ=255)
                                                     cTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCa                                                 >  1:1988223/1‑101 (MQ=255)
                                                      tGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAg                                                >  1:6279678/1‑101 (MQ=255)
                                                       gACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGt                                               >  1:7545805/1‑101 (MQ=255)
                                                       gACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAgc                                               >  1:3175293/1‑100 (MQ=255)
                                                       gACCACCGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAgc                                               <  1:7309071/101‑2 (MQ=255)
                                                          caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTc                                            <  1:2924729/101‑1 (MQ=255)
                                                          caccGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTc                                            <  1:5638410/101‑1 (MQ=255)
                                                           accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                           <  1:916770/101‑1 (MQ=255)
                                                           accGAACTGAGCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCa                                           <  1:10066557/101‑1 (MQ=255)
                                                             cGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGc                                         <  1:2536456/101‑1 (MQ=255)
                                                              gAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGca                                        <  1:6716866/101‑1 (MQ=255)
                                                               aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                       >  1:8827308/1‑101 (MQ=255)
                                                               aaCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcat                                       <  1:8015771/101‑1 (MQ=255)
                                                                aCTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatg                                      <  1:13119351/101‑1 (MQ=255)
                                                                  tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgta                                    <  1:3826718/101‑1 (MQ=255)
                                                                  tgtgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgta                                    <  1:7651485/101‑1 (MQ=255)
                                                                   gagCCTTACCACCCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGcatgtac                                   <  1:8720197/99‑1 (MQ=255)
                                                                    tgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtaca                                  <  1:296537/101‑1 (MQ=255)
                                                                    tgCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtaca                                  >  1:787690/1‑101 (MQ=255)
                                                                        ttACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGcatgtacatgta                              >  1:1285288/1‑101 (MQ=255)
                                                                         taacaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTaa                             >  1:13405817/4‑101 (MQ=255)
                                                                           ccaccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGcc                           >  1:179145/1‑101 (MQ=255)
                                                                           ccaccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAAcc                           >  1:8296005/1‑101 (MQ=255)
                                                                            caccCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccc                          <  1:2084606/101‑2 (MQ=255)
                                                                            caccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCa                          >  1:12397717/1‑101 (MQ=255)
                                                                            caccCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCa                          >  1:1576508/1‑101 (MQ=255)
                                                                             accCAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGccca                         <  1:5349165/101‑3 (MQ=255)
                                                                             accCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCaa                         >  1:6000827/1‑101 (MQ=255)
                                                                              cccAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAc                        >  1:6089372/1‑101 (MQ=255)
                                                                               ccAGCGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAcc                       <  1:13184071/101‑2 (MQ=255)
                                                                               ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                       <  1:7223894/101‑1 (MQ=255)
                                                                               ccAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACg                       >  1:3994807/1‑101 (MQ=255)
                                                                                cAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCTTGTACATGTAACCAACgg                      >  1:8031205/1‑101 (MQ=255)
                                                                                 aGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGt                     <  1:10508078/101‑1 (MQ=255)
                                                                                  gCGGCTGCTGAGTAACCAGGCTGTAAGAGCCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTGCAGCATGTACATGTAACCAACgggt                    <  1:1610817/101‑3 (MQ=255)
                                                                                  gCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGtt                    <  1:834280/101‑1 (MQ=255)
                                                                                   cGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                   >  1:633419/1‑96 (MQ=255)
                                                                                   cGGCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaa                   >  1:5074918/1‑96 (MQ=255)
                                                                                    ggCTGCTGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaac                  <  1:8185495/101‑7 (MQ=255)
                                                                                     gctgctGGGTAATCATGGAGTACGGACCGGTGGAACGAGCGTGAATCTTCTCGTCGACCAGGTGGTGCAGCTTCAGCATGTACATGTAGCCCAccgaaacc                 <  1:6999309/101‑8 (MQ=255)
                                                                                     gctgctGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTAcc                 <  1:11998878/101‑1 (MQ=255)
                                                                                      ctgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCg                <  1:11441400/101‑1 (MQ=255)
                                                                                       tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg               <  1:2872524/101‑1 (MQ=255)
                                                                                       tgctgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGTTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCgg               >  1:4734138/1‑101 (MQ=255)
                                                                                          tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg            <  1:2200984/101‑1 (MQ=255)
                                                                                          tgAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACg            <  1:13880702/101‑1 (MQ=255)
                                                                                            ggTAGCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCt          <  1:13142133/100‑1 (MQ=255)
                                                                                            aGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTGACCAACGGTTACCGGACGCt          <  1:11157402/101‑1 (MQ=255)
                                                                                              tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg        <  1:3914067/101‑1 (MQ=255)
                                                                                              tAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCg        >  1:7099078/1‑101 (MQ=255)
                                                                                                aCCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGaa      <  1:646893/101‑1 (MQ=255)
                                                                                                  cAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACt    <  1:9196108/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGCTGAGTAACCAGGCTGTAAGAACCGGTGGAACGCGCGTGCATCTTGTCGTCGACCAGGTGGTTCAGTTTCAGCATGTACATGTAACCAACGGTTACCGGACGCTCGAACTGT  >  NC_007779/3451580‑3451780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.