Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,026 T→C 6.1% intergenic (+72/‑125) PP_3009 → / → PP_3010 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,0260TC 6.1% 192.2 / 2.4 65intergenic (+72/‑125)PP_3009/PP_3010hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCTGGCA  >  NC_002947/3401889‑3402174
                                                                                                                                         |                                                                                                                                                    
aGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGcaaacaaa                                                                                                                                          >  1:1714076/1‑150 (MQ=255)
 gCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                                                         <  2:1872221/150‑1 (MQ=255)
       ttGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATg                                                                                                                                   >  2:1940152/1‑150 (MQ=255)
          ctgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGccc                                                                                                                                <  1:1940152/150‑1 (MQ=255)
           tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                                               >  2:1116278/1‑150 (MQ=255)
           tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                                               >  2:1904038/1‑150 (MQ=255)
           tgtgctgGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCt                                                                                                                               >  1:334271/1‑150 (MQ=255)
               ctgGTACTCGATCTGGCGGCGACTTATGTGGTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgttg                                                                                                                           >  1:1810001/1‑150 (MQ=255)
                   tACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGt                                                                                                                       <  1:1553529/150‑1 (MQ=255)
                          tCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAg                                                                                                                >  1:95004/1‑150 (MQ=255)
                                ggcgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACAt                                                                                                          >  1:1899835/1‑150 (MQ=255)
                                gacgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACAt                                                                                                          >  1:125284/3‑150 (MQ=255)
                                 gcgACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATg                                                                                                         <  2:334271/150‑1 (MQ=255)
                                     cTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                                                     <  2:1338461/150‑1 (MQ=255)
                                      ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                                                    <  1:1711036/150‑1 (MQ=255)
                                      ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                                                    >  1:297177/1‑150 (MQ=255)
                                      ttATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAg                                                                                                    >  2:1811033/1‑150 (MQ=255)
                                         tgtgTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgtt                                                                                                                            >  2:1594708/1‑123 (MQ=255)
                                         tgtgTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCtgtt                                                                                                                            <  1:1594708/123‑1 (MQ=255)
                                           tgtTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACaa                                                                                               >  1:791847/1‑150 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCtatata                                                                                                                                                   <  1:916230/92‑1 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCtatata                                                                                                                                                   >  2:916230/1‑92 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                                                         <  2:1903040/102‑1 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCaaacaaac                                                                                                                                         >  1:1903040/1‑102 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                                                     <  1:1228159/138‑1 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCa                                                                                                     >  2:1228159/1‑138 (MQ=255)
                                                 tcAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCa                                                                                         >  2:1032348/1‑150 (MQ=255)
                                                    aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTCGt                                                                                                                       <  1:1937478/117‑1 (MQ=255)
                                                    aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGt                                                                                                                       >  2:1937478/1‑117 (MQ=255)
                                                    aCGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaa                                                                                      >  1:1488950/1‑150 (MQ=255)
                                                     cGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaa                                                                                     <  1:1811033/150‑1 (MQ=255)
                                                       ccATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAg                                                                                   <  2:1488950/150‑1 (MQ=255)
                                                       ccATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAg                                                                                   <  2:125284/150‑1 (MQ=255)
                                                        cATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAgg                                                                                  >  2:1322237/1‑150 (MQ=255)
                                                         aTCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAggg                                                                                 >  2:1896442/1‑150 (MQ=255)
                                                         aTCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAggg                                                                                 <  1:1904038/150‑1 (MQ=255)
                                                            cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                                              >  1:1252357/1‑150 (MQ=255)
                                                            cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                                              >  1:1252714/1‑150 (MQ=255)
                                                            cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaaa                                                                                    <  1:1830572/144‑1 (MQ=255)
                                                            cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCaaaaaa                                                                                    >  2:1830572/1‑144 (MQ=255)
                                                            cccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCa                                                                              <  2:1714076/150‑1 (MQ=255)
                                                             ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCa                                                                                                                 <  2:1458017/114‑1 (MQ=255)
                                                             ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCa                                                                                                                 >  1:1458017/1‑114 (MQ=255)
                                                             ccTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAg                                                                             >  2:412966/1‑150 (MQ=255)
                                                                   cccGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGgctgct                                                                       >  2:216400/1‑150 (MQ=255)
                                                                     cGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTg                                                                     <  1:1322237/150‑1 (MQ=255)
                                                                        tAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCa                                                                  >  2:1040955/1‑150 (MQ=255)
                                                                                 ggCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTccc                                                         >  1:1745878/1‑150 (MQ=255)
                                                                                      tACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGttt                                                    <  1:412966/150‑1 (MQ=255)
                                                                                       aCTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                                                   <  2:1252357/150‑1 (MQ=255)
                                                                                       aCTCGTAATCGGGGGCAAGGCAGGTCGTGAGGGAAACGCCCCGAACCCTACATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                                                   <  2:1252714/150‑1 (MQ=255)
                                                                                       aCTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                                                   <  1:216400/150‑1 (MQ=255)
                                                                                       aCTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAAAAAAAAAAAAAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGtttt                                                   <  2:1810001/150‑1 (MQ=255)
                                                                                             aaTCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATgg                                             <  1:1040955/150‑1 (MQ=255)
                                                                                                            aGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTgg                                                                    >  1:663805/1‑112 (MQ=255)
                                                                                                            aGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTgg                                                                    <  2:663805/112‑1 (MQ=255)
                                                                                                             gggCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCgtgt                             >  2:1417036/1‑150 (MQ=255)
                                                                                                                  tGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGt                        <  1:1116278/150‑1 (MQ=255)
                                                                                                                         aaCGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTaa                                                                                                             >  2:1395030/1‑58 (MQ=255)
                                                                                                                         aaCGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTaa                                                                                                             <  1:1395030/58‑1 (MQ=255)
                                                                                                                            gcgccGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTc               >  1:1939112/4‑149 (MQ=255)
                                                                                                                            gCCCCGAACCCTCTATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCa              <  1:1363762/150‑1 (MQ=255)
                                                                                                                            gCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTc               <  2:1939112/149‑1 (MQ=255)
                                                                                                                             ccccGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCaa             <  1:1896442/150‑1 (MQ=255)
                                                                                                                              cccaaacccaaaATAGCAAACAAACCAGATGCCCTGTTGAAGTTGAAAAGTAAAATGGGCAAAAAAAAAACAAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCAAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAg            <  2:791847/150‑1 (MQ=255)
                                                                                                                              cccGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAg            >  1:1865356/1‑150 (MQ=255)
                                                                                                                                    cccTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCt      >  2:444356/1‑150 (MQ=255)
                                                                                                                                    cccTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGAAGAACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCt      <  2:95004/150‑1 (MQ=255)
                                                                                                                                       tatataGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCTg     <  1:978383/148‑1 (MQ=255)
                                                                                                                                       tatataGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCTg     >  2:978383/1‑148 (MQ=255)
                                                                                                                                        atataGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCTGGCa  >  1:64767/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
AGCGTGGTTGCTGTGCTGGTACTCGATCTGGCGGCGACTTATGTGTTTCTCAACGCCATCCCCTGACCCCGGTAATGGGTAGGCCTTACTCGTAATCGGGGGCAAGGCAGGGCGTGAGGGAAACGCCCCGAACCCTATATAGCAAACAAACCTGATGCCCTGTTGTAGTTGAGCAGTAACATGGGCAGGACAAAAACCAAAAAAGGGGCAGGCTGCTTGGCAACATTTCCCAGTTTTCGATGGGGCTGCAGGAGCGTGTCGCGTGAAAAACTCAAGATCGCTGGCA  >  NC_002947/3401889‑3402174

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: