Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,025 T→C 7.7% intergenic (‑218/+140) mscS ← / ← fbaA small conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,070,0250TC 7.7% 304.0 / 12.0 103intergenic (‑218/+140)mscS/fbaAsmall conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base T (47/48);  new base C (3/5);  total (50/53)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 104 total aligned reads displayed.

GGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTACAGAAC  >  NC_000913/3069879‑3070173
                                                                                                                                                  |                                                                                                                                                    
aggTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                                   >  1:421301/2‑150 (MQ=255)
   tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                >  1:1132089/1‑150 (MQ=255)
   tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                <  1:370019/150‑1 (MQ=255)
   tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                                >  1:121117/1‑150 (MQ=255)
   tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                <  2:325992/150‑1 (MQ=255)
   tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTCCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGAGCTTTGTCATGCCGGATGCGGCCTGAACGCCTTATCCgg                                                                                                                                                >  2:27673/1‑150 (MQ=255)
    aaaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                               <  1:945155/150‑1 (MQ=255)
    aaaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                               <  2:824862/150‑1 (MQ=255)
     aaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                              >  1:820148/1‑149 (MQ=255)
     aaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                              >  1:826114/1‑149 (MQ=255)
      aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                             <  2:820148/150‑3 (MQ=38)
      aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                             <  2:826114/150‑3 (MQ=38)
       cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1234873/148‑1 (MQ=255)
       cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              >  2:1441715/1‑148 (MQ=255)
       cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              >  2:1234873/1‑148 (MQ=255)
       cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                              <  1:1441715/148‑1 (MQ=255)
       cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatg                                                                                                                                            <  2:323002/150‑4 (MQ=39)
         aTAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                          <  1:367543/150‑1 (MQ=255)
         aTAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcatgac                                                                                                                                          <  1:1248115/150‑6 (MQ=38)
           aaaGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaa                                                                                                                                        >  2:820056/1‑150 (MQ=255)
              gCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAAc                                                                                                                                     >  2:844320/1‑150 (MQ=255)
                 tAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTg                                                                                                                                  <  2:1195790/150‑1 (MQ=255)
                 tAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTg                                                                                                                                  <  2:1189850/150‑1 (MQ=255)
                 tAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTg                                                                                                                                  >  1:1313862/1‑150 (MQ=255)
                   aTTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCa                                                                                                                                >  1:635298/1‑150 (MQ=255)
                   aTTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCa                                                                                                                                <  2:1066373/150‑1 (MQ=255)
                      ttGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAAt                                                                                                                             <  1:1156962/150‑1 (MQ=255)
                          cAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCaa                                                                                                                         >  1:1295954/1‑150 (MQ=255)
                          cAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCaa                                                                                                                         >  1:179849/1‑150 (MQ=255)
                                gATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATAtt                                                                                                                   <  2:1026004/150‑1 (MQ=255)
                                 aTCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGAGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGTCCTACAAAAACCTGCAAATTCAATATATTg                                                                                                                  >  2:130183/1‑150 (MQ=255)
                                  tCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGc                                                                                                                 >  2:1016876/1‑150 (MQ=255)
                                   cacaGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCa                                                                                                                >  2:1608519/1‑150 (MQ=255)
                                   cacaGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCa                                                                                                                >  2:1152668/1‑150 (MQ=255)
                                      aGATTTCTCATCTGAATCACCTTGCTTTTGTCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGa                                                                                                             <  2:1485729/150‑1 (MQ=255)
                                        aTTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATa                                                                                                           >  1:1318483/1‑150 (MQ=255)
                                          ttCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACg                                                                                                         <  2:739678/150‑1 (MQ=255)
                                            ctcATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACgtg                                                                                                       <  1:594903/150‑1 (MQ=255)
                                             tcATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACgtgt                                                                                                      <  2:1384212/150‑1 (MQ=255)
                                              cATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTa                                                                                                     >  2:549358/1‑150 (MQ=255)
                                              cATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTa                                                                                                     <  1:1612869/150‑1 (MQ=255)
                                              cATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTa                                                                                                     <  1:688302/150‑1 (MQ=255)
                                                tCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAgg                                                                                                   >  1:1665306/1‑150 (MQ=255)
                                                tCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAgg                                                                                                   <  2:1423158/150‑1 (MQ=255)
                                                tCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTACAGGATACGTGTAgg                                                                                                   <  1:43134/150‑1 (MQ=255)
                                                  tGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGcc                                                                                                 >  2:533057/1‑150 (MQ=255)
                                                  tGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGcc                                                                                                 >  1:1197299/1‑150 (MQ=255)
                                                         ccTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAg                                                                                          >  2:656567/1‑150 (MQ=255)
                                                           ttGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCg                                                                                        <  1:533057/150‑1 (MQ=255)
                                                              cTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAg                                                                                     >  2:1222623/1‑150 (MQ=255)
                                                               ttttGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAgc                                                                                    <  1:796333/150‑1 (MQ=255)
                                                               ttttGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAgc                                                                                    <  1:1581172/150‑1 (MQ=255)
                                                                tttGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATTTATTGCAGGATACGTGTAGGCCTGATAAGCGTAgcg                                                                                   >  1:676380/1‑150 (MQ=255)
                                                                    gCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATc                                                                               >  1:612481/1‑150 (MQ=255)
                                                                     cGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCa                                                                              >  1:1425826/1‑150 (MQ=255)
                                                                        aaaaaTCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                           <  1:1152668/150‑1 (MQ=255)
                                                                               cccGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTg                                                                    <  1:844320/150‑1 (MQ=255)
                                                                               cccGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTg                                                                    >  2:90367/1‑150 (MQ=255)
                                                                                 cGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCt                                                                  <  1:820056/150‑1 (MQ=255)
                                                                                   cAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCttt                                                                >  2:811500/1‑150 (MQ=255)
                                                                                    agagCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCtttt                                                               >  2:90499/1‑150 (MQ=255)
                                                                                              tCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGaaaaaa                                                     >  1:879549/1‑150 (MQ=255)
                                                                                               cagacagaTATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGaaaaaaa                                                    >  1:1594798/1‑150 (MQ=255)
                                                                                                     gaTATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCg                                              >  2:745730/1‑150 (MQ=255)
                                                                                                           aaTGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGc                                        >  2:942456/1‑150 (MQ=255)
                                                                                                             tgtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGTCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCgg                                      <  2:817494/150‑1 (MQ=255)
                                                                                                             tgtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCgg                                      <  2:1318483/150‑1 (MQ=255)
                                                                                                             tgtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCgg                                      <  1:1639579/150‑1 (MQ=255)
                                                                                                              gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                                     >  2:814727/1‑150 (MQ=255)
                                                                                                              gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                                     >  1:177886/1‑150 (MQ=255)
                                                                                                              gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                                     >  2:322907/1‑150 (MQ=255)
                                                                                                               tgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGc                                    >  2:1043307/1‑150 (MQ=255)
                                                                                                                 cgcTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCt                                  >  1:463416/1‑150 (MQ=255)
                                                                                                                  gcTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCtt                                 <  1:1320667/150‑1 (MQ=255)
                                                                                                                  gcTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCT‑‑AAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTga                               <  1:108046/150‑1 (MQ=255)
                                                                                                                  gcTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCT‑‑AAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTga                               <  1:114051/150‑1 (MQ=255)
                                                                                                                   cTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCATGATACGTGTAGTCCTGATAAGCGTAGCGCATCAGTCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTg                                <  2:1295954/150‑1 (MQ=255)
                                                                                                                     ttGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCTTAGATCATCAGGCAATTTTGCTTTTGGCTAAAAAAAGACCCGCAGATCTGGACTTgag                              <  2:179849/150‑1 (MQ=21)
                                                                                                                     ttGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTgag                              <  2:1665306/150‑1 (MQ=255)
                                                                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTgaga                             >  1:1500585/7‑150 (MQ=25)
                                                                                                                      tGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCTCTGCGGGTCtttttt                             <  1:1016876/150‑5 (MQ=18)
                                                                                                                      tGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTgaga                             <  1:1566419/150‑1 (MQ=255)
                                                                                                                       gTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGGCCCGCTCTGCGGGTCttttttt                            <  1:656567/150‑6 (MQ=18)
                                                                                                                       gTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGAt                            <  2:612481/150‑1 (MQ=255)
                                                                                                                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTATTCCTGATAAGCGTAGCGCATCAGTCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGc                        <  2:1611635/150‑1 (MQ=255)
                                                                                                                           tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGc                        <  2:421301/150‑1 (MQ=255)
                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa                       <  2:177886/150‑1 (MQ=255)
                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa                       >  1:740355/1‑150 (MQ=255)
                                                                                                                            gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa                       >  2:1527626/1‑150 (MQ=255)
                                                                                                                                    cggcgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGaa               <  2:740355/150‑1 (MQ=255)
                                                                                                                                      gcgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAata             >  1:1423888/1‑150 (MQ=255)
                                                                                                                                       cgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAatat            >  1:1614932/1‑150 (MQ=255)
                                                                                                                                       cgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAatat            >  1:1075887/1‑150 (MQ=255)
                                                                                                                                           aaCGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTa        <  1:1162735/150‑1 (MQ=255)
                                                                                                                                            aCGCCTTATCCGTCCTACAAAAACCTTCAAATTAAATATATTGCATGATACGTGTAGGACTGATAAGCGTAGCGAATCAGGCAATTTTGCTTTTTGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTAc       <  2:832089/150‑1 (MQ=38)
                                                                                                                                            aCGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGTATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTAc       <  2:832088/150‑1 (MQ=255)
                                                                                                                                            aCGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTAc       >  2:743430/1‑150 (MQ=255)
                                                                                                                                               ccTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTACAGa    >  1:503122/1‑150 (MQ=255)
                                                                                                                                               ccTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTACAGa    >  1:517676/1‑150 (MQ=255)
                                                                                                                                                 ttATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCTGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATTTCTTACAGAAc  >  1:555696/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATATCTTACAGAAC  >  NC_000913/3069879‑3070173

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