Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,915,493 G→A I14I (ATC→ATT proQ ← RNA chaperone ProQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,915,4930GA100.0% 1513.8 / NA 347I14I (ATC→ATTproQRNA chaperone ProQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (185/162);  total (185/162)

Only 100 of 350 total aligned reads displayed.

GCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCA  >  NC_000913/1915395‑1915590
                                                                                                  |                                                                                                 
gCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAAt                                                                                                 >  1:3506191/1‑101 (MQ=255)
 cAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACGGTGGGGAAAACGTTCGGCCAGAAACGCAAtt                                                                                                >  1:3506194/1‑101 (MQ=255)
 cAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAAtt                                                                                                >  2:3506193/1‑101 (MQ=255)
   acacGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAAtt                                                                                                <  2:3505493/99‑1 (MQ=255)
    cacGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTAct                                                                                             <  2:3505747/101‑1 (MQ=255)
    cacGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTAct                                                                                             >  1:3506201/1‑101 (MQ=255)
       gATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtct                                                                                          <  2:3505656/101‑1 (MQ=255)
       gATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtct                                                                                          <  2:3506009/101‑1 (MQ=255)
        aTCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtctt                                                                                         >  2:3506208/1‑101 (MQ=255)
        aTCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtctt                                                                                         >  1:3506207/1‑101 (MQ=255)
          cGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTa                                                                                       <  2:3505699/101‑1 (MQ=255)
           gACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTAc                                                                                      >  1:3506213/1‑101 (MQ=255)
            aCCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACt                                                                                     >  2:3506216/1‑101 (MQ=255)
            aCCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACt                                                                                     <  1:3505885/101‑1 (MQ=255)
             ccAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtctt                                                                                         >  1:3506225/1‑96 (MQ=255)
             ccAAATCCTGAAAAATACCGATTTTCAGCGGACGCGATTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTg                                                                                    >  2:3506222/1‑101 (MQ=255)
              cAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGc                                                                                   <  1:3505719/101‑1 (MQ=255)
               aaaTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCt                                                                                  <  1:3505860/101‑1 (MQ=255)
               aaaTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCt                                                                                  >  2:3506228/1‑101 (MQ=255)
                aaTCCTGAAAAATACCGATTTTCAGCGGACGCTCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTa                                                                                 <  2:3505586/101‑1 (MQ=255)
                aaTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTa                                                                                 <  1:3505872/101‑1 (MQ=255)
                 aTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTActtctt                                                                                         <  1:3505668/92‑1 (MQ=255)
                 aTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTAt                                                                                >  1:3506233/1‑101 (MQ=255)
                 aTCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTAt                                                                                <  1:3505791/101‑1 (MQ=255)
                   ccTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTc                                                                              <  1:3505626/101‑1 (MQ=255)
                       aaaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACt                                                                          >  1:3506239/1‑101 (MQ=255)
                        aaaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACtt                                                                         <  1:3505777/101‑1 (MQ=255)
                         aaaTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTa                                                                        >  2:3506243/1‑101 (MQ=255)
                           aTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAgg                                                                      >  2:3506249/1‑101 (MQ=255)
                           aTACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAgg                                                                      >  1:3506248/1‑101 (MQ=255)
                            tACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGt                                                                     >  1:3506252/1‑101 (MQ=255)
                             aCCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGtt                                                                    <  1:3505897/101‑1 (MQ=255)
                                 aTTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGAtt                                                                >  2:3506260/1‑101 (MQ=255)
                                    ttCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTc                                                             >  1:3506261/1‑101 (MQ=255)
                                     tCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTcc                                                            <  1:3505650/101‑1 (MQ=255)
                                      cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCa                                                           >  1:3506264/1‑101 (MQ=255)
                                      cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCa                                                           <  2:3506107/101‑1 (MQ=255)
                                      cAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCa                                                           <  1:3505832/101‑1 (MQ=255)
                                       aGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCAt                                                          <  1:3505941/101‑1 (MQ=255)
                                       aGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGGTGATTTTCCAt                                                          >  2:3506268/1‑101 (MQ=255)
                                        gCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATg                                                         >  2:3506272/1‑101 (MQ=255)
                                          ggACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGaa                                                       <  2:3505616/101‑1 (MQ=255)
                                          ggACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGaa                                                       >  1:3506281/1‑101 (MQ=255)
                                          ggACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGaa                                                       >  2:3506276/1‑101 (MQ=255)
                                           gACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGaaa                                                      >  1:3506285/1‑101 (MQ=255)
                                           gACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGaaa                                                      >  2:3506284/1‑101 (MQ=255)
                                            aCGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAAt                                                     >  2:3506288/1‑101 (MQ=255)
                                               cgcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTa                                                                        <  2:3505766/79‑1 (MQ=255)
                                               cgcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTc                                                  <  2:3505816/101‑1 (MQ=255)
                                                gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTcc                                                 >  2:3506293/1‑101 (MQ=255)
                                                gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTcc                                                 >  1:3506298/1‑101 (MQ=255)
                                                gcTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTcc                                                 >  2:3506299/1‑101 (MQ=255)
                                                  ttCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTg                                               >  1:3506303/1‑101 (MQ=255)
                                                       cTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTAc                                          <  2:3505913/101‑1 (MQ=255)
                                                        ttCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACa                                         <  1:3505758/101‑1 (MQ=255)
                                                         tCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACaa                                        >  1:3506310/1‑101 (MQ=255)
                                                          ccGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAAc                                       <  2:3505802/101‑1 (MQ=255)
                                                          ccGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAAc                                       <  2:3505962/101‑1 (MQ=255)
                                                           cGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACg                                      >  2:3506311/1‑101 (MQ=255)
                                                            gCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTa                                           <  1:3505449/95‑1 (MQ=255)
                                                            gCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACgg                                     <  1:3505906/101‑1 (MQ=255)
                                                              aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                   >  1:3506316/1‑101 (MQ=255)
                                                              aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                   <  2:3505887/101‑1 (MQ=255)
                                                              aCTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGAc                                   >  2:3506317/1‑101 (MQ=255)
                                                                 gAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGAtttt                                                              >  1:3506325/1‑71 (MQ=255)
                                                                 gAAACAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTa                                >  2:3506324/1‑101 (MQ=255)
                                                                    aCAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGcc                             >  1:3506328/1‑101 (MQ=255)
                                                                     cAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAg                               <  1:3506010/98‑1 (MQ=255)
                                                                     cAGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCca                            <  2:3505907/101‑1 (MQ=255)
                                                                      aGTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaa                           >  1:3506333/1‑101 (MQ=255)
                                                                       gTGGGGAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaac                          >  1:3506336/1‑101 (MQ=255)
                                                                         ggggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaacaa                        <  2:3505735/101‑1 (MQ=255)
                                                                         ggggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCcaacaa                        >  2:3506339/1‑101 (MQ=255)
                                                                          gtgAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTc                                                             <  1:3505801/61‑1 (MQ=255)
                                                                          gagAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTg                                               <  1:3505587/75‑1 (MQ=255)
                                                                           ggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGc                              <  2:3505826/93‑1 (MQ=255)
                                                                           ggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgc                      <  2:3505858/101‑1 (MQ=255)
                                                                           ggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgc                      <  1:3505501/101‑1 (MQ=255)
                                                                           ggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgc                      <  1:3505733/101‑1 (MQ=255)
                                                                           ggAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgc                      <  2:3505715/101‑1 (MQ=255)
                                                                            gAAAACGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAgcg                     >  1:3506343/1‑101 (MQ=255)
                                                                              aaaCGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCa                   >  2:3506346/1‑101 (MQ=255)
                                                                               aaCGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAg                  >  1:3506351/1‑101 (MQ=255)
                                                                                aCGTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAgg                 >  2:3506353/1‑101 (MQ=255)
                                                                                  gTTCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCa               >  1:3506358/1‑101 (MQ=255)
                                                                                    tCGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATg             >  2:3506362/1‑101 (MQ=255)
                                                                                     cGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGa            <  1:3505991/101‑1 (MQ=255)
                                                                                     cGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGa            >  1:3506363/1‑101 (MQ=255)
                                                                                     cGGCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGa            <  2:3505946/101‑1 (MQ=255)
                                                                                      ggCCAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGaa           >  2:3506366/1‑101 (MQ=255)
                                                                                         cAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg        >  1:3506371/1‑101 (MQ=255)
                                                                                         cAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg        <  1:3506029/101‑1 (MQ=255)
                                                                                         cAGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAg        <  1:3505820/101‑1 (MQ=255)
                                                                                          aGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCAt              <  2:3505673/94‑1 (MQ=255)
                                                                                          aGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAgg       >  1:3506378/1‑101 (MQ=255)
                                                                                          aGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAgg       <  1:3505970/101‑1 (MQ=255)
                                                                                          aGAAACGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAgg       >  2:3506376/1‑101 (MQ=255)
                                                                                             aaCGCAATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGt    >  2:3506383/1‑101 (MQ=255)
                                                                                               cgcaATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCa  <  2:3505933/101‑1 (MQ=255)
                                                                                               cgcaATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCa  >  1:3506387/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GCAACACGATCGACCAAATCCTGAAAAATACCGATTTTCAGCGGACGCGCTTCACCTTCCGCACTGAAACAGTGGGGAAAACGTTCGGCCAGAAACGCGATTACTTCTTTACTGCTATTCAACTTAGGTTGATTTTCCATGAAATTTCCTGATTACAACGGACGTAGCCAACAAGCGCAGGCATGAACAGGCGTCA  >  NC_000913/1915395‑1915590

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: