Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,010 G→A 8.0% intergenic (‑203/+155) 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,0100GA 8.0% 298.4 / 11.7 100intergenic (‑203/+155)mscS/fbaAsmall conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base G (44/48);  new base A (3/5);  total (47/53)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 108 total aligned reads displayed.

TGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCA  >  NC_000913/3069861‑3070152
                                                                                                                                                     |                                                                                                                                              
tGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATg                                                                                                                                                <  2:716374/150‑1 (MQ=255)
 gACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgc                                                                                                                                               >  1:72609/1‑150 (MQ=255)
   cGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgcgg                                                                                                                                             <  1:565978/150‑1 (MQ=255)
   cGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgcgg                                                                                                                                             >  1:1591541/1‑150 (MQ=255)
    gggCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgcggc                                                                                                                                            <  2:123498/150‑1 (MQ=255)
     ggCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgcggcg                                                                                                                                           >  2:1320667/1‑150 (MQ=255)
     ggCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATgcggcg                                                                                                                                           <  2:72609/150‑1 (MQ=255)
      gCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGt                                                                                                                                          <  1:963833/150‑1 (MQ=255)
        aGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGa                                                                                                                                        <  1:566152/150‑1 (MQ=255)
         gTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGaa                                                                                                                                       >  1:1163124/1‑150 (MQ=255)
           cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGTCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATAAATGTGCGCTTTTTCATGCCGGATGCGGCGTGAAcc                                                                                                                                     >  1:832089/1‑149 (MQ=255)
           cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACg                                                                                                                                     >  1:832088/1‑150 (MQ=255)
           cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACg                                                                                                                                     <  2:15598/150‑1 (MQ=255)
             gACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                   >  1:1611635/1‑150 (MQ=255)
              aCAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                                  <  2:640499/150‑1 (MQ=255)
                  aggTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                              >  1:421301/2‑150 (MQ=255)
                     tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                           <  1:370019/150‑1 (MQ=255)
                     tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                           >  1:1132089/1‑150 (MQ=255)
                     tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                           >  1:121117/1‑150 (MQ=255)
                     tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                           <  2:325992/150‑1 (MQ=255)
                     tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTCCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGAGCTTTGTCATGCCGGATGCGGCCTGAACGCCTTATCCgg                                                                                                                           >  2:27673/1‑150 (MQ=255)
                      aaaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                          <  2:824862/150‑1 (MQ=255)
                      aaaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                          <  1:945155/150‑1 (MQ=255)
                       aaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                         >  1:820148/1‑149 (MQ=255)
                       aaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                         >  1:826114/1‑149 (MQ=255)
                        aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                        <  2:826114/150‑3 (MQ=38)
                        aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                        <  2:820148/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)
                                 cATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAAcc                                                                                                               <  1:864656/150‑1 (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:179849/1‑150 (MQ=255)
                                            cAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCaa                                                                                                    >  1:1295954/1‑150 (MQ=255)
                                                  gATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATAtt                                                                                              <  2:1026004/150‑1 (MQ=255)
                                                    tCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGc                                                                                            >  2:1016876/1‑150 (MQ=255)
                                                     cacaGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCa                                                                                           >  2:1152668/1‑150 (MQ=255)
                                                     cacaGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCa                                                                                           >  2:1608519/1‑150 (MQ=255)
                                                        aGATTTCTCATCTGAATCACCTTGCTTTTGTCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGa                                                                                        <  2:1485729/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:688302/150‑1 (MQ=255)
                                                                cATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTa                                                                                <  1:1612869/150‑1 (MQ=255)
                                                                  tCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAgg                                                                              <  2:1423158/150‑1 (MQ=255)
                                                                  tCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAgg                                                                              >  1:1665306/1‑150 (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)
                                                                                      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                 <  1:1639579/150‑1 (MQ=255)
                                                                                                                               tgtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCgg                 <  2:1318483/150‑1 (MQ=255)
                                                                                                                                gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                >  2:322907/1‑150 (MQ=255)
                                                                                                                                gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                >  2:814727/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)
                                                                                                                                       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)
                                                                                                                                            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAg    <  2:310125/150‑1 (MQ=255)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTATTCCTGATAAGCGTAGCGCATCAGTCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGc   <  2:1611635/150‑1 (MQ=255)
                                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGc   <  2:421301/150‑1 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa  >  1:740355/1‑150 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa  <  2:177886/150‑1 (MQ=255)
                                                                                                                                              gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATATTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCa  >  2:1527626/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCA  >  NC_000913/3069861‑3070152

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: