Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,020,957 T→C Y140H (TAC→CAC)  ygfM → putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,9570TC100.0% 1721.6 / NA 411Y140H (TAC→CAC) ygfMputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (191/220);  total (191/220)

Only 100 of 415 total aligned reads displayed.

TCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTCAGGCAGG  >  NC_000913/3020859‑3021054
                                                                                                  |                                                                                                 
 cgccgcCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGaccacc                                                                                                >  2:5802526/1‑101 (MQ=255)
  gccgccCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCt                                                                                               <  1:5801869/101‑1 (MQ=255)
     gcccgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGc                                                                                            <  1:5801310/101‑1 (MQ=255)
     gcccgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGc                                                                                            >  1:5802530/1‑101 (MQ=255)
      cccgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGcc                                                                                           <  1:5801786/101‑1 (MQ=255)
       ccgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCt                                                                                          <  2:5801857/101‑1 (MQ=255)
       ccgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCt                                                                                          <  2:5801996/101‑1 (MQ=255)
        cgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTg                                                                                         <  1:5802519/101‑1 (MQ=255)
        cgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTg                                                                                         <  1:5802213/101‑1 (MQ=255)
        cgccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTg                                                                                         >  2:5802536/1‑101 (MQ=255)
         gccAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGc                                                                                        >  1:5802537/1‑101 (MQ=255)
           cAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGccc                                                                                      <  1:5801667/101‑1 (MQ=255)
            aGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCa                                                                                     >  2:5802544/1‑101 (MQ=255)
             ggAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCAt                                                                                    <  1:5801849/101‑1 (MQ=255)
             ggAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCAt                                                                                    <  2:5801766/101‑1 (MQ=255)
             ggAAGAGGCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCAt                                                                                    <  2:5801649/101‑1 (MQ=255)
              gAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATg                                                                                   <  2:5802276/101‑1 (MQ=255)
               aaGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGc                                                                                  >  1:5802547/1‑101 (MQ=255)
                agagTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCg                                                                                 <  1:5801909/101‑1 (MQ=255)
                  agTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGAt                                                                               <  2:5802337/101‑1 (MQ=255)
                   gTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATc                                                                              <  1:5801959/101‑1 (MQ=255)
                    tCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCg                                                                             >  2:5802554/1‑101 (MQ=255)
                    tCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCg                                                                             >  1:5802552/1‑101 (MQ=255)
                         gctgctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGc                                                                        >  2:5802559/1‑101 (MQ=255)
                          ctgctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCt                                                                       >  1:5802564/1‑101 (MQ=255)
                           tgctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTa                                                                      <  1:5801940/101‑1 (MQ=255)
                            gctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTaa                                                                     >  2:5802566/1‑101 (MQ=255)
                             ctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAAc                                                                    <  1:5801754/101‑1 (MQ=255)
                             ctTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAAc                                                                    >  2:5802573/1‑101 (MQ=255)
                              ttCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAAcc                                                                   >  1:5802575/1‑101 (MQ=255)
                              ttCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAAcc                                                                   <  1:5802021/101‑1 (MQ=255)
                                 ccGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGaa                                                                >  1:5802582/1‑101 (MQ=255)
                                   gTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTg                                                                                         >  2:5802587/1‑74 (MQ=255)
                                    tCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAAtt                                                             <  2:5802247/101‑1 (MQ=255)
                                    tCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAAtt                                                             <  2:5801947/101‑1 (MQ=255)
                                     ccTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCCCCTGGCCTGCCCATGc                                                                                  >  2:5802593/1‑79 (MQ=255)
                                     ccTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTa                                                            >  1:5802592/1‑101 (MQ=255)
                                       tGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATTGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAAcc                                                                   >  2:5802599/1‑92 (MQ=255)
                                        gTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCa                                                         >  1:5802601/1‑101 (MQ=255)
                                          tGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCAtt                                                       >  2:5802606/1‑101 (MQ=255)
                                           ggCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTa                                                      <  2:5801584/101‑1 (MQ=255)
                                            gCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTaa                                                     <  2:5801976/101‑1 (MQ=255)
                                             cACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTaaa                                                    <  2:5801723/101‑1 (MQ=255)
                                              aCTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAg                                                   <  2:5801498/101‑1 (MQ=255)
                                               cTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGa                                                  <  1:5802069/101‑1 (MQ=255)
                                               cTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGa                                                  >  1:5802610/1‑101 (MQ=255)
                                                tGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGAt                                                 <  1:5801927/101‑1 (MQ=255)
                                                tGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGAt                                                 <  1:5801967/101‑1 (MQ=255)
                                                tGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGAt                                                 >  2:5802614/1‑101 (MQ=255)
                                                 ggATGCTTAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATc                                                <  2:5802066/101‑1 (MQ=255)
                                                 ggATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATc                                                <  2:5802123/101‑1 (MQ=255)
                                                  gATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATcc                                               <  1:5802093/101‑1 (MQ=255)
                                                   aTGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCg                                              >  1:5802620/1‑101 (MQ=255)
                                                    tGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGt                                             >  2:5802622/1‑101 (MQ=255)
                                                     gCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTa                                            <  1:5801987/101‑1 (MQ=255)
                                                     gCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTa                                            <  1:5802162/101‑1 (MQ=255)
                                                      cTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTAt                                           <  1:5801890/101‑1 (MQ=255)
                                                       tGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATc                                          >  2:5802628/1‑101 (MQ=255)
                                                        gAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCg                                         <  2:5801895/101‑1 (MQ=255)
                                                         aaCTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGc                                        >  2:5802632/1‑101 (MQ=255)
                                                           cTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCAc                                      <  1:5802121/101‑1 (MQ=255)
                                                           cTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCAc                                      >  1:5802635/1‑101 (MQ=255)
                                                             ggTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCt                                    <  1:5802153/101‑1 (MQ=255)
                                                             ggTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCt                                    <  2:5802026/101‑1 (MQ=255)
                                                             ggTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCt                                    <  2:5802045/101‑1 (MQ=255)
                                                             ggTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCt                                    <  2:5802219/101‑1 (MQ=255)
                                                              gTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTa                                            <  2:5802013/92‑1 (MQ=255)
                                                              gTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCtg                                   <  1:5801811/101‑1 (MQ=255)
                                                               tttttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCtgt                                  <  1:5802049/101‑1 (MQ=255)
                                                               tttttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCtgt                                  >  2:5802641/1‑101 (MQ=255)
                                                                ttttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCtgtg                                 <  1:5802239/101‑1 (MQ=255)
                                                                ttttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCtgtg                                 <  2:5801929/101‑1 (MQ=255)
                                                                 tttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGc                                >  1:5802644/1‑101 (MQ=255)
                                                                  ttGGCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCg                               >  1:5802648/1‑101 (MQ=255)
                                                                    ggCAACGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGAc                             >  2:5802651/1‑101 (MQ=255)
                                                                       aaCGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCg                          >  2:5802657/1‑101 (MQ=255)
                                                                       aaCGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCg                          <  2:5802226/101‑1 (MQ=255)
                                                                        aCGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGc                         <  1:5802037/101‑1 (MQ=255)
                                                                         cGGCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCa                        >  1:5802660/1‑101 (MQ=255)
                                                                          ggCGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCaa                       <  2:5802083/101‑1 (MQ=255)
                                                                            cGAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCaaaa                     >  2:5802665/1‑101 (MQ=255)
                                                                             gAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCaaaa                     >  2:5802677/1‑100 (MQ=255)
                                                                             gAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAAt                    >  1:5802669/1‑101 (MQ=255)
                                                                             gAAACGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAAt                    >  1:5802675/1‑101 (MQ=255)
                                                                               aaCGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTa                                            >  1:5802683/1‑75 (MQ=255)
                                                                                aCGCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAg                 >  2:5802685/1‑101 (MQ=255)
                                                                                  gCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGcc               <  2:5802156/101‑1 (MQ=255)
                                                                                  gCTGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGcc               >  1:5802690/1‑101 (MQ=255)
                                                                                    tGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGt             >  1:5802698/1‑101 (MQ=255)
                                                                                    tGTCAATCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGt             >  2:5802696/1‑101 (MQ=255)
                                                                                        aaTCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTtctc         <  1:5802328/101‑1 (MQ=255)
                                                                                        aaTCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTtctc         >  2:5802703/1‑101 (MQ=255)
                                                                                         aTCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTca        <  2:5802203/101‑1 (MQ=255)
                                                                                         aTCGAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTca        >  1:5802705/1‑101 (MQ=255)
                                                                                            gAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTcaggc     <  1:5802289/101‑1 (MQ=255)
                                                                                            gAGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTcaggc     >  1:5802710/1‑101 (MQ=255)
                                                                                             aGGACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTcaggca    <  2:5802179/101‑1 (MQ=255)
                                                                                              ggACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTcaggcag   <  1:5802003/101‑1 (MQ=255)
                                                                                              ggACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTcaggcag   >  2:5802713/1‑101 (MQ=255)
                                                                                               gACCACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGAc                             <  2:5802104/74‑1 (MQ=255)
                                                                                                  |                                                                                                 
TCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCTGCCCATGCGATCGCCTGCTAACCGAAATTATCATTAAAGATCCGTATCGCACCTGTGCGACCCGCAAAATTAGCCGTTCTCAGGCAGG  >  NC_000913/3020859‑3021054

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