Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,299:1 +G 5.7% coding (875/897 nt) ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,2991.G 5.7% 784.2 / 41.8 176coding (875/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (83/83);  new base G (3/7);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 178 total aligned reads displayed.

CATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAC  >  NC_000913/3066152‑3066447
                                                                                                                                                    |                                                                                                                                                     
caTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATc                                                                                                                                                      <  4:1303795/149‑1 (MQ=255)
    cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                  >  2:1436874/1‑149 (MQ=255)
    cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                  <  3:1020850/149‑1 (MQ=255)
    cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                  <  3:40889/149‑1 (MQ=255)
    cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                                                                                                  <  1:268271/149‑1 (MQ=255)
    cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACCCGTCCCGTCACGATCTTCa                                                                                                                                                  <  3:1297999/149‑1 (MQ=255)
      aTGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGc                                                                                                                                                <  2:736940/149‑1 (MQ=255)
           tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAatat                                                                                                                                          >  1:1081546/1‑149 (MQ=255)
           tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAatat                                                                                                                                          >  2:160010/1‑149 (MQ=255)
           tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAatat                                                                                                                                          >  4:1196712/1‑149 (MQ=255)
           tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCATCACGATCTTCAGCAatat                                                                                                                                          >  4:1069868/1‑149 (MQ=255)
                 ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                    <  2:1323369/149‑1 (MQ=255)
                 ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                    <  2:1081546/149‑1 (MQ=255)
                 ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                    <  3:379469/149‑1 (MQ=255)
                    tCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAga                                                                                                                                  <  3:755930/149‑2 (MQ=255)
                       gCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCg                                                                                                                              >  3:64836/1‑149 (MQ=255)
                         aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAgatc                                                                                                                                <  1:1007809/146‑4 (MQ=255)
                            gCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCg                                                                                                                         >  2:867086/1‑149 (MQ=255)
                             cTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                                                        <  1:1511690/149‑1 (MQ=38)
                             cTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                                                        <  2:1154675/149‑1 (MQ=38)
                                tCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtat                                                                                                                     <  1:1190784/149‑1 (MQ=255)
                                  tGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                                                                                                   >  1:1110094/1‑149 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAata                                                                                                                                            >  3:892985/1‑124 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAata                                                                                                                                            <  4:892985/124‑1 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                                                                                                  <  1:1380194/149‑1 (MQ=255)
                                   ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                                                                                                  >  4:611089/1‑149 (MQ=255)
                                    gAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACg                                                                                                                 >  1:378427/1‑149 (MQ=255)
                                     aaaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                                <  1:1401954/149‑1 (MQ=255)
                                      aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgat                                                                                                               <  3:110655/149‑1 (MQ=255)
                                          gcgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAtt                                                                                                           <  1:229064/149‑1 (MQ=255)
                                               tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                                                      <  4:1533628/149‑1 (MQ=255)
                                                ggCGGACATATTCACCACGATCTTAGCCGCGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                     <  4:775635/149‑1 (MQ=255)
                                                ggCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                     <  1:1542580/149‑1 (MQ=255)
                                                   ggACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                                                  >  2:1192372/1‑149 (MQ=255)
                                                     aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                <  2:1186814/149‑1 (MQ=255)
                                                       atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                                              >  1:1349549/1‑149 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                         <  1:1030374/135‑1 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                         >  2:1030374/1‑135 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                           >  3:1337143/1‑149 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                           >  3:1089099/1‑149 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                           <  3:177226/149‑1 (MQ=255)
                                                          ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                           <  1:65756/149‑1 (MQ=255)
                                                               cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                             <  2:245256/142‑1 (MQ=255)
                                                               cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                      >  4:1096386/1‑149 (MQ=255)
                                                               cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                      >  1:144534/1‑149 (MQ=255)
                                                               cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                      >  4:691285/1‑149 (MQ=255)
                                                                acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                                                     >  3:1525473/1‑149 (MQ=255)
                                                                      ttAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                >  2:861752/1‑132 (MQ=255)
                                                                         gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                                                            <  1:911088/149‑1 (MQ=255)
                                                                           cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                                                            <  4:1159465/147‑1 (MQ=255)
                                                                            aCGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                         >  3:139292/1‑149 (MQ=255)
                                                                             cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                                                        >  1:1298606/1‑149 (MQ=255)
                                                                                cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                                                     <  1:867086/149‑1 (MQ=255)
                                                                                 ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                         <  4:1563144/144‑1 (MQ=255)
                                                                                 ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                                                    <  3:81833/149‑1 (MQ=255)
                                                                                 ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                                                    <  2:1249694/149‑1 (MQ=255)
                                                                                     aCTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                                                <  2:219936/149‑1 (MQ=255)
                                                                                     aCTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                                                <  4:71691/149‑1 (MQ=255)
                                                                                       ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCATGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                 <  1:740054/114‑1 (MQ=255)
                                                                                       ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                              >  1:115616/1‑149 (MQ=255)
                                                                                         acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                         <  2:1496657/136‑1 (MQ=255)
                                                                                         acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                         >  1:1496657/1‑136 (MQ=255)
                                                                                         acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                            <  4:362636/149‑1 (MQ=255)
                                                                                         acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                            <  4:139292/149‑1 (MQ=255)
                                                                                         acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                            <  3:1569396/149‑1 (MQ=255)
                                                                                          cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                                                           <  3:691285/149‑1 (MQ=255)
                                                                                           acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATGATGATGATTTCTGa                                                                                                      <  3:117563/105‑1 (MQ=255)
                                                                                           acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                          >  1:1534081/1‑149 (MQ=255)
                                                                                                 atttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                                                    <  2:732700/149‑1 (MQ=255)
                                                                                                  tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                                                   <  2:1349549/149‑1 (MQ=255)
                                                                                                     aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTATGCGGTGCATCAAAATAtt                                                >  2:553244/1‑149 (MQ=255)
                                                                                                      attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                                               <  2:150646/149‑1 (MQ=255)
                                                                                                      attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                                               <  3:132022/149‑1 (MQ=255)
                                                                                                        tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                                             <  3:1241843/149‑1 (MQ=255)
                                                                                                          cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                           <  4:1453422/149‑1 (MQ=255)
                                                                                                          cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                           <  3:1097171/149‑1 (MQ=255)
                                                                                                           aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                          >  4:2681/1‑149 (MQ=255)
                                                                                                           aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                          >  2:946293/1‑149 (MQ=255)
                                                                                                           aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                          >  2:449828/1‑149 (MQ=255)
                                                                                                            aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGt                                         >  1:54337/1‑149 (MQ=255)
                                                                                                                 cTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg                                    <  1:160010/149‑1 (MQ=255)
                                                                                                                   cGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGa                                  <  2:749466/149‑1 (MQ=255)
                                                                                                                       gCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCa                              >  2:361873/1‑149 (MQ=255)
                                                                                                                        cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                             <  1:946293/149‑1 (MQ=255)
                                                                                                                        cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                             >  2:215553/1‑149 (MQ=255)
                                                                                                                        cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                             <  2:1535228/149‑1 (MQ=255)
                                                                                                                             ttAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCCGTc                        <  4:1337143/149‑1 (MQ=255)
                                                                                                                             ttAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                        <  1:856664/149‑1 (MQ=255)
                                                                                                                             ttAATGACTCGTCCCGTCACGATCTTCAGCAAGATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                        <  2:815909/149‑1 (MQ=255)
                                                                                                                                 tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCa                    >  4:854214/1‑149 (MQ=255)
                                                                                                                                  gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAg                   <  3:121240/149‑1 (MQ=255)
                                                                                                                                   aCTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAgg                  >  4:790097/1‑149 (MQ=255)
                                                                                                                                       gTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCaaa              >  3:719900/1‑149 (MQ=255)
                                                                                                                                           cGTCAAGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGc          >  1:457370/1‑149 (MQ=255)
                                                                                                                                            gccACGATGCTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg         >  3:1474261/3‑149 (MQ=38)
                                                                                                                                             tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg        <  1:361873/149‑1 (MQ=255)
                                                                                                                                             tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg        <  3:694920/149‑1 (MQ=255)
                                                                                                                                              cACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCggg       >  2:1561537/1‑149 (MQ=255)
                                                                                                                                               aCGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGa      <  2:280383/149‑1 (MQ=255)
                                                                                                                                                   tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATCTTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc  >  1:77550/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAC  >  NC_000913/3066152‑3066447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: