Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,553,560 (T)6→5 coding (336/900 nt) hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,5550T.98.7% 314.0 / ‑2.0 76coding (331/900 nt)hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (44/31);  minor base G (1/0);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

AGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTG  >  NC_000913/2553427‑2553683
                                                                                                                                |                                                                                                                                
aGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCg                                                                                                                      >  1:1070438/1‑140 (MQ=255)
aGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCg                                                                                                                      >  2:9999/1‑140 (MQ=255)
 ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGa                                                                                                                     >  2:562314/1‑140 (MQ=255)
 ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGa                                                                                                                     >  2:413644/1‑140 (MQ=255)
 ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGa                                                                                                                     >  1:704911/1‑140 (MQ=255)
 ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGa                                                                                                                     >  1:793240/1‑140 (MQ=255)
  gcatcgCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGaa                                                                                                                    >  1:964196/6‑140 (MQ=255)
      tGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtttttat                                                                                                                            >  2:917288/1‑127 (MQ=255)
      tGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtttttat                                                                                                                            <  1:917288/129‑3 (MQ=255)
        ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCg                                                                                                              <  1:9999/140‑1 (MQ=255)
          tGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCggg                                                                                                            >  1:1208663/1‑140 (MQ=255)
            cTTCCGCCACCGCTCATCGCCCAGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTg                                                                                                          <  1:1606295/140‑1 (MQ=255)
                  ccACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATc                                                                                                    >  2:664786/1‑140 (MQ=255)
                  ccACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATc                                                                                                    <  1:1592248/140‑1 (MQ=255)
                  ccACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATCGCCGAAAAACCGGGTGCCGATc                                                                                                    >  1:452972/1‑140 (MQ=255)
                       gCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTc                                                                                               <  2:1208663/140‑1 (MQ=255)
                        cTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCt                                                                                              >  2:1270063/1‑140 (MQ=255)
                           aTCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtggt                                                                                           >  2:1599588/1‑140 (MQ=255)
                               cccGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGtt                                                                                       >  1:1156549/1‑140 (MQ=255)
                                ccGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGttt                                                                                      >  1:694567/1‑140 (MQ=255)
                                    aaCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCg                                                                                  <  1:185404/140‑1 (MQ=255)
                                    aaCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCg                                                                                  <  1:980321/140‑1 (MQ=255)
                                      cTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGt                                                                                >  2:1583193/1‑140 (MQ=255)
                                          ccGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCt                                                                            <  1:317944/140‑1 (MQ=255)
                                           cGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCtt                                                                           <  2:900760/140‑1 (MQ=255)
                                              gcgcAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGGATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAAGCGGGTGCAGATCCCGTCTGGTGGTTTGGCGGTGGATTCGa                                                                        >  2:1467498/1‑140 (MQ=255)
                                              gcgcAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGa                                                                        >  1:1212574/1‑140 (MQ=255)
                                              gcgcAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGa                                                                        >  2:833749/1‑140 (MQ=255)
                                              gcgcAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGa                                                                        >  1:1226672/1‑140 (MQ=255)
                                              gcgcAGTTTCGAGGCGATGAGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGa                                                                        >  2:316124/1‑140 (MQ=255)
                                                gcAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACt                                                                      <  2:694567/140‑1 (MQ=255)
                                                    tttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAc                                                                  >  2:1559362/1‑140 (MQ=255)
                                                     ttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAcc                                                                 <  2:1226672/140‑1 (MQ=255)
                                                     ttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAcc                                                                 <  2:704911/140‑1 (MQ=255)
                                                     ttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAcc                                                                 <  1:438421/140‑1 (MQ=255)
                                                     ttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAcc                                                                 <  1:967995/140‑1 (MQ=255)
                                                     ttCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAAcc                                                                 <  1:562314/140‑1 (MQ=255)
                                                        gAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCa                                                              <  2:1223526/140‑1 (MQ=255)
                                                             gATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTa                                                         <  2:1156549/140‑1 (MQ=255)
                                                                gggCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATgg                                                      >  2:135612/1‑140 (MQ=255)
                                                                gggCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGGGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATgg                                                      >  2:1030088/1‑140 (MQ=255)
                                                                         aCTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTgaaga                                             <  1:1599588/140‑1 (MQ=255)
                                                                         aCTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTgaaga                                             <  2:1212574/140‑1 (MQ=255)
                                                                           tGGTAGTGCATCCGCATAACCCGTATGTTCTCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGaagaag                                           >  1:978593/1‑140 (MQ=255)
                                                                           tGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTAt                                                        >  1:1647488/1‑127 (MQ=255)
                                                                           tGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTAt                                                        <  2:1647488/127‑1 (MQ=255)
                                                                           tGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGaagaag                                           <  1:197293/140‑1 (MQ=255)
                                                                             gTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTgaaga                                             <  2:156666/136‑1 (MQ=255)
                                                                             gTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTgaaga                                             >  1:156666/1‑136 (MQ=255)
                                                                              tAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATg                                        <  1:1583193/140‑1 (MQ=255)
                                                                              tAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATg                                        <  2:117966/140‑1 (MQ=255)
                                                                               aGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGc                                       >  1:861693/1‑140 (MQ=255)
                                                                                gTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAagaaga                                          <  1:1520550/136‑1 (MQ=255)
                                                                                gTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAagaaga                                          >  2:1520550/1‑136 (MQ=255)
                                                                                gTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCt                                      >  1:1289186/1‑140 (MQ=255)
                                                                                gTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCt                                      >  2:1124782/1‑140 (MQ=255)
                                                                                 tGCATCCGCATAACCCGTATGTTCCCACAAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTa                                     >  1:247285/1‑140 (MQ=255)
                                                                                   cATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTAtt                                   >  2:1194211/1‑140 (MQ=255)
                                                                                     tCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTa                                     >  1:53506/1‑136 (MQ=255)
                                                                                     tCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTa                                     <  2:53506/136‑1 (MQ=255)
                                                                                     tCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCa                                 >  2:229315/1‑140 (MQ=255)
                                                                                       cGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCAc                                >  2:906940/1‑139 (MQ=255)
                                                                                       cGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCAc                                <  1:906940/139‑1 (MQ=255)
                                                                                         cATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTgg                             >  2:1326674/1‑140 (MQ=255)
                                                                                           tAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCa                           >  2:202679/1‑140 (MQ=255)
                                                                                           tAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCa                           >  2:1008678/1‑140 (MQ=255)
                                                                                                 gTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     <  2:964196/140‑1 (MQ=255)
                                                                                                 gTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTATCGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     <  2:452972/140‑1 (MQ=255)
                                                                                                   aTGTTCCCACCAGCCACGCGAATGTGGGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTTAAGAAGATGCTATTCACTGGGATCGCAACg                   >  2:1018157/1‑140 (MQ=255)
                                                                                                     gTTCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGcc                 <  1:1194211/140‑1 (MQ=255)
                                                                                                       tCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     <  2:491910/134‑1 (MQ=255)
                                                                                                       tCCCACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     >  1:491910/1‑134 (MQ=255)
                                                                                                        cccACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     <  2:823747/133‑1 (MQ=255)
                                                                                                        cccACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCAc                     >  1:823747/1‑133 (MQ=255)
                                                                                                        cccACCAGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGt              >  2:1258834/1‑140 (MQ=255)
                                                                                                         ccaccaGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTg             >  1:724696/1‑140 (MQ=255)
                                                                                                            ccaGCCACGCGAATGTGCGGTGTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGAcc          >  2:1178332/1‑140 (MQ=255)
                                                                                                             caGCCACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCt         <  1:229315/140‑1 (MQ=255)
                                                                                                                ccACGCGAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCtgtg      >  2:872559/1‑140 (MQ=255)
                                                                                                                    gcgAATGTGCGGTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGtgcctg  >  1:715545/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                
AGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACTTAACCCCATTCTATGGTTTTGAAGAAGATGCTATTCACTGGCATCGCACCGCCCGTGACCTGTGCCTG  >  NC_000913/2553427‑2553683

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