Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,400,432 C→A N48K (AAC→AAA)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,400,4320CA98.5% 230.1 / ‑3.4 67N48K (AAC→AAA) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/32);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGG  >  NC_000913_3_bme_pgi/4400286‑4400577
                                                                                                                                                  |                                                                                                                                                 
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                <  3:154692/150‑1 (MQ=255)
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                <  1:9069/150‑1 (MQ=255)
agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACg                                                                                                                                                >  4:296859/1‑150 (MQ=255)
       cTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGcc                                                                                                                                         >  3:215163/1‑150 (MQ=255)
       cTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGcc                                                                                                                                         <  3:296859/150‑1 (MQ=255)
        tAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCa                                                                                                                                        >  4:46698/1‑150 (MQ=255)
          aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGa                                                                                                                                      >  1:564810/1‑150 (MQ=255)
          aGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGa                                                                                                                                      >  1:400591/1‑150 (MQ=255)
                 aTCTTTACAAGGTCCGTTCCTGAACGCACTGCGTCTGGTACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTa                                                                                                                               <  2:153632/150‑1 (MQ=255)
                 aTCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTa                                                                                                                               <  3:310248/150‑1 (MQ=255)
                  tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTAc                                                                                                                              <  1:615154/150‑1 (MQ=255)
                    tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACaa                                                                                                                            >  4:365328/1‑150 (MQ=255)
                    tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACaa                                                                                                                            <  3:46698/150‑1 (MQ=255)
                    tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACaa                                                                                                                            >  3:117404/1‑150 (MQ=255)
                        cAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGCTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAc                                                                                                                        >  3:17466/1‑150 (MQ=255)
                           gATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCAcgcg                                                                                                                     >  1:75150/1‑150 (MQ=255)
                              ccGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGAtt                                                                                                                  <  3:26068/150‑1 (MQ=255)
                                  tCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                                              <  1:77500/150‑1 (MQ=255)
                                  tCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa                                                                                                              <  4:17466/150‑1 (MQ=255)
                                   ccTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                                                                                                             <  3:12601/150‑1 (MQ=255)
                                   ccTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                                                                                                             <  1:608648/150‑1 (MQ=255)
                                     tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  2:693968/1‑150 (MQ=255)
                                     tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg                                                                                                           >  2:693956/1‑150 (MQ=255)
                                       aaCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACaa                                                                                                                            <  2:632462/131‑1 (MQ=255)
                                       aaCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACaa                                                                                                                            >  1:632462/1‑131 (MQ=255)
                                            aCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTATTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTccc                                                                                                    <  2:681327/150‑1 (MQ=255)
                                              tGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAACAAACGGTCACCCAACTGGTTTACAAGACCGCGATTTCTACTGTTGTCCCGt                                                                                                  >  3:198588/1‑150 (MQ=255)
                                               gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttgt                                                                                                       >  2:291053/1‑144 (MQ=255)
                                               gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgttgt                                                                                                       <  1:291053/144‑1 (MQ=255)
                                                      gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                          <  2:564810/150‑1 (MQ=255)
                                                          gtgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGttt                                                                                      <  3:435901/150‑1 (MQ=255)
                                                           tgtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtc                                                                                     >  4:385208/1‑150 (MQ=255)
                                                            gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                    <  1:693968/150‑1 (MQ=255)
                                                            gtTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                                                                    <  1:693956/150‑1 (MQ=255)
                                                             ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAA‑ACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtca                                                                                  >  2:312865/1‑150 (MQ=255)
                                                                cAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatc                                                                                >  1:689894/1‑150 (MQ=255)
                                                                   tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                             >  2:251663/1‑150 (MQ=255)
                                                                   tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                             >  2:251574/1‑150 (MQ=255)
                                                                   tttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                                                             <  2:689894/150‑1 (MQ=255)
                                                                       tatttattTGGTGAATGGTATTAAGCTGCGAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAg                                                                                                                           <  1:433972/100‑1 (MQ=255)
                                                                       tatttattTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAg                                                                                                                           >  2:433972/1‑100 (MQ=255)
                                                                         tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaaca                                                                       >  2:496217/1‑150 (MQ=255)
                                                                           tatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaac                                                                     <  2:400591/150‑1 (MQ=255)
                                                                              ttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGcc                                                                  <  4:215163/150‑1 (MQ=255)
                                                                                  tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTg                                                              >  1:630065/1‑150 (MQ=255)
                                                                                   aaaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                                                             <  2:390313/149‑1 (MQ=255)
                                                                                    aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGc                                                            >  3:12387/1‑150 (MQ=255)
                                                                                     aTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCg                                                           <  4:117404/150‑1 (MQ=255)
                                                                                      tGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                                                          >  4:29123/1‑150 (MQ=255)
                                                                                        gTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAaaac                                                                                                                                                 <  2:385175/61‑1 (MQ=255)
                                                                                         tATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTAc                                                       >  2:124615/1‑150 (MQ=255)
                                                                                              aGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca                                                  <  4:227839/150‑1 (MQ=255)
                                                                                                  gCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa                                              >  4:262064/1‑150 (MQ=255)
                                                                                                   cAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAACCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAAc                                             <  2:301057/150‑1 (MQ=255)
                                                                                                    aaGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACt                                            >  1:167771/1‑150 (MQ=255)
                                                                                                          aaaTCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACcatc                                      >  1:522208/1‑150 (MQ=255)
                                                                                                             tCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                         >  1:68642/1‑112 (MQ=255)
                                                                                                             tCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaa                                                                         <  2:68642/112‑1 (MQ=255)
                                                                                                             tCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATg                                   >  2:305121/1‑150 (MQ=255)
                                                                                                               gAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGt                                 >  1:587906/1‑150 (MQ=255)
                                                                                                                         aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                          >  2:287543/1‑83 (MQ=255)
                                                                                                                         aTCAGTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCg                                                                                          <  1:287543/83‑1 (MQ=255)
                                                                                                                             gTTCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGaa                   >  3:235332/1‑150 (MQ=255)
                                                                                                                              ttCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATg                                   >  2:134163/1‑133 (MQ=255)
                                                                                                                              ttCGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATg                                   <  1:134163/133‑1 (MQ=255)
                                                                                                                                cGTGATCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATAc                >  2:277015/1‑150 (MQ=255)
                                                                                                                                    aTCCTGTTGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTcc            <  2:630065/150‑1 (MQ=255)
                                                                                                                                          ttGAAAAAAACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCaa      <  1:496217/150‑1 (MQ=255)
                                                                                                                                             aaaaaaACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAg   >  3:324383/1‑150 (MQ=255)
                                                                                                                                              aaaaaaCGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAgg  >  3:62432/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTACCATCATGGTAGCAGCGCGCAGAATACTTCCGCGCAACAGG  >  NC_000913_3_bme_pgi/4400286‑4400577

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