Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 3,966,742 G→T 22.8% R109L (CGC→CTC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,966,7420GT22.8% 90.8 / 26.8 67R109L (CGC→CTC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (32/18);  new base T (8/7);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

GATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAA  >  NC_000913_3_bme_pgi/3966594‑3966888
                                                                                                                                                     |                                                                                                                                                  
gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTCATACCATCTTTGGTAAGATTCGCCCGCCGAAAAAAGGTA‑AACg                                                                                                                                                    >  1:157242/1‑149 (MQ=255)
    gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAACGTG‑ACCGCtct                                                                                                                                                >  1:441009/1‑147 (MQ=255)
     atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTAtt                                                                                                                                               >  2:262529/1‑149 (MQ=255)
     atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTGCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTAtt                                                                                                                                               <  1:863224/149‑1 (MQ=255)
      tttGGTTTCCTCCGTTCAGCAGACAGATCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGACAAATCCGCCGTTTCAACCTCCGCACTGGTGATCCCATTCCTGGAAAAATACGCCAGCCGAAAGAAGGGG‑AACCCAAttc                                                                                                                                              >  2:897394/1‑148 (MQ=255)
       ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCGTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCACGACAAGAAGGTG‑AAGGCTAtttt                                                                                                                                             >  1:792101/1‑149 (MQ=255)
        tGGTTGCCTCCTTTCCGCAACCCGCTCCTACCACGCCGTGCCTTATTACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTg                                                                                                                                            <  1:244372/149‑1 (MQ=255)
           tttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTGTCCCCTCGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTgggc                                                                                                                                         >  1:173164/1‑149 (MQ=255)
            ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCCCTGGTAAGATTCACCCGCCGAAAGAAGACG‑AAGACGATTTTGCgct                                                                                                                                        >  1:156210/1‑149 (MQ=255)
               ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCgctgct                                                                                                                                     >  1:561194/1‑149 (MQ=255)
               ctccGTTCCGCAGACAGCTCATACCTCGTCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCGCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCgctgct                                                                                                                                     <  1:386401/149‑1 (MQ=255)
                tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATAACATCTACGTTTCACCTAGCCAAATCCGCCGTTTCAACCTCCGCACAGCAGGACCCGTGTTGGGTAAGATTCGCCCGCCGCAGGAAGGTG‑AACGCTGTTTTACgctgcta                                                                                                                                    >  1:429109/1‑148 (MQ=255)
                       agcggACAGTGCCGACCTCGACGGTCCGGATGACATCTACGTTTCACCTAGCCAAAGCCGCTGTTTCAACCTACGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGGCGAAAGCAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTa                                                                                                                             <  2:517851/145‑1 (MQ=255)
                        gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTaa                                                                                                                            >  1:411998/1‑149 (MQ=255)
                           gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                                                                                         <  1:694635/149‑1 (MQ=255)
                                  cctccctcGCCGGTCCTTATTACATTTTCGTTTCCCTTATTCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                                                                                  <  1:140090/149‑1 (MQ=255)
                                   cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTTTCAGTACTTACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑TACGTTATTTTGGGCTGCTGAAAGTTAACGAAGTTAACt                                                                                                                 <  1:429977/149‑1 (MQ=255)
                                       cTCGCCGGTCCTGATGACATATACGGTTCCCCGAACGAAATCCGCAGTTGCAACATCCGCACTGGTGAGACCAGCTCTGGTAAGATTCGCCCGCCGAAAGCAGTTG‑AACGCAATTTTGATCTACTGCAAGTTACCAAAGTTAACTTctc                                                                                                             >  2:843085/1‑147 (MQ=255)
                                          gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                                                                                          >  1:476683/1‑149 (MQ=255)
                                              gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                                                                                      >  1:50798/1‑149 (MQ=255)
                                                  tgatgaCATCTACGCTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGCTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGCTAACTCCGACAAACCTGaaa                                                                                                  >  2:628451/1‑149 (MQ=255)
                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAGGGTG‑AACGCTATTTTGAGCGGCTGAAAGTTAACGAAGTTAACTTCGACAGACCTGacaa                                                                                                 >  2:495237/1‑149 (MQ=255)
                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                                                                                                 <  1:463157/149‑1 (MQ=255)
                                                     tgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATGCCCTCTCTGGTAAGATTCGCCCGCCGAAAGACGGTG‑AAGGCTATTTGGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCGGaaaaag                                                                                               >  2:726729/1‑147 (MQ=255)
                                                         aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                                           >  2:471966/1‑149 (MQ=255)
                                                          tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTTTAGTAGGATTTGCCTTCTGAAA‑ATGGAGTAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAATACTGAAAACGCCCga                                                                                          >  2:676913/1‑148 (MQ=255)
                                                            tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGATGCTGAAAGTTAACGAAGTTAACTTCGACAACCCTGAAAAAGccaccac                                                                                        >  2:48309/1‑144 (MQ=255)
                                                                tttCCCCTAGCCAAATCCGGCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTATGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaccaac                                                                                    <  1:854402/149‑2 (MQ=255)
                                                                  tCCCCTAGCCAAATCCGGCGTTTCAACCTCCTCCCTGGTGATACCATCTCTTGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaagg                                                                                  <  2:855573/149‑3 (MQ=255)
                                                                     ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATCCGCCCGCCGAAAGAAGGTA‑AACGCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCCACAAACCTGCAAACGCCCGCAACACACTCct                                                                               >  2:727731/1‑149 (MQ=255)
                                                                      cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                              >  1:580467/1‑149 (MQ=255)
                                                                      cTAGCCAAATCCGCCGTTTAAACCTCCGCACTGGTGATACCATATCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCGCTTTTGCGCTGCTGAAAGTTAACGACGTTAACGTCGACAAACCGGAAGACGCCCGCAACAAAATCctc                                                                              >  2:861527/1‑149 (MQ=255)
                                                                       tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct                                                                             >  2:609808/1‑149 (MQ=255)
                                                                          ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                          >  1:159508/1‑149 (MQ=255)
                                                                          ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                          >  1:849953/1‑149 (MQ=255)
                                                                           cAAACCCGCCGTTTCAACCTCCGAACAGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga                                                                         <  2:827044/149‑1 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTGACTTCGACAAACCTGAAAACGTCCGCAACAAAATCCTCTTTGAGAAc                                                                     >  1:116790/1‑149 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                     <  1:539150/149‑1 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                     >  1:493120/1‑149 (MQ=255)
                                                                                ccgccgTTTGATAGTCCGCGTGGGTGCTACCATCTCTGGTAAGATTGGCCGGGCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCGGAAAGTTAACGAAGGTAACTTCGACAAACCTGAAAACGGCCGCAACAAAATCCTCTGTGAGAACt                                                                    <  2:868138/149‑1 (MQ=255)
                                                                                  gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCAGCAACAAAATACTCTTTAAGAACTTa                                                                  >  1:743630/1‑149 (MQ=255)
                                                                                      tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACAAATTTAACTACGACAAACCTGACAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAcac                                                              >  2:738806/1‑147 (MQ=255)
                                                                                          aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACGACGACAAACCGAAAAACGCCCGCAACAAAAACCTCTTTGAGAACTTAACCCCGCt                                                          >  1:718962/1‑149 (MQ=255)
                                                                                            ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAGGTTAACGAAGTTAACTTCGACAAACCCAAAAACGCCCGCAACAAAATCCCCTTTGGGAACTTATCCCCGCGGc                                                        >  2:186804/1‑149 (MQ=255)
                                                                                             cTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGTG‑AACTCTATTTTGAGCTGCTGAAAGTCAACGAAGTTAACTTCGACAAACCTGAAAACGCCAGCAACAAAATCCTCTTTGAGATCTGAACCCCGCTGCa                                                       >  1:79315/1‑149 (MQ=255)
                                                                                                          aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                                          <  2:561194/149‑1 (MQ=255)
                                                                                                           tACCATTTCTGGTAAGATTCGCCGGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                         <  2:131270/149‑1 (MQ=255)
                                                                                                            aCCATTCCTGGTAAGATTCGCCCGCCTAAAGAAGGTG‑AACTCTATTTTGAGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                        <  1:135525/149‑1 (MQ=255)
                                                                                                            aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGATG‑AACGTTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                        <  1:135689/149‑1 (MQ=255)
                                                                                                                    cggTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTAAGGGGTAACGGAGTTAACTTCGCCAAACCTGAAAATGCCAGCAACGAAATCTTCTTTGAGAACTTAACCCAGCTGCACGCAACATCTCTTTTGCGTATgg                                >  2:680707/2‑149 (MQ=255)
                                                                                                                      gTAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa                              <  1:471966/149‑1 (MQ=255)
                                                                                                                      gTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAATAAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGCGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCATATGGaa                              >  1:569241/1‑149 (MQ=255)
                                                                                                                       tAAGATTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc                             >  1:275684/1‑149 (MQ=255)
                                                                                                                          gATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTg                          >  2:572509/1‑149 (MQ=255)
                                                                                                                           aTTCGCCCTACGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                         <  1:262529/149‑1 (MQ=255)
                                                                                                                           aTTCGCCCGCCGAAAGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                         <  1:613573/149‑1 (MQ=255)
                                                                                                                                 ccgccgAAAGAAGGTG‑AACTCTATTTTTCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg                   <  1:103708/149‑1 (MQ=255)
                                                                                                                                   gccgACAGAAGGTG‑AACGCTATTTTGCGCAGCTGAAAGTTAACGAAGTTAACTTCGACAACCCTGAAAACGCCCGCCCGAAAATCACCTTTGAGAGCTTAACCCCGCTTCACGCAAACACACGTCTGCGTATGGAAGGTGGTAACgggt                 >  2:779583/1‑147 (MQ=255)
                                                                                                                                   gccgAAAGAAGGTG‑AACTCTATTTTGCGATGGTGAATGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAAAAATTAACACCGATGACCACAAACTCTCGTATGCGTATGGAAAGTCGTAACCGtt                 >  2:503331/1‑149 (MQ=255)
                                                                                                                                       aaaGAAGGTG‑AACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGTGAACGTGTCTAc             >  2:623897/1‑149 (MQ=255)
                                                                                                                                        aagaagGTT‑AACGTTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAGCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt            <  2:550414/149‑1 (MQ=255)
                                                                                                                                           aagGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAATCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa         <  1:824702/149‑1 (MQ=255)
                                                                                                                                            agGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTAAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAg        >  1:632014/1‑149 (MQ=255)
                                                                                                                                              gTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGACCGGGGTAACGGTTCTACTAAAGAt      >  2:251205/1‑149 (MQ=255)
                                                                                                                                              gTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt      >  1:536797/1‑149 (MQ=255)
                                                                                                                                                   aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTaa  <  1:556758/149‑1 (MQ=255)
                                                                                                                                                   aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTaa  <  1:649296/149‑1 (MQ=255)
                                                                                                                                                   aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTCCGTATGGAACGTGGTAACGGTTCTACTAAAGATTaaa  >  1:397862/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTG‑AACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAA  >  NC_000913_3_bme_pgi/3966594‑3966888

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