Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,011 C→A 8.0% intergenic (‑204/+154) 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,0110CA 8.0% 326.5 / 12.8 100intergenic (‑204/+154)mscS/fbaAsmall conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base C (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 = 9.57e-01

Only 100 of 108 total aligned reads displayed.

GACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAA  >  NC_000913/3069862‑3070160
                                                                                                                                                     |                                                                                                                                                     
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:121117/1‑150 (MQ=255)
                    tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                   >  1:1132089/1‑150 (MQ=255)
                    tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                   <  1:370019/150‑1 (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:826114/150‑3 (MQ=38)
                       aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                <  2:820148/150‑3 (MQ=38)
                        cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:1441715/1‑148 (MQ=255)
                        cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 <  1:1441715/148‑1 (MQ=255)
                        cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 >  2:1234873/1‑148 (MQ=255)
                        cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                 <  1:1234873/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)
                                                           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: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)
                                                                                                                                                     |                                                                                                                                                     
GACGGGCAGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAA  >  NC_000913/3069862‑3070160

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: