Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,904,595 T→G 21.2% I234S (ATT→AGT)  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,5950TG21.2% 257.7 / 76.9 127I234S (ATT→AGT) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (46/54);  new base G (15/12);  total (61/66)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 129 total aligned reads displayed.

GGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGG  >  NC_000913/1904497‑1904694
                                                                                                   |                                                                                                   
ggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTAtt                                                                                                    >  1:309801/1‑100 (MQ=255)
    ttGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTgg                                                                                                <  1:806252/100‑1 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGAc                                                                                              >  2:1451931/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGAc                                                                                              >  2:1109673/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGAc                                                                                              >  2:1929008/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGAc                                                                                              <  2:265989/100‑1 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACAGCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGa                                                                                               <  1:3922781/99‑1 (MQ=255)
       gTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAca                                                                                             <  2:1305958/100‑2 (MQ=255)
        ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCa                                                                                            >  1:1287180/1‑100 (MQ=255)
        ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcag                                                                                            <  1:3564726/100‑3 (MQ=255)
         tAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAg                                                                                           <  2:628408/100‑1 (MQ=255)
          aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagta                                                                                          <  1:3030032/100‑5 (MQ=255)
          aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagta                                                                                          <  2:2361136/100‑5 (MQ=255)
          aaCAAGTGGACACACGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagt                                                                                           <  2:3900311/100‑4 (MQ=255)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagta                                                                                          <  1:733076/100‑5 (MQ=38)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagtag                                                                                         <  2:837701/100‑6 (MQ=255)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagtag                                                                                         <  1:1805640/100‑6 (MQ=255)
             aaGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagtagta                                                                                       <  2:1082408/100‑8 (MQ=255)
              aGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAAt                                                                                      <  2:768860/100‑1 (MQ=255)
               gTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGAATAGTCTGGAcagtagtaac                                                                                     <  1:2466319/100‑10 (MQ=38)
                 ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagtagtaacg                                                                                    <  1:1522143/100‑11 (MQ=37)
                 ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGcc                                                                                   >  2:1501381/1‑100 (MQ=255)
                     acaTGTCGACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGcc                                                                                   <  2:2711133/96‑1 (MQ=255)
                       aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCTAGGc                                                                               >  2:3482157/1‑98 (MQ=255)
                       aTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCt                                                                             <  1:1945542/100‑1 (MQ=255)
                         gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTgg                                                                           <  2:3223724/100‑1 (MQ=255)
                           cAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                         <  2:2800965/100‑1 (MQ=255)
                             aCATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTAcc                                                                       >  2:1732009/1‑100 (MQ=255)
                               aTCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCAc                                                                     <  2:2148118/100‑1 (MQ=255)
                                 cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATG‑CAGGCCTGGTACCActgc                                                                  <  2:2072504/100‑1 (MQ=255)
                                  ccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgc                                                                  >  2:1252411/1‑100 (MQ=255)
                                    gCTGGTTGTCTCTCGCATTACTGACTATACGGGCAAATAAAATGTTACTACT‑GTCCAGACTATTCTGGACTAGTTAATGCAAGGCCTGGTACCAct                                                                    <  2:1866872/96‑1 (MQ=255)
                                    gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                                >  1:3721150/1‑100 (MQ=255)
                                    gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                                <  1:3871167/100‑1 (MQ=255)
                                    gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                                <  1:639045/100‑1 (MQ=255)
                                     cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                               <  2:2883703/100‑1 (MQ=255)
                                     cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                               <  1:1829235/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:1731879/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:1731867/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:2749846/100‑1 (MQ=255)
                                           gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCt                                                         >  2:1952295/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  1:1208926/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  2:2248797/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  2:1452779/1‑100 (MQ=255)
                                               ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                     >  1:1813845/1‑100 (MQ=255)
                                               ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                     >  2:568608/1‑100 (MQ=255)
                                                tcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                    <  1:2848294/100‑1 (MQ=255)
                                                 cGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                    <  1:1318200/99‑1 (MQ=255)
                                                  gCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAAGGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                  <  1:2957941/100‑1 (MQ=255)
                                                        cTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgtg                                            >  1:1597789/1‑100 (MQ=255)
                                                          gACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                          <  2:2447263/100‑1 (MQ=255)
                                                          gACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                          >  1:2002693/1‑100 (MQ=255)
                                                          gACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATGCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                           >  2:2678037/1‑99 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         <  1:1664643/100‑1 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         >  2:759045/1‑100 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         >  1:2212172/1‑100 (MQ=255)
                                                             cAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAc                                       >  1:2089430/1‑100 (MQ=255)
                                                               gACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCATGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                     <  2:1562162/100‑1 (MQ=255)
                                                               gACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                     >  2:165572/1‑100 (MQ=255)
                                                                 cGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                   >  2:905718/1‑100 (MQ=255)
                                                                 cGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                   >  1:2519404/1‑100 (MQ=255)
                                                                   ggCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCa                                 >  1:3259251/1‑100 (MQ=255)
                                                                    gCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                                >  1:2404348/1‑100 (MQ=255)
                                                                       aaGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaa                             >  1:1646461/1‑100 (MQ=255)
                                                                         gAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:1745460/100‑1 (MQ=255)
                                                                         gAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:1746442/100‑1 (MQ=255)
                                                                         gAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:2605649/100‑1 (MQ=255)
                                                                           acacGTTACTACT‑GTCCAGACTATTCTGGACGAGTTAATGCTAGGCTTGGTACCATTGCTGCTGACCTTTGCTTGGATCTGGCTACTGCGCAAAAAAg                           >  2:1761528/1‑98 (MQ=255)
                                                                           acacGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  2:3640534/98‑1 (MQ=255)
                                                                           acacGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  2:3590531/98‑1 (MQ=255)
                                                                            ctcGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGGCAGGCCTGGTACCACTTCTGCTGACCTTTGCTTGTATGTGGCAATTGCGCAAAAAAGTTa                        >  2:3421573/3‑100 (MQ=255)
                                                                              cGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                      <  1:63861/100‑1 (MQ=255)
                                                                              cGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                      >  2:2369102/1‑100 (MQ=255)
                                                                               gTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                     <  1:1347730/100‑1 (MQ=255)
                                                                                  ttaactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCg                   >  1:3830737/8‑100 (MQ=255)
                                                                                  actact‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc                  <  1:2344441/100‑1 (MQ=255)
                                                                                   ctact‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                 >  1:4077369/1‑100 (MQ=255)
                                                                                    aactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                 >  2:2697097/6‑100 (MQ=255)
                                                                                     actggTCCAGAATAGTCTGGAGCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTGGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtg                >  2:2938851/5‑100 (MQ=255)
                                                                                      ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt               >  1:2906977/4‑100 (MQ=255)
                                                                                      ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt               >  2:823215/4‑100 (MQ=255)
                                                                                      ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt               >  2:3187250/4‑100 (MQ=255)
                                                                                      ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt               >  1:2361804/4‑100 (MQ=255)
                                                                                       tggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg              >  1:1693133/3‑100 (MQ=255)
                                                                                       t‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAACTTTGATTGTATGTGGCTACTGCGCAAGAAAGTTAACCCGATGTgg             >  2:3464708/1‑100 (MQ=255)
                                                                                        ggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTgg             >  1:3616907/2‑100 (MQ=255)
                                                                                         gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa            >  2:3952457/1‑100 (MQ=255)
                                                                                         gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa            >  2:3819148/1‑100 (MQ=255)
                                                                                           ccAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc          <  1:3247940/100‑1 (MQ=255)
                                                                                           ccAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc          >  1:592151/1‑100 (MQ=255)
                                                                                           ccAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatc          <  1:3330129/100‑1 (MQ=255)
                                                                                            cAGACTATGATGGACCAGTTAATGCCAGGCCTCATAAAACTGCTGCTGACCTTTGCTTGTATGTGAAAAAAAAAAAAAAAAGTTAACCCGCTGTGGatc          <  2:1936805/99‑1 (MQ=255)
                                                                                             aGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatcat        >  1:76375/1‑100 (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)
                                                                                                  aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg   <  2:3160720/100‑1 (MQ=255)
                                                                                                  aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCCTTGCTTGTATGTGAAAAATAAACAAAAAAGTTAACCCGCTGTGGATCATCGTTg   <  2:1530618/100‑1 (MQ=255)
                                                                                                   ttCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTgg  <  2:3438990/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACT‑GTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGG  >  NC_000913/1904497‑1904694

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