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.

GGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGG  >  NC_000913/1904497‑1904694
                                                                                                   |                                                                                                   
ggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAtt                                                                                                    >  1:309801/1‑100 (MQ=255)
    ttGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTgg                                                                                                <  1:806252/100‑1 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                              >  2:1451931/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                              >  2:1109673/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                              >  2:1929008/1‑100 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                              <  2:265989/100‑1 (MQ=255)
      ggTTAACAAGTGGACACATGTCAACAGCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGa                                                                                               <  1:3922781/99‑1 (MQ=255)
       gTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAca                                                                                             <  2:1305958/100‑2 (MQ=255)
        ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCa                                                                                            >  1:1287180/1‑100 (MQ=255)
        ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                            <  1:3564726/100‑3 (MQ=255)
         tAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAg                                                                                           <  2:628408/100‑1 (MQ=255)
          aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagta                                                                                          <  1:3030032/100‑5 (MQ=255)
          aaCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagta                                                                                          <  2:2361136/100‑5 (MQ=255)
          aaCAAGTGGACACACGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagt                                                                                           <  2:3900311/100‑4 (MQ=255)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagta                                                                                          <  1:733076/100‑5 (MQ=38)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                         <  2:837701/100‑6 (MQ=255)
           aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                         <  1:1805640/100‑6 (MQ=255)
             aaGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtagta                                                                                       <  2:1082408/100‑8 (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)
                     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)
                                 cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCAGGCCTGGTACCActgc                                                                  <  2:2072504/100‑1 (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:3871167/100‑1 (MQ=255)
                                    gCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctg                                                                <  1:639045/100‑1 (MQ=255)
                                     cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                               <  2:2883703/100‑1 (MQ=255)
                                     cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                               <  1:1829235/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:1731879/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGATGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:1731867/100‑1 (MQ=255)
                                          tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                          <  1:2749846/100‑1 (MQ=255)
                                           gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCt                                                         >  2:1952295/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  1:1208926/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  2:2248797/1‑100 (MQ=255)
                                             ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                       >  2:1452779/1‑100 (MQ=255)
                                               ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                     >  1:1813845/1‑100 (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                                          <  2:2447263/100‑1 (MQ=255)
                                                          gACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                          >  1:2002693/1‑100 (MQ=255)
                                                          gACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATGCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                           >  2:2678037/1‑99 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         <  1:1664643/100‑1 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         >  2:759045/1‑100 (MQ=255)
                                                           aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCt                                         >  1:2212172/1‑100 (MQ=255)
                                                             cAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAc                                       >  1:2089430/1‑100 (MQ=255)
                                                               gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCATGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                     <  2:1562162/100‑1 (MQ=255)
                                                               gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                     >  2:165572/1‑100 (MQ=255)
                                                                 cGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                   >  2:905718/1‑100 (MQ=255)
                                                                 cGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                   >  1:2519404/1‑100 (MQ=255)
                                                                   ggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCa                                 >  1:3259251/1‑100 (MQ=255)
                                                                    gCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                                >  1:2404348/1‑100 (MQ=255)
                                                                       aaGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaa                             >  1:1646461/1‑100 (MQ=255)
                                                                         gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:1745460/100‑1 (MQ=255)
                                                                         gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:1746442/100‑1 (MQ=255)
                                                                         gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  1:2605649/100‑1 (MQ=255)
                                                                           acacGTTACTACTGTCCAGACTATTCTGGACGAGTTAATGCTAGGCTTGGTACCATTGCTGCTGACCTTTGCTTGGATCTGGCTACTGCGCAAAAAAg                           >  2:1761528/1‑98 (MQ=255)
                                                                           acacGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  2:3640534/98‑1 (MQ=255)
                                                                           acacGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                           <  2:3590531/98‑1 (MQ=255)
                                                                            ctcGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGGCAGGCCTGGTACCACTTCTGCTGACCTTTGCTTGTATGTGGCAATTGCGCAAAAAAGTTa                        >  2:3421573/3‑100 (MQ=255)
                                                                              cGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                      <  1:63861/100‑1 (MQ=255)
                                                                              cGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                      >  2:2369102/1‑100 (MQ=255)
                                                                               gTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                     <  1:1347730/100‑1 (MQ=255)
                                                                                  ttaactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCg                   >  1:3830737/8‑100 (MQ=255)
                                                                                  actactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc                  <  1:2344441/100‑1 (MQ=255)
                                                                                   ctactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                 >  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)
                                                                                       tGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAACTTTGATTGTATGTGGCTACTGCGCAAGAAAGTTAACCCGATGTgg             >  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)
                                                                                                   |                                                                                                   
GGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGG  >  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: