Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,877,557 G→C 8.1% intergenic (‑93/+147) yidB ← / ← gyrB DUF937 family protein/DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,877,5570GC 8.1% 268.8 / 13.8 88intergenic (‑93/+147)yidB/gyrBDUF937 family protein/DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (41/38);  new base C (3/4);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATT  >  NC_000913/3877409‑3877705
                                                                                                                                                    |                                                                                                                                                    
tatTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTAGAGATTTGCGGTATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGGgatg                                                                                                                                                      >  1:41293/1‑149 (MQ=255)
 atTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAgatga                                                                                                                                                     <  2:374642/149‑1 (MQ=255)
     ccccGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCAtctc                                                                                                                                                 <  1:279455/149‑1 (MQ=255)
         gccacGCCTATCAAAAAGACCCCGCAAACTACATCAACCAGCCCCATGATAACCACCCATAAAGCCTCACGTGTCGACCCAAACCGTATCGCTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGcc                                                                                                                                             <  1:765856/145‑1 (MQ=255)
         gCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGcc                                                                                                                                             <  1:536431/149‑1 (MQ=255)
         gCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGGAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATCATCTCCGcc                                                                                                                                             <  2:75978/149‑1 (MQ=255)
           aTCGCCTTTCAGAAAGGCACCGACAACTTCCTCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCt                                                                                                                                           <  1:22409/149‑1 (MQ=255)
              gCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCtttt                                                                                                                                        <  1:758840/149‑1 (MQ=255)
               ccgttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTc                                                                                                                                       <  1:280118/146‑1 (MQ=255)
                cTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGGTGATCTCCGCCCGTTTCt                                                                                                                                      >  2:256079/1‑149 (MQ=255)
                cTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCt                                                                                                                                      <  2:237594/149‑1 (MQ=255)
                cTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCt                                                                                                                                      <  1:422389/149‑1 (MQ=255)
                 tttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCtt                                                                                                                                     >  2:421779/1‑149 (MQ=255)
                  ttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                    <  1:25013/149‑1 (MQ=255)
                  ttCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTta                                                                                                                                    <  2:482856/149‑1 (MQ=255)
                     aGAAAGGCACCGACACCTTCATCAAAAGGCCCCATTATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCTCTTTTCTTATAg                                                                                                                                 <  1:421910/149‑1 (MQ=255)
                     aGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAg                                                                                                                                 >  1:374799/1‑149 (MQ=255)
                              ccGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGGGATGCTCTCCGCCCTTTTCTTATAGCTTCTTGcc                                                                                                                        >  1:249664/1‑149 (MQ=255)
                                gACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCCGATATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCgg                                                                                                                      >  1:224000/1‑149 (MQ=255)
                                      ttCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATgggg                                                                                                                >  1:43677/1‑149 (MQ=255)
                                            aaaaaGCCCCATGATAATCACTTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGaa                                                                                                          <  1:602932/149‑1 (MQ=255)
                                            aaaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGaa                                                                                                          <  1:356597/149‑1 (MQ=255)
                                            aaaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGaa                                                                                                          <  2:789185/149‑1 (MQ=255)
                                            aaaaaGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGaa                                                                                                          <  1:470568/149‑1 (MQ=255)
                                                aGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGcc                                                                                                      <  1:89451/149‑1 (MQ=255)
                                                     cATTATAATCACCTGTCAAGCGATACGTAACGACCAAAAAAGTAGATATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 <  1:543357/149‑1 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 >  1:56714/1‑149 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 >  2:139493/1‑149 (MQ=255)
                                                     cATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCAGATGCGGCGTGAACGCCTTATc                                                                                                 >  2:658853/1‑149 (MQ=255)
                                                     cATGATAAACGCCTGAAAAGCGTCACGTGTGGACCCAAACAGGATAGCGTTGCGGGATGATAATTGCAGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATCGCTTCTTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                 <  1:452922/149‑1 (MQ=255)
                                                        gATAATCACCTGTAAAGCGTTAGGTGTTGACCCAAAACGTATAGATTTGCGGGATGACAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCATTTTCTTCTAGCTTCTTGCCGGCTGCGGCGTGAACGCCTTATCCgg                                                                                              <  2:598558/149‑1 (MQ=255)
                                                          tAATCACCTGTAAAGCGTTACGTTTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                            <  2:238853/149‑1 (MQ=255)
                                                            aTCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATTATCTCCGCCCTTTTCTTATAGCTTATTGCCCGAGGGGGCGTGAACGCCTTATCCGGCCTa                                                                                          <  2:569675/149‑1 (MQ=255)
                                                            aTCACCTGTAAAGCGTGACGGGTTGACCCAAGAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCAGCGTGAACGCCTTATCCGGCCTa                                                                                          <  1:613051/149‑1 (MQ=255)
                                                             tCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                         >  1:214864/1‑149 (MQ=255)
                                                             tCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                         >  1:98979/1‑149 (MQ=255)
                                                               aCCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATACTTGCGGGTTGAAAAAATAGAAAGGGCGGGGATGCACCGTGAAGTGGGGTGAGGGCAGCGTGCCGGAGGCGGGGTGAACGCCTTATCCGGCCTAcca                                                                                       >  2:366823/1‑147 (MQ=255)
                                                                 cTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                     <  1:444386/149‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATGGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  >  2:224180/1‑149 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  <  2:151325/149‑1 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  >  2:608962/1‑149 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCGTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATGCGGCCTACCAAATCa                                                                                  >  2:747419/1‑149 (MQ=255)
                                                                    tAAAGCGTTACGTGTTGACCCAAAAAGGATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCa                                                                                  <  2:796538/149‑1 (MQ=255)
                                                                                    gtcccAAAAAGCATAGATGTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                                                  <  1:287508/147‑1 (MQ=255)
                                                                                    gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAt                                                                  <  1:608083/149‑1 (MQ=255)
                                                                                     aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                 >  1:714829/1‑149 (MQ=255)
                                                                                       ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAGTACATTGc                                                               >  1:245218/1‑149 (MQ=255)
                                                                                       ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAATCAATACATTGc                                                               >  1:327086/1‑149 (MQ=255)
                                                                                        cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCa                                                              >  1:514397/1‑149 (MQ=255)
                                                                                        cAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGACGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATTCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCa                                                              <  2:442069/149‑1 (MQ=255)
                                                                                              gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGGGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAttt                                                        >  2:402332/1‑149 (MQ=255)
                                                                                               tataGATTTGCGGGATGATAATTGCTGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                                                       >  1:401625/1‑149 (MQ=255)
                                                                                                  tgATTTGTGGGATGATGATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGACGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGt                                                    <  2:672387/148‑1 (MQ=255)
                                                                                                    aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCAAAATCATGAAAATTCAATACATTTCAAGATTTACGTAg                                                  >  1:55246/1‑149 (MQ=255)
                                                                                                    aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACAATCATGAAAATTCAATACATTGCAAGATTTTCTTAg                                                  >  2:353491/1‑149 (MQ=255)
                                                                                                    aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAg                                                  >  1:607852/1‑149 (MQ=255)
                                                                                                    aTTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAg                                                  >  1:607841/1‑149 (MQ=255)
                                                                                                     tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTGTCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                                                 >  2:399050/1‑149 (MQ=255)
                                                                                                     tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAgg                                                 >  2:657477/1‑149 (MQ=255)
                                                                                                      ttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                <  2:566799/149‑1 (MQ=255)
                                                                                                      gtGCGGGATGATAATGGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                <  1:140086/148‑1 (MQ=255)
                                                                                                         cGGGTTGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACACAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                             >  2:373282/1‑149 (MQ=255)
                                                                                                         cGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                             >  1:180791/1‑149 (MQ=255)
                                                                                                         cGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTAGTAGGCCTg                                             >  1:555442/1‑149 (MQ=255)
                                                                                                          gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                            <  1:421779/149‑1 (MQ=255)
                                                                                                                atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACCTTGCAAGATTTTCGTAGGCCTGATAAGCg                                      >  1:322683/1‑149 (MQ=255)
                                                                                                                 tAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                     <  1:195874/149‑1 (MQ=255)
                                                                                                                 tAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                     <  1:192533/149‑1 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                  >  2:235457/1‑149 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                  >  1:608258/1‑149 (MQ=255)
                                                                                                                    ttGCGGATTGCAGAAATAAAAAGGGCCGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCAGGATGCGGCGTGAACGCCTTATCCGGCCTACAACATCATGAAAATTCAATAAATTGCACCCTTTTATGACGCCTGATAAGCGTCgc                                  >  2:362259/1‑149 (MQ=25)
                                                                                                                       aggATTGCAGAAATAAAAAGGGCGGAGATGACCTCCGCTCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCa                               <  1:797091/148‑1 (MQ=255)
                                                                                                                         gactGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCTGGCCTACAAAATCCTGAAAATTCAATACATTGCAAGAGTTTCGTAGGCCTGATAAGCGTAGCGCATc                             <  1:489103/146‑1 (MQ=255)
                                                                                                                           ttGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:436520/1‑149 (MQ=255)
                                                                                                                           ttGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:482020/1‑149 (MQ=255)
                                                                                                                            tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAgg                          <  2:736411/149‑1 (MQ=255)
                                                                                                                              cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:512407/1‑149 (MQ=255)
                                                                                                                                tataagaaaaGGGCGGAGATCATCTCCGCGCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCACTACATTGCAAGATTTTCGTAGGCGGGATAAGCGAAGCGCATCAGGCACg                      >  1:487087/7‑149 (MQ=21)
                                                                                                                                tataagaaaaGGGCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACg                      >  1:729356/7‑149 (MQ=255)
                                                                                                                                      aaaaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCTGAGGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCa                <  1:150384/149‑1 (MQ=255)
                                                                                                                                         aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATgg             >  2:413832/1‑149 (MQ=255)
                                                                                                                                         aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATgg             >  2:264704/1‑149 (MQ=255)
                                                                                                                                          gggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCTTGTACTCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGt            <  1:224180/149‑1 (MQ=255)
                                                                                                                                          ccacgctgatgatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGt            >  2:135463/8‑149 (MQ=255)
                                                                                                                                           ggCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAGAATTCAATACATTGCAAGATTTTCGTAGGCCGGATAAGCGTCGCGCATCGGGCACGCCACCATGGtt           >  2:620615/1‑149 (MQ=25)
                                                                                                                                            gCGGAGATCATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGAAGCGCATCAGGCACGCTCGCATGGTTa          >  1:722021/1‑149 (MQ=255)
                                                                                                                                             cGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCCGGCACGCTCGCATGGTTAg         >  1:293309/1‑149 (MQ=255)
                                                                                                                                              ggggatgatCTCCGCCCTTTTCTTATAGCTGCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAGAAGCGTAGCGAATCAGGCACGCTCGCATGGTTAgc        <  2:606523/146‑1 (MQ=255)
                                                                                                                                              ggAGATTATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCCGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgc        <  2:352342/149‑1 (MQ=255)
                                                                                                                                              ggAGATGATCTCCGCCCTTTTCTTCTAGCTTCTTGCCGTGTGCGTCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgc        <  1:617242/149‑1 (MQ=255)
                                                                                                                                              ggAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgc        <  2:482263/149‑1 (MQ=255)
                                                                                                                                                agaTGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAgcgc      <  2:539259/149‑1 (MQ=255)
                                                                                                                                                    gatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTGGGCTGGATAATCGTAGCGCATCAGGCAAGCTCGCAGGGTTAGGGCCAtt  >  2:767668/1‑149 (MQ=25)
                                                                                                                                                    |                                                                                                                                                    
TATTTCCCCGCATCGCCTTTCAGAAAGGCACCGACAACTTCATCAAAAAGCCCCATGATAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATT  >  NC_000913/3877409‑3877705

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