Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,298 A→T 5.4% H293L (CAT→CTT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2980AT 5.4% 219.2 / 2.4 84H293L (CAT→CTT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base T (3/2);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCG  >  NC_000913/3181152‑3181446
                                                                                                                                                  |                                                                                                                                                    
ctcAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcatc                                                                                                                                                    >  2:1147379/1‑149 (MQ=255)
      cAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGACCCTGCTGCCCGCTTATTTTGCGGAACATGATCATCCGCaaa                                                                                                                                              >  2:759285/1‑149 (MQ=255)
      cAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGAGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTACCGAATCACCCGATCCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCata                                                                                                                                              >  1:549504/1‑147 (MQ=255)
       aaGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAAt                                                                                                                                             <  1:1147379/149‑1 (MQ=255)
          cTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataa                                                                                                                                          <  1:1292046/149‑5 (MQ=255)
             gaggaTGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCCCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGgaaa                                                                                                                                       >  1:1284988/1‑149 (MQ=255)
                gaTGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATTTTACCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAATATGATCATCCGCAAATGGAAAAatt                                                                                                                                    >  1:365140/1‑148 (MQ=255)
                 aTGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCCCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAatat                                                                                                                                   >  1:97855/1‑149 (MQ=255)
                   gCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                                                                 <  2:1376166/149‑1 (MQ=255)
                   gCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                                                                 >  1:1275605/1‑149 (MQ=255)
                           aCGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                                                                         <  1:950362/149‑1 (MQ=255)
                           aCGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                                                                         <  1:426584/149‑1 (MQ=255)
                           aCGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGAGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                                                                         <  1:429673/149‑1 (MQ=255)
                             gctggctgGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                                                                                       >  1:225486/1‑149 (MQ=255)
                               tggctggAACCGTCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTTCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAt                                                                                                                     <  1:140058/149‑1 (MQ=255)
                               tggctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAt                                                                                                                     <  2:245414/149‑1 (MQ=255)
                               gggctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAt                                                                                                                     <  1:1413431/148‑1 (MQ=255)
                                ggctggAAGCGGCGTGGAAGAGCATTATCTCCTACAAGGCACTTCTACCGCTACTGTGGGTGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                                                                                    <  2:109525/149‑1 (MQ=255)
                                 gctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGCACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtct                                                                                                                   >  1:1212603/1‑149 (MQ=255)
                                 gctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCACGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtct                                                                                                                   >  2:489162/1‑149 (MQ=255)
                                  ctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtctt                                                                                                                  >  1:402912/1‑149 (MQ=255)
                                   tggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTATGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtctta                                                                                                                 >  1:632583/1‑148 (MQ=255)
                                     gAACCGGCGTGGAAGAGCAATATCTCCCACAAGGCACGTCTACCGCTACGGTGGGAGATGGTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTtctc                                                                                                               <  2:887402/149‑1 (MQ=255)
                                        ccGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCGACCGCTTCTGTGGGAGATGTTCCCGAATCCCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTGCTCCCg                                                                                                            >  1:397956/1‑149 (MQ=255)
                                           gCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAATAATATGTGGTTAAACCGATCTTCTCCCGTGa                                                                                                         >  1:1055871/1‑149 (MQ=255)
                                             gTGGAAGTGCATTATCTACAACATGGCACTTCTAGCTCTACTGTTGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTTCGGAAGATGATCATCCGCTAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAg                                                                                                       <  2:980406/149‑1 (MQ=255)
                                             gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAg                                                                                                       >  2:1450564/1‑149 (MQ=255)
                                               ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAAAGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                                                                     >  2:412983/1‑149 (MQ=255)
                                                gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                                                                                    >  2:366347/1‑149 (MQ=255)
                                                gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTACCCGAATCAACCGACCCTGCTGCCCGCTTATTTTGCGGCAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                                                                                    >  1:1365089/1‑149 (MQ=255)
                                                gAAGAGCAATATCTCCATCTAGGCTCGTCAACCGCTACTGTGGGAGTAGTTCCCGAATCACCAGAACCTGCTTCGCGCTTTATTTGCGGAATATGTTCATCCGCATATGGAAATATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                                                                                    <  2:433377/149‑1 (MQ=255)
                                                     gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAAc                                                                                               <  2:1239887/149‑1 (MQ=255)
                                                     gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAAc                                                                                               >  2:1430378/1‑149 (MQ=255)
                                                       aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCTTTAAGGCGCAAACgt                                                                                             >  1:447466/1‑149 (MQ=255)
                                                       aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgt                                                                                             <  1:563617/149‑1 (MQ=255)
                                                         tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                                                                           >  1:481784/1‑149 (MQ=255)
                                                          aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAGCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                                                                          <  2:1203076/149‑1 (MQ=255)
                                                           tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                                                                         >  1:1080725/1‑149 (MQ=255)
                                                              ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCCCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAta                                                                                      >  1:651926/1‑148 (MQ=255)
                                                              ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATc                                                                                      <  2:34955/149‑1 (MQ=255)
                                                               caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCa                                                                                     >  1:1210417/1‑149 (MQ=255)
                                                               caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCa                                                                                     >  2:1305330/1‑149 (MQ=255)
                                                                aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAt                                                                                    <  1:1450564/149‑1 (MQ=255)
                                                                 ccaaGCCACCTCTACCGCTACTGGGGGTGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGTAAGATGATCATCCGCAAATGGAAAAATATGTGGTTACACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAtt                                                                                   <  2:757684/148‑1 (MQ=255)
                                                                 acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAtt                                                                                   <  2:823936/149‑1 (MQ=255)
                                                                 acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGGTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAtt                                                                                   <  2:1162739/149‑1 (MQ=255)
                                                                   aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                                                                                 <  1:489162/149‑1 (MQ=255)
                                                                    atgCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACTTGTCGATCATTgag                                                                                <  2:1002327/147‑1 (MQ=255)
                                                                       cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAAc                                                                             <  2:1303091/149‑1 (MQ=255)
                                                                         cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCCCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGAGTCGATCATTGAGAACgg                                                                           >  1:259565/1‑149 (MQ=255)
                                                                              aCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                                                                      >  1:1563518/1‑149 (MQ=255)
                                                                                    aCTGTGGGAGTTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTGCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg                                                                <  2:1020471/149‑1 (MQ=255)
                                                                                     cTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGa                                                               >  2:967070/1‑149 (MQ=255)
                                                                                      tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTCaa                                                              >  2:125478/1‑149 (MQ=255)
                                                                                       gtgGGAGATGTTCCCGAATCCCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag                                                             >  1:489996/1‑149 (MQ=255)
                                                                                       gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCGGCTTTTTTTGCGGTAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCGCGTGATGGCGCAAACGTGTCGATCATTGAGATCGGCAAAACCATTGAag                                                             <  1:415871/149‑1 (MQ=255)
                                                                                       gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag                                                             <  1:1451025/149‑1 (MQ=255)
                                                                                          ggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcag                                                          <  1:1275310/149‑1 (MQ=255)
                                                                                           gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcagc                                                         <  1:1430378/149‑1 (MQ=255)
                                                                                                  tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAACACCAGTGAAGCAGCGGAAggg                                                  >  1:785223/1‑148 (MQ=255)
                                                                                                        aTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTAt                                            <  2:1563518/149‑1 (MQ=255)
                                                                                                           aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAACATTGCAGCAGCTGAAGGTCCGTATGGc                                         >  1:788426/1‑149 (MQ=255)
                                                                                                           aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGGGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGc                                         <  1:751181/149‑1 (MQ=255)
                                                                                                            cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTATACCGATCTTCTCCCGTGAAGGCGCAAGCGTGTCGATCATTGAGAGCGGCAAAACCATTGAAGGAGCGGAAGGTCCGTATGGCg                                        >  2:241306/1‑149 (MQ=255)
                                                                                                            cccGAACCTGCTGCCCGCTTATTTTGAGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGAGTCGATCAGTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTCTGGCg                                        >  2:693779/1‑149 (MQ=255)
                                                                                                            cccGAACCTGCTGCCAGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAGCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCg                                        <  1:759285/149‑1 (MQ=255)
                                                                                                             ccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCga                                       >  1:1225940/1‑149 (MQ=255)
                                                                                                              cGAACCTGCTGCCCTCTTATTTTTCGGAAGATTATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCgaa                                      <  2:178438/149‑1 (MQ=255)
                                                                                                                   ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAgg                                 >  2:832308/1‑149 (MQ=255)
                                                                                                                    tgctgcCCGCTTATTTTGCGGAAGATGAGCATCCGCAAATGGAAAAATATGTGGTTACACCGATCTTCTCCCGTGAACGCGCAAACGTGGCGATCATTGAGAACGCCAAAACCATTGAAGCAGCGGAAGGGCCGTATGGCGAAGAAggg                                <  2:1466484/149‑1 (MQ=255)
                                                                                                                     gctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGATCGGCAAAACCATGGAAGCAGCGGAAGGTCCGTATGGCGAAGAAgggt                               >  1:1044746/1‑148 (MQ=255)
                                                                                                                      ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgat                              >  1:607826/1‑149 (MQ=255)
                                                                                                                        gcCCCCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgatga                            >  2:837271/1‑149 (MQ=255)
                                                                                                                         cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgatgat                           >  1:777166/1‑149 (MQ=255)
                                                                                                                         cccGCTTATTTTGAGGAAGATGAGCTTCCTTAAATGGAAAAATATGGGGTTAACCCGATCTTCTCCCGTGAAGGCGCCAACGTGTCGATCATTGAGAACGGCTAAACCATTGTAGCAGCGGTAGGTCCGATTGGCGAAGAAGGgatgat                           <  2:1252217/149‑1 (MQ=255)
                                                                                                                            gCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGt                        >  1:850482/1‑149 (MQ=255)
                                                                                                                              ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAAGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTc                      <  2:791950/149‑1 (MQ=255)
                                                                                                                                catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAg                    >  2:148845/3‑149 (MQ=255)
                                                                                                                                 ttttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGc                   <  1:366347/149‑1 (MQ=255)
                                                                                                                                 atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAATCCGATCTTCTCCCGTGATGGCGCAACCGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGc                   >  1:1000255/2‑149 (MQ=255)
                                                                                                                                 atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGc                   >  1:413340/2‑149 (MQ=255)
                                                                                                                                  tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAACCCGATCTTCTCCCGTGAAGGCGCAAAGGTGTCGATCATTGAGAACGGCAAAGCCATTGAAGCAGCGGAAGGTCGGTATGGCGAAGAAGGGATGATTGTTCAGCa                  >  1:176390/1‑149 (MQ=255)
                                                                                                                                         aaGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCAc           <  2:795670/149‑1 (MQ=255)
                                                                                                                                         aaGATGATAATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGGGAAGGCGCAAACGTGTCGATCATAGAGAACCGCAAAACCAGTGAAGCAGCGGAAGGTCCGTATCGCGAAGAAGGGATGATTGTTCAGCAATTCCAc           >  2:914150/1‑149 (MQ=255)
                                                                                                                                          aGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGGCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCAcc          <  2:481784/149‑1 (MQ=255)
                                                                                                                                              tatcatcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGtt      <  1:1502738/148‑1 (MQ=255)
                                                                                                                                               atcatcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTa     <  2:52648/149‑1 (MQ=255)
                                                                                                                                               atcatcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTa     >  1:647258/1‑149 (MQ=255)
                                                                                                                                                  atcCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGACCCGCAAACCCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGTATGATTCTTCAGCAATTCCACCCCTTAACg  >  2:1149045/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTACCG  >  NC_000913/3181152‑3181446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑