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

GATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAA  >  NC_000913_3_bme_pgi/3966594‑3966888
                                                                                                                                                     |                                                                                                                                                  
gatggatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTCATACCATCTTTGGTAAGATTCGCCCGCCGAAAAAAGGTAAACg                                                                                                                                                    >  1:157242/1‑149 (MQ=255)
    gatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAACGTGACCGCtct                                                                                                                                                >  1:441009/1‑147 (MQ=255)
     atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                                                                                               >  2:262529/1‑149 (MQ=255)
     atTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTGCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                                                                                               <  1:863224/149‑1 (MQ=255)
      tttGGTTTCCTCCGTTCAGCAGACAGATCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGACAAATCCGCCGTTTCAACCTCCGCACTGGTGATCCCATTCCTGGAAAAATACGCCAGCCGAAAGAAGGGGAACCCAAttc                                                                                                                                              >  2:897394/1‑148 (MQ=255)
       ttGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCGTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCACGACAAGAAGGTGAAGGCTAtttt                                                                                                                                             >  1:792101/1‑149 (MQ=255)
        tGGTTGCCTCCTTTCCGCAACCCGCTCCTACCACGCCGTGCCTTATTACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                                                                                            <  1:244372/149‑1 (MQ=255)
           tttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTGTCCCCTCGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTgggc                                                                                                                                         >  1:173164/1‑149 (MQ=255)
            ttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCCCTGGTAAGATTCACCCGCCGAAAGAAGACGAAGACGATTTTGCgct                                                                                                                                        >  1:156210/1‑149 (MQ=255)
               ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                                                                                                     >  1:561194/1‑149 (MQ=255)
               ctccGTTCCGCAGACAGCTCATACCTCGTCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCGCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCgctgct                                                                                                                                     <  1:386401/149‑1 (MQ=255)
                tccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATAACATCTACGTTTCACCTAGCCAAATCCGCCGTTTCAACCTCCGCACAGCAGGACCCGTGTTGGGTAAGATTCGCCCGCCGCAGGAAGGTGAACGCTGTTTTACgctgcta                                                                                                                                    >  1:429109/1‑148 (MQ=255)
                       agcggACAGTGCCGACCTCGACGGTCCGGATGACATCTACGTTTCACCTAGCCAAAGCCGCTGTTTCAACCTACGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGGCGAAAGCAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                                                                                                                             <  2:517851/145‑1 (MQ=255)
                        gCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTaa                                                                                                                            >  1:411998/1‑149 (MQ=255)
                           gACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                                                                                         <  1:694635/149‑1 (MQ=255)
                                  cctccctcGCCGGTCCTTATTACATTTTCGTTTCCCTTATTCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCGGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                                                                                                  <  1:140090/149‑1 (MQ=255)
                                   cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTTTCAGTACTTACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGTACGTTATTTTGGGCTGCTGAAAGTTAACGAAGTTAACt                                                                                                                 <  1:429977/149‑1 (MQ=255)
                                       cTCGCCGGTCCTGATGACATATACGGTTCCCCGAACGAAATCCGCAGTTGCAACATCCGCACTGGTGAGACCAGCTCTGGTAAGATTCGCCCGCCGAAAGCAGTTGAACGCAATTTTGATCTACTGCAAGTTACCAAAGTTAACTTctc                                                                                                             >  2:843085/1‑147 (MQ=255)
                                          gCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                                                                                          >  1:476683/1‑149 (MQ=255)
                                              gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                                                                                                      >  1:50798/1‑149 (MQ=255)
                                                  tgatgaCATCTACGCTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGCTGATACCATCTCTGGTAAGCTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGCTAACTCCGACAAACCTGaaa                                                                                                  >  2:628451/1‑149 (MQ=255)
                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAGGGTGAACGCTATTTTGAGCGGCTGAAAGTTAACGAAGTTAACTTCGACAGACCTGacaa                                                                                                 >  2:495237/1‑149 (MQ=255)
                                                   gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                                                                                                 <  1:463157/149‑1 (MQ=255)
                                                     tgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATGCCCTCTCTGGTAAGATTCGCCCGCCGAAAGACGGTGAAGGCTATTTGGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCGGaaaaag                                                                                               >  2:726729/1‑147 (MQ=255)
                                                         aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                                                                                           >  2:471966/1‑149 (MQ=255)
                                                          tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTTTAGTAGGATTTGCCTTCTGAAAATGGAGTAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAATACTGAAAACGCCCga                                                                                          >  2:676913/1‑148 (MQ=255)
                                                            tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGATGCTGAAAGTTAACGAAGTTAACTTCGACAACCCTGAAAAAGccaccac                                                                                        >  2:48309/1‑144 (MQ=255)
                                                                tttCCCCTAGCCAAATCCGGCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTATGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaccaac                                                                                    <  1:854402/149‑2 (MQ=255)
                                                                  tCCCCTAGCCAAATCCGGCGTTTCAACCTCCTCCCTGGTGATACCATCTCTTGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaagg                                                                                  <  2:855573/149‑3 (MQ=255)
                                                                     ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATCCGCCCGCCGAAAGAAGGTAAACGCTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCCACAAACCTGCAAACGCCCGCAACACACTCct                                                                               >  2:727731/1‑149 (MQ=255)
                                                                      cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc                                                                              >  1:580467/1‑149 (MQ=255)
                                                                      cTAGCCAAATCCGCCGTTTAAACCTCCGCACTGGTGATACCATATCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCGCTTTTGCGCTGCTGAAAGTTAACGACGTTAACGTCGACAAACCGGAAGACGCCCGCAACAAAATCctc                                                                              >  2:861527/1‑149 (MQ=255)
                                                                       tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct                                                                             >  2:609808/1‑149 (MQ=255)
                                                                          ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                          >  1:159508/1‑149 (MQ=255)
                                                                          ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg                                                                          >  1:849953/1‑149 (MQ=255)
                                                                           cAAACCCGCCGTTTCAACCTCCGAACAGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga                                                                         <  2:827044/149‑1 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTGACTTCGACAAACCTGAAAACGTCCGCAACAAAATCCTCTTTGAGAAc                                                                     >  1:116790/1‑149 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                     <  1:539150/149‑1 (MQ=255)
                                                                               tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                     >  1:493120/1‑149 (MQ=255)
                                                                                ccgccgTTTGATAGTCCGCGTGGGTGCTACCATCTCTGGTAAGATTGGCCGGGCGAAAGAAGGTGAACTCTATTTTGCGCTGCGGAAAGTTAACGAAGGTAACTTCGACAAACCTGAAAACGGCCGCAACAAAATCCTCTGTGAGAACt                                                                    <  2:868138/149‑1 (MQ=255)
                                                                                  gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCAGCAACAAAATACTCTTTAAGAACTTa                                                                  >  1:743630/1‑149 (MQ=255)
                                                                                      tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACAAATTTAACTACGACAAACCTGACAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAcac                                                              >  2:738806/1‑147 (MQ=255)
                                                                                          aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACGACGACAAACCGAAAAACGCCCGCAACAAAAACCTCTTTGAGAACTTAACCCCGCt                                                          >  1:718962/1‑149 (MQ=255)
                                                                                            ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAGGTTAACGAAGTTAACTTCGACAAACCCAAAAACGCCCGCAACAAAATCCCCTTTGGGAACTTATCCCCGCGGc                                                        >  2:186804/1‑149 (MQ=255)
                                                                                             cTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGACGGTGAACTCTATTTTGAGCTGCTGAAAGTCAACGAAGTTAACTTCGACAAACCTGAAAACGCCAGCAACAAAATCCTCTTTGAGATCTGAACCCCGCTGCa                                                       >  1:79315/1‑149 (MQ=255)
                                                                                                          aTACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGt                                          <  2:561194/149‑1 (MQ=255)
                                                                                                           tACCATTTCTGGTAAGATTCGCCGGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTc                                         <  2:131270/149‑1 (MQ=255)
                                                                                                            aCCATTCCTGGTAAGATTCGCCCGCCTAAAGAAGGTGAACTCTATTTTGAGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                        <  1:135525/149‑1 (MQ=255)
                                                                                                            aCCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGATGAACGTTATTTTGCGCTGCTGAAAGGTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCt                                        <  1:135689/149‑1 (MQ=255)
                                                                                                                    cggTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTAAGGGGTAACGGAGTTAACTTCGCCAAACCTGAAAATGCCAGCAACGAAATCTTCTTTGAGAACTTAACCCAGCTGCACGCAACATCTCTTTTGCGTATgg                                >  2:680707/2‑149 (MQ=255)
                                                                                                                      gTAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGaa                              <  1:471966/149‑1 (MQ=255)
                                                                                                                      gTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAATAAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGCGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCATATGGaa                              >  1:569241/1‑149 (MQ=255)
                                                                                                                       tAAGATTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAAc                             >  1:275684/1‑149 (MQ=255)
                                                                                                                          gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTg                          >  2:572509/1‑149 (MQ=255)
                                                                                                                           aTTCGCCCTACGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                         <  1:262529/149‑1 (MQ=255)
                                                                                                                           aTTCGCCCGCCGAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                         <  1:613573/149‑1 (MQ=255)
                                                                                                                                 ccgccgAAAGAAGGTGAACTCTATTTTTCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg                   <  1:103708/149‑1 (MQ=255)
                                                                                                                                   gccgACAGAAGGTGAACGCTATTTTGCGCAGCTGAAAGTTAACGAAGTTAACTTCGACAACCCTGAAAACGCCCGCCCGAAAATCACCTTTGAGAGCTTAACCCCGCTTCACGCAAACACACGTCTGCGTATGGAAGGTGGTAACgggt                 >  2:779583/1‑147 (MQ=255)
                                                                                                                                   gccgAAAGAAGGTGAACTCTATTTTGCGATGGTGAATGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAAAAATTAACACCGATGACCACAAACTCTCGTATGCGTATGGAAAGTCGTAACCGtt                 >  2:503331/1‑149 (MQ=255)
                                                                                                                                       aaaGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGTGAACGTGTCTAc             >  2:623897/1‑149 (MQ=255)
                                                                                                                                        aagaagGTTAACGTTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAGCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt            <  2:550414/149‑1 (MQ=255)
                                                                                                                                           aagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAATCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGaa         <  1:824702/149‑1 (MQ=255)
                                                                                                                                            agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTAAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAg        >  1:632014/1‑149 (MQ=255)
                                                                                                                                              gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGACCGGGGTAACGGTTCTACTAAAGAt      >  2:251205/1‑149 (MQ=255)
                                                                                                                                              gTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGAt      >  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)
                                                                                                                                                     |                                                                                                                                                  
GATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAA  >  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: