Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,460,552 G→C 7.7% A307G (GCC→GGC)  algJ ← alginate O‑acetylase AlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,460,5520GC 7.7% 243.7 / 10.3 78A307G (GCC→GGC) algJalginate O‑acetylase AlgJ
Reads supporting (aligned to +/- strand):  ref base G (46/26);  new base C (2/4);  total (48/30)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTT  >  NC_002947/1460405‑1460700
                                                                                                                                                    |                                                                                                                                                    
tttttCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGcc                                                                                                                                                    >  1:301497/1‑150 (MQ=255)
       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                                                                                                                                           >  1:1051267/1‑127 (MQ=255)
                                ccACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCAGCAGGCCCTGCGCAgc                                                                                                                                           <  2:1051267/127‑1 (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:1579569/150‑1 (MQ=255)
                                         cacTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCt                                                                                                           >  1:1082930/1‑150 (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:368383/1‑150 (MQ=255)
                                                                                    ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                >  2:1923282/1‑150 (MQ=255)
                                                                                    ttGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCg                                                                >  2:764658/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                                                      >  2:761188/1‑150 (MQ=255)
                                                                                              gcagcaGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGt                                                      >  1:1731058/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)
                                                                                                                                                    |                                                                                                                                                    
TTTTTCATTGGCAGATAACGTTCGGGGAATTCCCACACCACCACTTGTGGTGCGGCGTTCTTGAAGGCATCGCTTTGCAGGTACTTGAGCATCGGCAGCAGCGGGCCATGGCCGTCTTCGGCGTAGTTGGCGACGTCGCTGCGCAGGGCCTGCTGCAGCGCGCCAAGGAAGTTCCAGTGCGGGTTGGCGCTGTAGCTGGTGCCTACCAGTGCAACCGGGATCTGCTTGTCGGCGAACAACGCGTCGCCGGCGTCACCCTCGGCATCGACCGGCCGGGTGGTGCGTTTCTGCAGGTT  >  NC_002947/1460405‑1460700

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