Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,904,592 C→A 21.4% T233N (ACT→AAT)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,5920CA21.4% 288.9 / 72.3 126T233N (ACT→AAT) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base C (46/53);  new base A (15/12);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 128 total aligned reads displayed.

TGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGT  >  NC_000913/1904493‑1904691
                                                                                                    |                                                                                                   
  ggggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTa                                                                                                   <  2:231949/100‑1 (MQ=255)
    ggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAtt                                                                                                 >  1:309801/1‑100 (MQ=255)
        ttGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTgg                                                                                             <  1:806252/100‑1 (MQ=255)
         tGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGa                                                                                            <  1:3377933/100‑1 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                           <  2:968692/100‑1 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                           <  2:265989/100‑1 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                           >  2:1929008/1‑100 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                           >  2:1451931/1‑100 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                           >  2:1109673/1‑100 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAc                                                                                           <  1:3729058/100‑1 (MQ=255)
          ggTTAACAAGTGGACACATGTCAACAGCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGa                                                                                            <  1:3922781/99‑1 (MQ=255)
            ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCa                                                                                         >  1:1287180/1‑100 (MQ=255)
             tAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAg                                                                                        <  2:628408/100‑1 (MQ=255)
              aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGt                                                                                       >  2:2265528/1‑100 (MQ=255)
              aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagta                                                                                       <  1:3030032/100‑5 (MQ=255)
              aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagta                                                                                       <  2:2361136/100‑5 (MQ=255)
               aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagta                                                                                       <  1:733076/100‑5 (MQ=38)
               aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                      <  1:1805640/100‑6 (MQ=255)
               aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                      <  2:837701/100‑6 (MQ=255)
                  aGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAAt                                                                                   <  2:768860/100‑1 (MQ=255)
                   gTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtagtaac                                                                                  <  1:2466319/100‑10 (MQ=38)
                     ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagtagtaacg                                                                                 <  1:1522143/100‑11 (MQ=37)
                     ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGcc                                                                                >  2:1501381/1‑100 (MQ=255)
                        cacaTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAg                                                                              >  2:2791911/1‑99 (MQ=255)
                         acaTGTCGACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGcc                                                                                <  2:2711133/96‑1 (MQ=255)
                           aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCTAGGc                                                                            >  2:3482157/1‑98 (MQ=255)
                           aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCt                                                                          <  1:1945542/100‑1 (MQ=255)
                             gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTgg                                                                        <  2:3223724/100‑1 (MQ=255)
                               cAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                      <  2:2800965/100‑1 (MQ=255)
                                 aCATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTAcc                                                                    >  2:1732009/1‑100 (MQ=255)
                                   aTCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCAc                                                                  <  2:2148118/100‑1 (MQ=255)
                                    tCCCGCTGGTTGTCTCTCGCATTACTTCCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTCTTTTTGACCCGTTAATGCCAAGCCTGGTACCAct                                                                 >  2:1905069/1‑100 (MQ=255)
                                     cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCAGGCCTGGTACCActgc                                                               <  2:2072504/100‑1 (MQ=255)
                                      ccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgc                                                               >  1:3205516/1‑100 (MQ=255)
                                      ccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgc                                                               >  2:1252411/1‑100 (MQ=255)
                                        gCTGGTTGTCTCTCGCATTACTGACTATACGGGCAAATAAAATGTTACTACTGTCCAGACTATTCTGGACTAGTTAATGCAAGGCCTGGTACCAct                                                                 <  2:1866872/96‑1 (MQ=255)
                                        gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                             >  1:3721150/1‑100 (MQ=255)
                                        gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                             <  1:639045/100‑1 (MQ=255)
                                         cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                            <  1:1829235/100‑1 (MQ=255)
                                         cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                            <  2:2883703/100‑1 (MQ=255)
                                         cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                            <  2:1568589/100‑1 (MQ=255)
                                           ggTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgctg                                                          >  2:3107421/1‑100 (MQ=255)
                                             ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAc                                                        >  1:1807829/1‑100 (MQ=255)
                                              tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                       <  1:1731879/100‑1 (MQ=255)
                                              tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                       <  1:1731867/100‑1 (MQ=255)
                                              tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                       <  1:2749846/100‑1 (MQ=255)
                                                tctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                       <  1:436042/98‑1 (MQ=255)
                                                 ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                    >  1:1208926/1‑100 (MQ=255)
                                                 ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                    >  2:1452779/1‑100 (MQ=255)
                                                   ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                  <  1:2352757/100‑1 (MQ=255)
                                                   ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                  >  2:568608/1‑100 (MQ=255)
                                                    tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                 <  1:2848294/100‑1 (MQ=255)
                                                     cGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                 <  1:1318200/99‑1 (MQ=255)
                                                      gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAAGGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                               <  1:2957941/100‑1 (MQ=255)
                                                            cTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgtg                                         >  1:1597789/1‑100 (MQ=255)
                                                              gACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                       >  1:2002693/1‑100 (MQ=255)
                                                              gACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATGCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                        >  2:2678037/1‑99 (MQ=255)
                                                               aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                      >  1:2212172/1‑100 (MQ=255)
                                                               aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                      >  2:759045/1‑100 (MQ=255)
                                                                 cAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAc                                    >  1:2089430/1‑100 (MQ=255)
                                                                  aGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACt                                   <  1:987561/100‑1 (MQ=255)
                                                                   gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAc                                    <  2:1117839/98‑1 (MQ=255)
                                                                   gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                  >  2:165572/1‑100 (MQ=255)
                                                                     cGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                >  1:2519404/1‑100 (MQ=255)
                                                                     cGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                >  2:905718/1‑100 (MQ=255)
                                                                      gggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcgc                               <  2:807548/100‑1 (MQ=255)
                                                                       ggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCa                              >  1:3259251/1‑100 (MQ=255)
                                                                        gCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                             >  1:2404348/1‑100 (MQ=255)
                                                                        gCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                             >  2:2666367/1‑100 (MQ=255)
                                                                           aaGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaa                          >  1:1646461/1‑100 (MQ=255)
                                                                             gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  1:1746442/100‑1 (MQ=255)
                                                                             gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  1:2605649/100‑1 (MQ=255)
                                                                             gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  1:1745460/100‑1 (MQ=255)
                                                                             gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  1:165572/100‑1 (MQ=255)
                                                                               acacGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  2:3640534/98‑1 (MQ=255)
                                                                               acacGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                        <  2:3590531/98‑1 (MQ=255)
                                                                                ctcGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGGCAGGCCTGGTACCACTTCTGCTGACCTTTGCTTGTATGTGGCAATTGCGCAAAAAAGTTa                     >  2:3421573/3‑100 (MQ=255)
                                                                                  cGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                   >  2:2369102/1‑100 (MQ=255)
                                                                                  cGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                   <  1:63861/100‑1 (MQ=255)
                                                                                   gTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                  <  1:1347730/100‑1 (MQ=255)
                                                                                    cattaactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                  >  2:1967680/10‑100 (MQ=38)
                                                                                      ttaactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCg                >  1:3830737/8‑100 (MQ=255)
                                                                                      actactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc               <  1:2344441/100‑1 (MQ=255)
                                                                                       ctactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt              >  1:4077369/1‑100 (MQ=255)
                                                                                         actggTCCAGGATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtg             >  2:1462807/5‑100 (MQ=255)
                                                                                          ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt            >  2:3187250/4‑100 (MQ=255)
                                                                                          ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt            >  1:2906601/4‑100 (MQ=255)
                                                                                           tggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg           >  1:1693133/3‑100 (MQ=255)
                                                                                           tGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAACTTTGATTGTATGTGGCTACTGCGCAAGAAAGTTAACCCGATGTgg          >  2:3464708/1‑100 (MQ=255)
                                                                                            ggTCCAGAATAGTTTGGACCCGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAGTGCGCAAAAAAGTTAACCCGCtgtg           >  2:3899357/2‑99 (MQ=255)
                                                                                            ggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTgg          >  1:3616907/2‑100 (MQ=255)
                                                                                             gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa         >  2:3952457/1‑100 (MQ=255)
                                                                                             gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa         >  1:688644/1‑100 (MQ=255)
                                                                                             gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa         >  2:3819148/1‑100 (MQ=255)
                                                                                               ccAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc       <  1:3330129/100‑1 (MQ=255)
                                                                                                cAGACTATGATGGACCAGTTAATGCCAGGCCTCATAAAACTGCTGCTGACCTTTGCTTGTATGTGAAAAAAAAAAAAAAAAGTTAACCCGCTGTGGatc       <  2:1936805/99‑1 (MQ=255)
                                                                                                   aCTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCg   <  2:499256/100‑1 (MQ=255)
                                                                                                   aCTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCg   <  1:89816/100‑1 (MQ=255)
                                                                                                   aCTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCg   <  2:1083618/100‑1 (MQ=255)
                                                                                                    cTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGt  >  2:3037544/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
TGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGT  >  NC_000913/1904493‑1904691

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