Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,460,558 Δ1 bp 7.7% coding (914/1167 nt) algJ ← alginate O‑acetylase AlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,460,5580T. 7.7% 237.1 / 13.6 78coding (914/1167 nt)algJalginate O‑acetylase AlgJ
Reads supporting (aligned to +/- strand):  ref base T (45/27);  new base . (2/4);  total (47/31)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGG  >  NC_002947/1460412‑1460707
                                                                                                                                                   |                                                                                                                                                     
ttGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATAGCTTTGCAGGTACTTGAGCAGCGTCGGTAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCtgctgcc                                                                                                                                                    >  2:825842/1‑149 (MQ=255)
   gCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgcg                                                                                                                                                 >  1:752937/1‑150 (MQ=255)
   gCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCAGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgcg                                                                                                                                                 >  1:264811/1‑150 (MQ=255)
      gATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgcgcgc                                                                                                                                              >  1:354027/1‑150 (MQ=255)
      gATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgcgcgc                                                                                                                                              >  1:1556323/1‑150 (MQ=255)
          aCGTGCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCATCAGGCCCTGCGCAgcgacgtcgc                                                                                                                                          <  2:942562/150‑8 (MQ=255)
            gTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCa                                                                                                                                                    >  1:1076883/1‑138 (MQ=255)
            gTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCa                                                                                                                                                    <  2:1076883/138‑1 (MQ=255)
                 gggAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTc                                                                                                                                   >  2:1764856/1‑150 (MQ=255)
                 gggAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTc                                                                                                                                   >  2:1111381/1‑150 (MQ=255)
                     aTTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAAGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGt                                                                                                                               >  2:1828804/1‑150 (MQ=255)
                         ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGCGCAgc                                                                                                                                                  <  2:1051267/127‑1 (MQ=255)
                         ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGCGCAgc                                                                                                                                                  >  1:1051267/1‑127 (MQ=255)
                          cacaCCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgc                                                                                                                                                  >  1:465957/1‑126 (MQ=255)
                          cacaCCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAgc                                                                                                                                                  <  2:465957/126‑1 (MQ=255)
                          cacaCCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCggg                                                                                                                          >  1:722164/1‑150 (MQ=255)
                           acacCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGt                                                                                                                         >  1:1857778/1‑150 (MQ=255)
                            caccaccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGtt                                                                                                                        <  1:1892014/150‑1 (MQ=255)
                               caccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCGGGGCCTGCTGCAGCGCGCCAATTAAGGTCCAGTGCGGGTTGgc                                                                                                                     >  2:207752/1‑150 (MQ=255)
                               caccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAACATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCCAGGAAGTTCCAGTGCTGGTTGgc                                                                                                                     >  1:221948/1‑150 (MQ=255)
                                 ccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGgcgc                                                                                                                   >  1:1873917/1‑150 (MQ=255)
                                 ccacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGgcgc                                                                                                                   >  2:1223291/1‑150 (MQ=255)
                                  cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCt                                                                                                                  >  1:1082930/1‑150 (MQ=255)
                                  cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCt                                                                                                                  <  1:1579569/150‑1 (MQ=255)
                                          gtgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTgg                                                                                                          <  1:1223291/150‑1 (MQ=255)
                                           tgCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGt                                                                                                         >  1:501315/1‑150 (MQ=255)
                                             cggcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGc                                                                                                       <  2:136267/150‑1 (MQ=255)
                                               gcgTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCt                                                                                                     >  1:952976/1‑150 (MQ=255)
                                                  ttCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTAcc                                                                                                  >  1:1270046/1‑150 (MQ=255)
                                                   tCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCa                                                                                                 <  2:667718/150‑1 (MQ=255)
                                                         aGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCaa                                                                                           >  2:1231617/1‑150 (MQ=255)
                                                           gCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAcc                                                                                         >  1:1107145/1‑150 (MQ=255)
                                                           gCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAcc                                                                                         >  1:741210/1‑150 (MQ=255)
                                                             aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGTAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCgg                                                                                       >  2:980561/1‑150 (MQ=255)
                                                             aTCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCgg                                                                                       <  1:1415256/150‑1 (MQ=255)
                                                                         gTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGTAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTc                                                                           <  2:542692/150‑1 (MQ=255)
                                                                          tACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCg                                                                          <  1:1231617/150‑1 (MQ=255)
                                                                          tACTTGACCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGCTGGCGCTGTAGCTGGGTGGTACCAGGGCAACCGGGATCTGCTTGTCg                                                                          >  2:1753931/1‑150 (MQ=255)
                                                                            cTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGc                                                                        <  2:741210/150‑1 (MQ=255)
                                                                             ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAATGAAGTTCCAGTGCAGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                       >  2:1035322/1‑150 (MQ=255)
                                                                             ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                       >  2:764658/1‑150 (MQ=255)
                                                                             ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                       >  2:368383/1‑150 (MQ=255)
                                                                             ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                       >  2:1923282/1‑150 (MQ=255)
                                                                              tGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGa                                                                      >  2:1473208/1‑150 (MQ=255)
                                                                                 gCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGaaca                                                                   >  1:208535/1‑150 (MQ=255)
                                                                                 gCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGaaca                                                                   >  1:140062/1‑150 (MQ=255)
                                                                                    tCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAAcg                                                                >  2:717395/1‑150 (MQ=255)
                                                                                    tCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAAcg                                                                >  2:172493/1‑150 (MQ=255)
                                                                                       gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGCGCAgc                                                                                                                                                  >  1:165451/1‑65 (MQ=25)
                                                                                       gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGCGCAgc                                                                                                                                                  <  2:165451/65‑1 (MQ=25)
                                                                                       gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGt                                                             >  1:1731058/1‑150 (MQ=255)
                                                                                       gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGt                                                             >  2:761188/1‑150 (MQ=255)
                                                                                       gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAAACGGGATCTGCTTGTCGGCGAACAACGCGt                                                             >  1:39355/1‑150 (MQ=255)
                                                                                        cagcagCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTc                                                            <  1:717395/150‑1 (MQ=255)
                                                                                          gcagcGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTg                                                                                                        >  1:714384/1‑104 (MQ=255)
                                                                                          gcagcGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTg                                                                                                        <  2:714384/104‑1 (MQ=255)
                                                                                                 gCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTc                                                   <  2:1082930/150‑1 (MQ=255)
                                                                                                   cATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCAc                                                 >  2:1623771/1‑150 (MQ=255)
                                                                                                       gCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTc                                             >  1:533267/1‑150 (MQ=255)
                                                                                                       gCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTc                                             >  1:1202123/1‑150 (MQ=255)
                                                                                                         cGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCgg                                           <  1:139157/150‑1 (MQ=255)
                                                                                                         cGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCgg                                           <  2:752937/150‑1 (MQ=255)
                                                                                                                ggCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAc                                    <  1:1111381/150‑1 (MQ=255)
                                                                                                                 gCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAcc                                   >  2:247225/1‑150 (MQ=255)
                                                                                                                    tAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGAccggc                                <  2:1731058/150‑1 (MQ=255)
                                                                                                                         ggCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGggt                           <  1:1764856/150‑1 (MQ=255)
                                                                                                                            gACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGggtggt                        <  2:1880536/150‑1 (MQ=255)
                                                                                                                             aCGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGgtggtg                       <  2:952976/150‑1 (MQ=255)
                                                                                                                             aCGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGgtggtg                       >  1:337366/1‑150 (MQ=255)
                                                                                                                                  gCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTAGCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGttt                  <  1:1771817/150‑1 (MQ=255)
                                                                                                                                   cTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTc                 <  2:264811/150‑1 (MQ=255)
                                                                                                                                     gcgcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTg               <  1:761188/150‑1 (MQ=255)
                                                                                                                                       agcAGGCCCTGCGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCa             <  1:474043/149‑1 (MQ=255)
                                                                                                                                        gcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCa             <  2:337366/150‑1 (MQ=255)
                                                                                                                                        gcAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCa             <  2:1243240/150‑1 (MQ=255)
                                                                                                                                          aGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAgg           >  2:437004/1‑150 (MQ=255)
                                                                                                                                          aGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAgg           >  2:1066830/1‑150 (MQ=255)
                                                                                                                                           gggCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTGTCTGCAGGt          >  1:248502/1‑150 (MQ=255)
                                                                                                                                            ggCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGtt         <  1:1923282/150‑1 (MQ=255)
                                                                                                                                               ctgctgCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTc      >  2:1148420/1‑150 (MQ=255)
                                                                                                                                                  ctgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCgg   <  2:85688/150‑1 (MQ=255)
                                                                                                                                                   tgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTTGTGCGTTTCTGCAGGTTGTCCggg  <  2:1282288/150‑1 (MQ=255)
                                                                                                                                                   tgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCggg  <  1:764658/150‑1 (MQ=255)
                                                                                                                                                   tgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCggg  <  1:368383/150‑1 (MQ=255)
                                                                                                                                                   tgcAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCggg  >  2:163776/1‑150 (MQ=255)
                                                                                                                                                   tgcAGCGCGCCAAGGAAGTCCCAGTGCGGGATGGCGCTGTAGCTTGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCggg  <  2:208535/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTTGTCCGGG  >  NC_002947/1460412‑1460707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: