Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,407,947 G→A 5.3% E432K (GAA→AAA)  rnr → exoribonuclease R, RNase R

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,407,9470GA 5.3% 335.7 / 2.2 110E432K (GAA→AAA) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base G (55/49);  new base A (4/2);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCA  >  NC_000913/4407798‑4408096
                                                                                                                                                     |                                                                                                                                                     
      ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTTAGAGCGaagaag                                                                                                                                                 >  2:1526131/1‑150 (MQ=255)
      ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                                                                                                 >  1:978043/1‑150 (MQ=255)
      ccAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGaagaag                                                                                                                                                 >  2:439926/1‑150 (MQ=255)
          aGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGaa                                                                                                                                             <  1:1511469/150‑1 (MQ=255)
          aGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCATGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGaa                                                                                                                                             <  1:886170/150‑1 (MQ=255)
           gTCTGGCATATTCTTGAGGGCGATCTGGATCTGCGCGATGAGTACGCCCCGGTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATATAGTGATGGATATAGCCCGTGATGATCTCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAg                                                                                                                                            <  1:82837/150‑1 (MQ=255)
           gTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCTTAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGAGCGCAGTTGAGAGCGAAGAGGCGAAg                                                                                                                                            >  1:963103/1‑150 (MQ=255)
           gTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCAAAGAAGCGAAg                                                                                                                                            <  2:547471/150‑1 (MQ=255)
           gTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCAAAGAAGCGAAg                                                                                                                                            <  2:540317/150‑1 (MQ=255)
               ggCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCa                                                                                                                                        >  1:1445936/1‑150 (MQ=255)
                   tatTCTGCAGGGCGATCAGGATCTTCTCGAGCAGTACGCCCCGCGGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAAGGCGGTGGGATCTCATTTGAGAGCGAAGAAGTGAAGTTCAtttt                                                                                                                                    <  1:281697/150‑1 (MQ=255)
                     ttCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCa                                                                                                                                  <  2:1219302/150‑1 (MQ=255)
                      tCTGCAGGGCGATGAGGATCTGCGCTAGCAGTACGCCCCGCTGGTTAAGCATCTCGACGAGTTGCATAACCTCTATAAAGTGCTGGTTAAAGCCCGTGAAGAACGCGGTGTGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCaa                                                                                                                                 <  2:389854/150‑1 (MQ=255)
                      tCTGCAGGGCGATCAGGATCTGCGCGATCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCaa                                                                                                                                 <  2:1511006/150‑1 (MQ=255)
                       cTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAAc                                                                                                                                <  1:1259865/150‑1 (MQ=255)
                        tGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACg                                                                                                                               >  1:434353/1‑150 (MQ=255)
                        tGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACg                                                                                                                               <  1:453861/150‑1 (MQ=255)
                            gggCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGa                                                                                                                           >  2:356366/1‑150 (MQ=255)
                               cGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGATGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACg                                                                                                                        >  1:1078439/1‑150 (MQ=255)
                               cGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACg                                                                                                                        >  1:1821770/1‑150 (MQ=255)
                                      gATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAtt                                                                                                                 <  2:84863/150‑1 (MQ=255)
                                      gATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAtt                                                                                                                 >  1:1800689/1‑150 (MQ=255)
                                       aTCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGGATCTCATTTGATAACGAAGAAGCGAAGTTCATTTTCAACGCTGCACGCCGTATTg                                                                                                                >  2:1630693/1‑150 (MQ=255)
                                        tCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGa                                                                                                               <  1:1373122/150‑1 (MQ=255)
                                          tGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACAATAAATAAAAAAAAACAGAAGTTCATTTTCAACGCTGAACGCCGTATTGAAc                                                                                                             >  2:24302/1‑150 (MQ=255)
                                          tGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAAc                                                                                                             >  1:633570/1‑150 (MQ=255)
                                          tGCGCGAGCAGTAAGCCCCGCTGGTTAAGCAGCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAAc                                                                                                             <  1:28331/150‑1 (MQ=255)
                                           gcgcGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACg                                                                                                            >  1:1348066/1‑150 (MQ=255)
                                              cgAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTAt                                                                                                         <  2:314249/150‑1 (MQ=255)
                                               gAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTAGTGAACGTATc                                                                                                        >  1:581284/1‑150 (MQ=255)
                                                aGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCg                                                                                                       <  2:721827/150‑1 (MQ=255)
                                                 gCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGa                                                                                                      >  1:876012/1‑150 (MQ=255)
                                                   aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAAc                                                                                                    <  2:779968/150‑1 (MQ=255)
                                                   aGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAAc                                                                                                    <  2:1445936/150‑1 (MQ=255)
                                                    gTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCGGGATTAAGCCCGTGAAGAACGCGGTGGGTTGTCATGTGAGAGCGCCGCCGCGAAGTTCAGTATCACCGCTGACCGCCGTATAGAACGAATCGAACa                                                                                                   >  1:153384/1‑150 (MQ=255)
                                                      aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGCACGCGGTGGGATCTCATTTGAGAGCGAAGCAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGCACGTATCGAACAGa                                                                                                 >  1:595878/1‑150 (MQ=255)
                                                      aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCCACGCTGACCGCCGTATTGAACGTATCgtccagc                                                                                                 >  2:269754/1‑144 (MQ=255)
                                                      aCGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGa                                                                                                 >  1:1823869/1‑150 (MQ=255)
                                                            cGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGc                                                                                           <  2:392454/150‑1 (MQ=255)
                                                             gCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAACGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCg                                                                                          <  2:1328770/150‑1 (MQ=255)
                                                                   tAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAacga                                                                                    >  1:1770685/1‑150 (MQ=255)
                                                                     aGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGAGAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCCGCGTAacgacg                                                                                  >  2:319705/1‑150 (MQ=255)
                                                                      gCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGCCACAGCGTAACGAcgc                                                                                 >  1:147691/1‑150 (MQ=255)
                                                                      gCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGAcgc                                                                                 >  2:382703/1‑150 (MQ=255)
                                                                           tcGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGcaca                                                                            >  1:555537/1‑150 (MQ=255)
                                                                            cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGCCCCAGCGTAACGACGCGCACaa                                                                           >  2:442386/1‑150 (MQ=255)
                                                                            cGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTAAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACaa                                                                           >  1:143207/1‑150 (MQ=255)
                                                                             gAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCg                                                                                          <  1:937286/134‑1 (MQ=255)
                                                                             gAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACACAGCg                                                                                          >  2:937286/1‑134 (MQ=255)
                                                                                 tgTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaa                                                                      <  2:1559176/149‑1 (MQ=255)
                                                                                   ttGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaatt                                                                    <  1:78257/150‑1 (MQ=255)
                                                                                   ttGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaatt                                                                    <  2:1117269/150‑1 (MQ=255)
                                                                                   ttGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCAGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaatt                                                                    <  1:350004/150‑1 (MQ=255)
                                                                                   ttGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACACTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAaattaatt                                                                    <  1:1520813/150‑1 (MQ=255)
                                                                                    tGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTg                                                                   >  1:1328695/1‑150 (MQ=255)
                                                                                     gCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCACATTTCAGAGCAAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGCACAGACCCAGCGTAACGCCGCGCACAAATTAATTGa                                                                  >  2:1046470/1‑150 (MQ=255)
                                                                                      cATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGaa                                                                 >  1:986675/1‑150 (MQ=255)
                                                                                         aaCCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAagag                                                              <  2:896192/150‑1 (MQ=255)
                                                                                           ccTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAGCAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTg                                                            >  1:555242/1‑150 (MQ=255)
                                                                                            ctctATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGc                                                           <  1:1526131/150‑1 (MQ=255)
                                                                                                ataAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGCACGCCGTATTGAACGTATCGAACATCCCCAGCGTAACGACGCGCAAAAATTAATTGAAGATTGCATGa                                                       >  1:111541/1‑150 (MQ=255)
                                                                                                   aaGTGCTGGATAAAGCCCGTGAAGAACTCGGTGGGATCACATTTGAGAGCAAAAAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTc                                                    >  2:1345926/1‑150 (MQ=255)
                                                                                                   aaGTGCTGGATAAAGCCCGTGAAAAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGCACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGCCGCGCACAAATTAATTGAAGAGTGCATGATTc                                                    >  2:1353322/1‑150 (MQ=255)
                                                                                                    aGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCt                                                   >  1:1664246/1‑150 (MQ=255)
                                                                                                       gCTGGATAAAGCCCGTGAAGAACTCGGTGGTATCACAATAAATAAAAAAAAAAAGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGCCGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                                >  2:21318/1‑150 (MQ=255)
                                                                                                       gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                                <  1:941891/150‑1 (MQ=255)
                                                                                                       gCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGc                                                >  1:365153/1‑150 (MQ=255)
                                                                                                        cTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCg                                               <  2:1348066/150‑1 (MQ=255)
                                                                                                          ggATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGaa                                             >  1:879635/1‑150 (MQ=255)
                                                                                                           gATATAGCCAGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACACTGAACGCCGTATTGAACGTTTCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAat                                            <  1:894456/150‑1 (MQ=255)
                                                                                                            aTAAAGCCCGTGAAGAACGCGGTGGTATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTTTTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAata                                           <  2:137992/150‑1 (MQ=255)
                                                                                                            aTAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCGGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAata                                           <  2:1176724/150‑1 (MQ=255)
                                                                                                             tAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGAGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAatat                                          >  1:865835/1‑150 (MQ=255)
                                                                                                                aGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAATGCATGATTCTGGCGAATAtctc                                       >  2:1275384/1‑150 (MQ=255)
                                                                                                                 gCCCGTGAAGAACGCGGTGGGATCTAATTTGAGATCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCTAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGTCGAAAATCTcg                                      >  2:1781516/1‑150 (MQ=255)
                                                                                                                    cGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcg                                   >  2:827937/1‑150 (MQ=255)
                                                                                                                    agTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcg                                   >  1:756500/2‑150 (MQ=255)
                                                                                                                     gTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGCACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTcggcgg                                  <  1:590778/150‑1 (MQ=255)
                                                                                                                       gaagaaCGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGgcggcg                                >  1:943490/1‑150 (MQ=255)
                                                                                                                            aCGCGGTGGGAACTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATAGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAAATAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGttt                           >  1:1176435/1‑150 (MQ=255)
                                                                                                                              gcgGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCg                         <  2:1577848/150‑1 (MQ=255)
                                                                                                                               cgGTGGCATCTCATTTGAGATCGACGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGTACGTATCGAACAGACCCAGCTTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGt                        <  1:319705/150‑1 (MQ=255)
                                                                                                                                 gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTg                      <  2:1821770/150‑1 (MQ=255)
                                                                                                                                 gTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTg                      <  1:1411474/150‑1 (MQ=255)
                                                                                                                                  tGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGGGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTga                     <  1:1670253/150‑1 (MQ=255)
                                                                                                                                      aTCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCCACGCTGAACGCCGTATTGAACGTATCTCTCAGCCCCCGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTGGAGaaa                 >  2:863527/1‑150 (MQ=255)
                                                                                                                                        ttcttttGAGAGCGAAGTAGCGAAGTTGATTTTCATCGCTGATCGTCGTATTGAACGTATCGAACAGACCCAGCGTAACGTCGCGCACATATTAATTGAAGAGTGCATGATTCTGGCTAATATCTCGGCGGCGCGTTTCGTTGAGAAAGc               <  1:24302/146‑1 (MQ=255)
                                                                                                                                        ctcATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGc               <  1:356366/150‑1 (MQ=255)
                                                                                                                                         tcATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCg              >  2:69472/1‑150 (MQ=255)
                                                                                                                                         tcATTAGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTTAACGTATCGATCAGACCCAGCGTAACGACGCGCACAAATAAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCg              <  1:323728/150‑1 (MQ=255)
                                                                                                                                          cATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGa             >  1:1085567/1‑150 (MQ=255)
                                                                                                                                          aaTTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGa             <  1:1098779/149‑1 (MQ=255)
                                                                                                                                              tGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGATGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGa         >  2:1517115/1‑150 (MQ=255)
                                                                                                                                               gagaGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTAGAGAAAGCGAAAGaa        >  1:1658338/1‑150 (MQ=255)
                                                                                                                                                  agCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCg     >  1:974874/1‑150 (MQ=255)
                                                                                                                                                   gCGAAGAAGCGAAGTTCATTTTCACCGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATAACTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCgg    >  1:130015/1‑150 (MQ=255)
                                                                                                                                                   gCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCgg    >  2:1790336/1‑150 (MQ=255)
                                                                                                                                                    cGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGc   >  1:326534/1‑150 (MQ=255)
                                                                                                                                                    cGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGc   >  2:1470007/1‑150 (MQ=255)
                                                                                                                                                     gaagaaGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCa  >  2:1756087/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTACACCAAAGTCTGGCATATTCTGCAGGGCGATCAGGATCTGCGCGAGCAGTACGCCCCGCTGGTTAAGCATCTCGAAGAGTTGCATAACCTCTATAAAGTGCTGGATAAAGCCCGTGAAGAACGCGGTGGGATCTCATTTGAGAGCGAAGAAGCGAAGTTCATTTTCAACGCTGAACGCCGTATTGAACGTATCGAACAGACCCAGCGTAACGACGCGCACAAATTAATTGAAGAGTGCATGATTCTGGCGAATATCTCGGCGGCGCGTTTCGTTGAGAAAGCGAAAGAACCGGCA  >  NC_000913/4407798‑4408096

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