Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,319 Δ1 bp 5.6% coding (855/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,3190C. 5.6% 448.1 / 19.1 124coding (855/897 nt)ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (51/66);  new base . (4/3);  total (55/69)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 125 total aligned reads displayed.

ATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGC  >  NC_000913/3066171‑3066466
                                                                                                                                                      |                                                                                                                                                   
aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                                     <  3:295769/149‑1 (MQ=255)
 tCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                                                                                                                    >  3:694993/1‑149 (MQ=255)
   tGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                                  >  1:227877/1‑149 (MQ=255)
    gCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCg                                                                                                                                                 >  2:571981/1‑149 (MQ=255)
     cATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATGCTTCAGCATATTGCTGAAGATCGt                                                                                                                                                <  3:407105/149‑1 (MQ=255)
      aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTg                                                                                                                                               >  3:302291/1‑149 (MQ=255)
       tCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                                                                                                                              <  2:750110/149‑1 (MQ=255)
         gCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCg                                                                                                                                            <  2:479113/149‑1 (MQ=255)
          cTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGCTCTTCATCAATATGCTGAAGCATCGTGGCg                                                                                                                                            >  2:1151066/1‑148 (MQ=255)
          cTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCg                                                                                                                                            <  1:1151066/148‑1 (MQ=255)
            aTCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                                                                                                                         <  2:862630/149‑1 (MQ=255)
               tGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                                                                                                                      <  2:227877/149‑1 (MQ=255)
                   aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgat                                                                                                                                  <  3:384844/149‑1 (MQ=255)
                      tGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                                                                                                                               <  4:373163/149‑1 (MQ=255)
                            tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                                                                         <  1:887822/149‑1 (MQ=255)
                            tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                                                                         <  3:1057511/149‑1 (MQ=255)
                            tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                                                                                                                         <  1:1050128/149‑1 (MQ=255)
                              gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                                       <  4:465227/149‑1 (MQ=255)
                              gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                                       <  4:694993/149‑1 (MQ=255)
                              gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAAAATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                                       <  2:971481/149‑1 (MQ=255)
                              gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATATAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCACCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                                       <  4:705797/149‑1 (MQ=255)
                                  aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTgac                                                                                                                                             >  4:866894/1‑121 (MQ=37)
                                  aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTgac                                                                                                                                             <  3:866894/123‑3 (MQ=37)
                                     tatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                                <  1:6198/149‑1 (MQ=255)
                                      atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGc                                                                                                                                             <  4:579523/119‑1 (MQ=255)
                                      atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGc                                                                                                                                             >  3:579523/1‑119 (MQ=255)
                                       ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCTGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                                  >  2:870733/1‑113 (MQ=255)
                                       ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                                  <  1:870733/113‑1 (MQ=255)
                                       ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                                              <  1:682503/149‑1 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGTTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                                              >  4:631208/1‑144 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGt                                                                                                                                                >  4:210617/1‑110 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                                                                                                                      >  1:297521/1‑120 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                                                                                                                      <  2:297521/120‑1 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                                                              <  3:631208/144‑1 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                         >  2:304950/1‑149 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                         >  4:367166/1‑149 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                         <  2:496843/149‑1 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                                         >  2:996296/1‑149 (MQ=255)
                                            cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCACGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                                  >  2:1091929/1‑108 (MQ=255)
                                                   ttAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACGCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCata                                                                                                  >  3:715651/1‑149 (MQ=255)
                                                      gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAc                                                                                               <  4:302291/149‑1 (MQ=255)
                                                       ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                                                                              >  1:513927/1‑149 (MQ=255)
                                                       ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                                                                              <  1:849226/149‑1 (MQ=255)
                                                       ccACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACa                                                                                              >  4:859790/1‑149 (MQ=255)
                                                          cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                                                                           >  1:1124177/1‑149 (MQ=255)
                                                            tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCGTCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                                         <  1:936681/149‑1 (MQ=255)
                                                            tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                                         <  3:431850/149‑1 (MQ=255)
                                                            tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                                         <  1:730518/149‑1 (MQ=255)
                                                            tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                                         <  3:658514/149‑1 (MQ=255)
                                                            tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                                         >  3:795832/1‑149 (MQ=255)
                                                              ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                                                                       >  1:1019419/1‑149 (MQ=255)
                                                                    ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                                                                      <  2:713012/112‑1 (MQ=255)
                                                                    ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                                                                      >  1:713012/1‑112 (MQ=255)
                                                                      acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                                                               <  1:538688/149‑1 (MQ=255)
                                                                        acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                                                   <  4:251324/95‑1 (MQ=255)
                                                                        acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                                                   >  3:251324/1‑95 (MQ=255)
                                                                           gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                                                                          >  4:1013012/1‑149 (MQ=255)
                                                                               tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAat                                                                      >  4:412139/1‑149 (MQ=255)
                                                                                  aattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                                                                   >  4:881611/1‑149 (MQ=255)
                                                                                   attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                                                 <  4:44853/134‑1 (MQ=255)
                                                                                   attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                                                 >  3:44853/1‑134 (MQ=255)
                                                                                   attGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAttt                                                                  <  1:571981/149‑1 (MQ=255)
                                                                                      gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGTGGTGCATCAAAATATTTCAg                                                               <  2:1019419/149‑1 (MQ=255)
                                                                                       cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                                              <  1:1140770/149‑1 (MQ=255)
                                                                                       cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                                              >  3:731715/1‑149 (MQ=255)
                                                                                          ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGtt                                                           <  4:795832/149‑1 (MQ=255)
                                                                                                 gATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt                                                    >  3:114833/1‑149 (MQ=255)
                                                                                                  aTGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAtt                                                   <  3:412139/149‑1 (MQ=255)
                                                                                                     cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAc                                                <  1:304950/149‑1 (MQ=255)
                                                                                                      tGGGTTAATGTCTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTCCATCAAAATATTTCAGCGTTGCCGGATTCAcc                                               <  4:731715/149‑1 (MQ=255)
                                                                                                      tGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAcc                                               >  2:1020305/1‑149 (MQ=255)
                                                                                                          ttAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                                                                       <  3:21346/89‑1 (MQ=255)
                                                                                                          ttAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                                           <  1:844431/149‑1 (MQ=255)
                                                                                                          ttAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                                           <  2:513927/149‑1 (MQ=255)
                                                                                                            aaTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTc                                           <  3:1051687/147‑1 (MQ=255)
                                                                                                             aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGcc                                        <  2:181258/149‑1 (MQ=255)
                                                                                                                  tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtat                                                                                                                                        <  2:150093/48‑1 (MQ=255)
                                                                                                                  tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtat                                                                                                                                        >  1:150093/1‑48 (MQ=255)
                                                                                                                  tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                                              >  1:1001401/1‑122 (MQ=255)
                                                                                                                  tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                                                              <  2:1001401/122‑1 (MQ=255)
                                                                                                                          tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                           <  3:859790/149‑1 (MQ=255)
                                                                                                                          tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                           <  1:1020305/149‑1 (MQ=255)
                                                                                                                          tCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                           <  2:303128/149‑1 (MQ=255)
                                                                                                                            acgatgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGa                         >  2:276991/7‑149 (MQ=255)
                                                                                                                            aCGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                           >  4:119932/1‑147 (MQ=255)
                                                                                                                            aCGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                           <  3:119932/147‑1 (MQ=255)
                                                                                                                             cgatgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGAt                        <  1:276991/144‑1 (MQ=255)
                                                                                                                              gatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATc                       >  4:11763/5‑149 (MQ=255)
                                                                                                                              gATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATc                       <  1:78178/149‑1 (MQ=255)
                                                                                                                                tgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc                     >  3:636108/3‑149 (MQ=255)
                                                                                                                                tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc                     >  1:552133/1‑149 (MQ=255)
                                                                                                                                tCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc                     <  1:784289/149‑1 (MQ=255)
                                                                                                                                   ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg                   >  4:207881/1‑149 (MQ=255)
                                                                                                                                   ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg                   >  2:695944/1‑149 (MQ=255)
                                                                                                                                   ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg                   >  1:815552/1‑149 (MQ=255)
                                                                                                                                   ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg                   >  2:717739/1‑149 (MQ=255)
                                                                                                                                    tCAGTAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGt                  <  4:715651/149‑1 (MQ=255)
                                                                                                                                      aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCa                >  4:290617/1‑149 (MQ=255)
                                                                                                                                                   aaGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCcgc    >  4:572958/1‑149 (MQ=255)
                                                                                                                                                     gCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCgcgc  >  4:584505/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAAATGTCCGCGC  >  NC_000913/3066171‑3066466

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