Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,017 G→T 7.7% intergenic (‑210/+148) 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,0170GT 7.7% 343.7 / 11.5 104intergenic (‑210/+148)mscS/fbaAsmall conductance mechanosensitive channel MscS/fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base G (48/48);  new base T (3/5);  total (51/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.03e-01

Only 100 of 104 total aligned reads displayed.

AGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATAT  >  NC_000913/3069869‑3070163
                                                                                                                                                    |                                                                                                                                                  
aGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGa                                                                                                                                                   <  1:566152/150‑1 (MQ=255)
 gTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGaa                                                                                                                                                  >  1:1163124/1‑150 (MQ=255)
   cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGTCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATAAATGTGCGCTTTTTCATGCCGGATGCGGCGTGAAcc                                                                                                                                                >  1:832089/1‑149 (MQ=255)
   cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACg                                                                                                                                                <  2:15598/150‑1 (MQ=255)
   cTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACg                                                                                                                                                >  1:832088/1‑150 (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:1132089/1‑150 (MQ=255)
             tAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                      >  1:121117/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:826114/1‑149 (MQ=255)
               aaCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGca                                                                                                                                    >  1:820148/1‑149 (MQ=255)
                aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                   <  2:820148/150‑3 (MQ=38)
                aCAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGGCAGATATCAATGTGCGCTTTGTCATGCCGGATAAGGCGTTCACGCCGCATCCGGcat                                                                                                                                   <  2:826114/150‑3 (MQ=38)
                 cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  1:1441715/148‑1 (MQ=255)
                 cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  1:1234873/148‑1 (MQ=255)
                 cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    >  2:1441715/1‑148 (MQ=255)
                 cAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    >  2:1234873/1‑148 (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)
                                           aTCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGAGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGTCCTACAAAAACCTGCAAATTCAATATATTg                                                                                                        >  2:130183/1‑150 (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)
                                                  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                                                                                         <  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:1581172/150‑1 (MQ=255)
                                                                         ttttGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAgc                                                                          <  1:796333/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                                                          >  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:322907/1‑150 (MQ=255)
                                                                                                                        gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                           >  2:814727/1‑150 (MQ=255)
                                                                                                                        gtgCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCggg                           >  1:177886/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)
                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGaa     <  2:740355/150‑1 (MQ=255)
                                                                                                                                                gcgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAata   >  1:1423888/1‑150 (MQ=255)
                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAatat  >  1:1620628/1‑150 (MQ=255)
                                                                                                                                                 cgTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAatat  >  1:1075887/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTCTGACAAGGGTAAACAATAAAGCATAATTTTGCCAAATAGATCACAGATTTCTCATCTGAATCACCTTGCTTTTGGCGAAAAAATCCCCGCAGAGCGGGGCTCAGACAGATATCAATGTGCGCTTTGTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCTTTTGGCGAAAAAAAGACCCGCAGAGCGGGCCTTGAGATAAGCAGAAAGGAATAT  >  NC_000913/3069869‑3070163

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: