Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,165,154 T→A 7.5% M274L (ATG→TTG)  parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,165,1540TA 7.5% 265.9 / 8.4 88M274L (ATG→TTG) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base A (4/3);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCG  >  NC_000913/3165005‑3165301
                                                                                                                                                     |                                                                                                                                                   
ctTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCAt                                                                                                                                                     <  1:1282259/150‑1 (MQ=255)
    tCCAGCTCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGGTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgc                                                                                                                                                 <  1:317100/150‑1 (MQ=255)
    tCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgc                                                                                                                                                 >  1:252460/1‑150 (MQ=255)
      cAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgcgc                                                                                                                                               >  1:1057708/1‑150 (MQ=255)
      cAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgcgc                                                                                                                                               >  1:421464/1‑150 (MQ=255)
         aTCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCgcagc                                                                                                                                            <  1:1284825/150‑1 (MQ=255)
           cGGTGGTAGCGAAGAGGTTTTTAATAAACTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCCTCGGCAGCTTTTTGTTGGGCATTTGCGCAGCaa                                                                                                                                          >  2:938135/1‑150 (MQ=255)
             gtggtAGCGAAGAGGTGGTTCATCACCTGATGCATATCCACGCGGTTGGAACGCGGCACAATCACCTGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATGGGCAGCTTTTTGTTGCGCAAATGCgcaagaaaa                                                                                                                                        <  2:235168/150‑7 (MQ=255)
             gtggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAAtt                                                                                                                                        <  1:867594/150‑1 (MQ=255)
               ggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTg                                                                                                                                      >  1:1467227/1‑150 (MQ=255)
               ggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGCTGCGCAAATGCGcaacaaaaag                                                                                                                                      <  1:461152/150‑5 (MQ=255)
               ggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGCTGCGCAAATGCGcaacaaaaag                                                                                                                                      <  2:1260306/150‑5 (MQ=255)
                gtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGc                                                                                                                                     <  1:373683/150‑1 (MQ=255)
                gtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGc                                                                                                                                     <  1:319126/150‑1 (MQ=255)
                 tAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTATCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCCCATTTGCGCAGCAATTTGCt                                                                                                                                    >  2:694601/1‑150 (MQ=255)
                   gCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTcc                                                                                                                                  <  1:101297/150‑1 (MQ=255)
                    cGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCa                                                                                                                                 >  2:1273343/1‑150 (MQ=255)
                     gAAGAGGTGGTTCATCACCTGCTCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAg                                                                                                                                >  2:817382/1‑150 (MQ=255)
                      aaGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCACCCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGt                                                                                                                               >  2:513477/1‑150 (MQ=255)
                          ggtggtTCAGCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTAcgc                                                                                                                           >  1:699920/1‑150 (MQ=255)
                            tggtTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCCCCTGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTAcgcgc                                                                                                                         >  2:712229/1‑150 (MQ=255)
                            tggtTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTAcgcgc                                                                                                                         <  2:379867/150‑1 (MQ=255)
                               ttCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCa                                                                                                                      >  1:1183455/1‑150 (MQ=255)
                                  atcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCActt                                                                                                                   >  1:709014/1‑148 (MQ=255)
                                   tcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTg                                                                                                                  <  2:252460/150‑1 (MQ=255)
                                    aaCCTGATCCATATCCACGCGGTTGGACCGCGGCACAATCACCAGGCGCGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCCCCTGa                                                                                                                 >  1:277418/2‑150 (MQ=255)
                                      ccTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGaaa                                                                                                               <  1:1228739/150‑1 (MQ=255)
                                        tGATCCATATCCACGCGGTTGGAACTCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAAcc                                                                                                             >  2:1340485/1‑150 (MQ=255)
                                        tGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAAcc                                                                                                             >  1:308901/1‑150 (MQ=255)
                                          aTCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCCCCTGAAACCTg                                                                                                           >  2:960415/1‑150 (MQ=255)
                                           tCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGa                                                                                                          >  2:1418578/1‑150 (MQ=255)
                                                atCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCgg                                                                                                     <  2:119748/150‑1 (MQ=255)
                                                   cACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCaa                                                                                                  >  1:470328/1‑150 (MQ=255)
                                                         gTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCt                                                                                            <  1:754588/150‑1 (MQ=255)
                                                         gTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCt                                                                                            <  2:470328/150‑1 (MQ=255)
                                                          ttGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTg                                                                                           <  2:458656/150‑1 (MQ=255)
                                                          ttGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTg                                                                                           <  1:1273343/150‑1 (MQ=255)
                                                           tGGAACGCGGCACAATCACAAGGCGGGTCGGGTTCTCGTGGTCAGAATCATCGCGCAGATCGTCAACCATCGGCAGCTATTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTCAAACCTGATGCGGCAATGCGCTGa                                                                                          >  2:816550/1‑150 (MQ=255)
                                                             gAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATa                                                                                        >  2:790793/1‑150 (MQ=255)
                                                                      acaATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACcgcg                                                                               <  2:169507/150‑1 (MQ=255)
                                                                       caATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCgcgc                                                                              <  2:1265214/150‑1 (MQ=255)
                                                                        aaTCACCAGTCGGGTAGGGTTCTCGTGGTCAGATTCATAGAGCAGATCGTCAACAATCGGCAGCTTTTAGTTGAGCAATTGCGCAGAAATTTGCTACACTACGCACGCACCTGAAACCAGATGCGGCAATGAGCTGATAACCACCGCGcc                                                                             >  2:1386055/1‑150 (MQ=255)
                                                                         aTCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCa                                                                            >  2:973417/1‑150 (MQ=255)
                                                                          tCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCACCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCAt                                                                           >  2:957683/1‑150 (MQ=255)
                                                                           cACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCGACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCAtc                                                                          >  1:199964/1‑150 (MQ=255)
                                                                                ggCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATcttctt                                                                     >  2:221667/1‑150 (MQ=255)
                                                                                    ggtcgggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCttctttcta                                                                 >  2:4919/1‑149 (MQ=255)
                                                                                    ggtcgggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTtctttctt                                                                 <  1:1001687/150‑1 (MQ=255)
                                                                                      tcgggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTcc                                                               <  1:168198/150‑1 (MQ=255)
                                                                                        gggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCcac                                                             >  2:451832/1‑150 (MQ=255)
                                                                                        gggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCcac                                                             >  2:373301/1‑150 (MQ=255)
                                                                                           ttCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCACATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACAcc                                                          <  2:710535/150‑1 (MQ=255)
                                                                                             ctcGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGTGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACcgc                                                        >  1:981661/1‑150 (MQ=255)
                                                                                               cGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATcttctt                                                                     >  1:983547/1‑135 (MQ=255)
                                                                                               cGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATcttctt                                                                     <  2:983547/135‑1 (MQ=255)
                                                                                                gttgTCAGATTCATGGCGCAGATCGTCAACCATCGGGAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACATGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACcgcgcg                                                     <  2:942028/147‑1 (MQ=255)
                                                                                                  ggTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCa                                                   <  2:995841/150‑1 (MQ=255)
                                                                                                  ggTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCa                                                   >  2:477465/1‑150 (MQ=255)
                                                                                                    tCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATa                                                 >  2:1060650/1‑150 (MQ=255)
                                                                                                       gATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGc                                              <  1:221667/150‑1 (MQ=255)
                                                                                                       gATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGAGTCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGc                                              >  1:816971/1‑150 (MQ=255)
                                                                                                          tCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCCCCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACt                                           >  2:1449577/1‑150 (MQ=255)
                                                                                                           cATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTg                                          >  1:938368/1‑150 (MQ=255)
                                                                                                             tCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGaa                                        >  2:452835/1‑150 (MQ=255)
                                                                                                                cgcAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCa                                     >  1:1363687/1‑150 (MQ=255)
                                                                                                                cgcAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCa                                     >  2:942670/1‑150 (MQ=255)
                                                                                                                cgcAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCa                                     <  2:421464/150‑1 (MQ=255)
                                                                                                                     aTCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTcc                                <  1:1186420/150‑1 (MQ=255)
                                                                                                                       cGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGt                              >  1:920073/1‑150 (MQ=255)
                                                                                                                       cGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGt                              <  1:477465/150‑1 (MQ=255)
                                                                                                                       cGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGt                              >  2:630797/1‑150 (MQ=255)
                                                                                                                            aCCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCg                         <  2:1467227/150‑1 (MQ=255)
                                                                                                                             ccATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGt                        <  1:1340485/150‑1 (MQ=255)
                                                                                                                              cATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTa                       <  1:973417/150‑1 (MQ=255)
                                                                                                                                 cGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGAt                    >  1:602300/1‑150 (MQ=255)
                                                                                                                                   gCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGAttt                  <  1:630797/150‑1 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAACCCTGATTACGAAATGCGCTAAAAACACCGGCGCACTCTTCTTTCTTCAAGACCGCGCGCATACGAACTCAACACAGTACGGTCTCGTAGCTTTTAc               >  2:688426/1‑150 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAAcc                                      >  1:1443272/1‑127 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAAcc                                      <  2:1443272/127‑1 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTAc               >  1:190100/1‑150 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTAc               <  2:308901/150‑1 (MQ=255)
                                                                                                                                       caaattgctgCGCAAATGCGCAACAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACg              >  1:939795/5‑150 (MQ=255)
                                                                                                                                       cTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACTTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACg              <  2:709014/150‑1 (MQ=255)
                                                                                                                                       cTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACg              <  1:790793/150‑1 (MQ=255)
                                                                                                                                        aaattgctgCGCAAATGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg             >  2:73762/4‑150 (MQ=255)
                                                                                                                                             gttgCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTTTCGTAGATTTTACGGAtctc        >  1:1331701/1‑150 (MQ=255)
                                                                                                                                               tgCGCAAATGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCgg      >  1:316149/1‑150 (MQ=255)
                                                                                                                                                gcgcATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTGCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgc     >  2:1027238/1‑150 (MQ=255)
                                                                                                                                                   cATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgcgcg  <  1:1418578/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCG  >  NC_000913/3165005‑3165301

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