Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,109 A→T 8.1% intergenic (+157/+20) yajD → / ← tsx HNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913431,1090AT 8.1% 336.0 / 16.6 99intergenic (+157/+20)yajD/tsxHNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base A (44/47);  new base T (4/4);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 103 total aligned reads displayed.

AAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGAC  >  NC_000913/430963‑431257
                                                                                                                                                  |                                                                                                                                                    
aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATaaa                                                                                                                                                    <  4:839433/149‑1 (MQ=255)
aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATaaa                                                                                                                                                    <  2:563947/149‑1 (MQ=255)
aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATaaa                                                                                                                                                    <  3:102087/149‑1 (MQ=255)
 aaTTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAAt                                                                                                                                                   <  4:339130/149‑1 (MQ=255)
      ccTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCgg                                                                                                                                              <  1:214244/149‑1 (MQ=255)
      ccTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCgg                                                                                                                                              <  3:222551/149‑1 (MQ=255)
       cTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                               >  4:250310/1‑147 (MQ=255)
       cTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                               <  3:250310/147‑1 (MQ=255)
              gccgggattATGAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGATTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCa                                                                                                                                      >  2:118612/8‑149 (MQ=255)
               cTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttctga                                                                                                                                     <  1:319383/149‑5 (MQ=37)
               cTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttctga                                                                                                                                     <  3:416162/149‑5 (MQ=255)
               cTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAGATCCCGGCATTTTCAt                                                                                                                                     <  2:285472/149‑1 (MQ=255)
                tACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATa                                                                                                                                    >  3:49420/1‑149 (MQ=255)
                   gcgcaTCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATc                                                                                                                                 <  2:59934/145‑1 (MQ=255)
                    cTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCa                                                                                                                                <  4:206195/149‑1 (MQ=255)
                     ttATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAg                                                                                                                               >  1:40282/1‑149 (MQ=255)
                          aGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                          <  3:121152/149‑1 (MQ=255)
                          aGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                          <  3:677423/149‑1 (MQ=255)
                          aGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                          <  1:932868/149‑1 (MQ=255)
                             ccTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTa                                                                                                                       >  3:785294/1‑149 (MQ=255)
                              cTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCATTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                                      >  1:816951/1‑149 (MQ=255)
                                  gTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                                  >  1:784882/1‑149 (MQ=255)
                                  gTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                                  <  4:878329/149‑1 (MQ=255)
                                   ttATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACctac                                                                                                                 >  3:540947/1‑149 (MQ=255)
                                         cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                           <  3:738431/149‑1 (MQ=255)
                                         cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                           <  2:754920/149‑1 (MQ=255)
                                         cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                           <  2:805243/149‑1 (MQ=255)
                                         cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                           <  1:123168/149‑1 (MQ=255)
                                          aGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAgg                                                                                                          >  1:322963/1‑149 (MQ=255)
                                             aaTATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTaa                                                                                                       <  2:189003/149‑1 (MQ=255)
                                              atatatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAac                                                                                                      <  2:691908/149‑1 (MQ=255)
                                                   ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 >  3:699266/1‑149 (MQ=255)
                                                   ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 >  4:498032/1‑149 (MQ=255)
                                                   ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 >  1:505180/1‑149 (MQ=255)
                                                      aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGt                                                                                                                        <  3:746154/123‑1 (MQ=255)
                                                      aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGt                                                                                                                        >  4:746154/1‑123 (MQ=255)
                                                      aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAg                                                                                              >  2:237209/1‑149 (MQ=255)
                                                        tttACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGcc                                                                                            >  3:178610/1‑149 (MQ=255)
                                                            cGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCTCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTa                                                                                        <  1:765641/149‑1 (MQ=255)
                                                                cTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGc                                                                                    >  2:118037/1‑149 (MQ=255)
                                                                 ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 <  3:183299/135‑1 (MQ=255)
                                                                 ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                                 >  4:183299/1‑135 (MQ=255)
                                                                 ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCg                                                                                   >  3:121485/1‑149 (MQ=255)
                                                                   ttgttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTCTTAACAATCAGAAATGCCGGGAATAAAACCCAGCACTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCTCATCCGGTAGAGTGaa                                                                                 <  2:328915/149‑1 (MQ=38)
                                                                       tgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGtt                                                                             >  4:844646/1‑149 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  2:646212/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  3:59981/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  4:446961/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  2:502379/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  2:416835/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  4:692994/149‑1 (MQ=255)
                                                                          ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  2:322963/149‑1 (MQ=255)
                                                                           cGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAg                                                                         >  1:627264/1‑149 (MQ=255)
                                                                            ggACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGaa                                                                                                                             <  4:236419/96‑1 (MQ=255)
                                                                            ggACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGaa                                                                                                                             >  3:236419/1‑96 (MQ=255)
                                                                             gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgtt                                                                       >  1:919227/1‑149 (MQ=255)
                                                                             gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgtt                                                                       <  3:72632/149‑1 (MQ=255)
                                                                                aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTAcc                                                                                                             <  3:460722/108‑1 (MQ=255)
                                                                                aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTAcc                                                                                                             >  4:460722/1‑108 (MQ=255)
                                                                                aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgcc                                                                    >  3:85377/1‑149 (MQ=255)
                                                                                    cGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttg                                                                >  1:657039/1‑149 (MQ=255)
                                                                                     gTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgc                                                               <  4:85377/149‑1 (MQ=255)
                                                                                         aCGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGaa                                                           <  4:78713/149‑1 (MQ=255)
                                                                                         aCGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCAGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGaa                                                           >  2:277659/1‑149 (MQ=255)
                                                                                              gcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttca                                                      <  1:237209/149‑1 (MQ=255)
                                                                                                aTCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGGATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagt                                                    >  3:632438/1‑149 (MQ=255)
                                                                                                     gCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagttcc                                                 <  3:269454/147‑2 (MQ=255)
                                                                                                     gCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagttcc                                                 >  4:269454/1‑146 (MQ=255)
                                                                                                       aTGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          <  1:914742/120‑1 (MQ=255)
                                                                                                       aTGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                          >  2:914742/1‑120 (MQ=255)
                                                                                                        tGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtc                                            <  4:785294/149‑1 (MQ=255)
                                                                                                          aaCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgt                                          >  1:311225/1‑149 (MQ=255)
                                                                                                          aaCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgt                                          >  4:908808/1‑149 (MQ=255)
                                                                                                           aCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgtc                                         >  4:531916/1‑149 (MQ=255)
                                                                                                           aCAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCAtcgtc                                         >  4:871969/1‑149 (MQ=255)
                                                                                                            cAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                                  <  2:403699/75‑1 (MQ=255)
                                                                                                            cAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                                  >  1:403699/1‑75 (MQ=255)
                                                                                                             aaaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATcgtcgt                                       <  4:82043/149‑1 (MQ=255)
                                                                                                              aaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGtt                                      <  3:310928/149‑1 (MQ=255)
                                                                                                               aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTc                                     >  1:845661/1‑149 (MQ=255)
                                                                                                                gCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTcc                                    >  4:912987/1‑149 (MQ=255)
                                                                                                                  acacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCAc                                  >  3:791238/1‑149 (MQ=255)
                                                                                                                    acGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTg                                >  4:608855/1‑149 (MQ=255)
                                                                                                                     cGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGa                               <  1:277659/149‑1 (MQ=255)
                                                                                                                                aTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTATCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGTCTTTGCCGAAGTTCAGTTATGCATTGTCGTTCCTCTGACCACCGTCGTg                    >  1:802571/1‑149 (MQ=255)
                                                                                                                                  cAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGcc                  <  2:311225/149‑1 (MQ=255)
                                                                                                                                   aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa                 <  1:697654/149‑1 (MQ=255)
                                                                                                                                   aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa                 >  1:747681/1‑149 (MQ=255)
                                                                                                                                   aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa                 >  3:599419/1‑149 (MQ=255)
                                                                                                                                     aaaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGt               <  2:784882/149‑1 (MQ=255)
                                                                                                                                      aaTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTa              >  2:93880/1‑149 (MQ=255)
                                                                                                                                       aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCa                                   >  4:147892/1‑127 (MQ=255)
                                                                                                                                       aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCa                                   <  3:147892/127‑1 (MQ=255)
                                                                                                                                       aTGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTaa             >  2:141421/1‑149 (MQ=255)
                                                                                                                                        tGCCGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCATTCCACTGACCACCGTCGTGCCAGTAAc            >  1:524212/1‑149 (MQ=255)
                                                                                                                                          ccGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGa          <  4:49420/149‑1 (MQ=255)
                                                                                                                                           cGGGAATAAATCCCGGCATTTTCATATTCAGAAGTTGTAAGCAACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTACCGAAGTTCAGTTCTGCATCGTCGTTTCACGGACCCCCGTCGGGCCAGTAACGAg         <  2:604599/149‑1 (MQ=255)
                                                                                                                                            gggAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGc        >  1:553272/1‑149 (MQ=255)
                                                                                                                                               aaTAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTAc     <  2:51309/149‑1 (MQ=255)
                                                                                                                                                  aaaTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGaa  >  3:80386/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAGCTACGAC  >  NC_000913/430963‑431257

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