Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,244,543 T→C F61L (TTC→CTC)  malK → ATP‑binding component of transport system for maltose

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,244,5430TC96.0% 448.0 / 5.7 149F61L (TTC→CTC) malKATP‑binding component of transport system for maltose
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base C (61/82);  total (63/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 151 total aligned reads displayed.

CGTGGTGTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCACCT  >  NC_000913/4244455‑4244631
                                                                                        |                                                                                        
 gtggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     <  1:429125/93‑1 (MQ=255)
 gtggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     <  1:429124/93‑1 (MQ=255)
 gtggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     >  2:215742/1‑93 (MQ=255)
 gtggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     <  2:4300587/93‑1 (MQ=255)
 gtggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     >  2:3073394/1‑93 (MQ=255)
 gtggtgTTGGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATc                                                                                     <  1:1074154/93‑1 (MQ=255)
   ggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCgg                                                                                   >  2:2112692/1‑93 (MQ=255)
   ggtgTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCgg                                                                                   <  1:3194615/93‑1 (MQ=255)
    gtgtTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGt                                                                                  <  2:2525021/93‑1 (MQ=255)
     tgtTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTg                                                                                 >  1:2194151/1‑93 (MQ=255)
      gtTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATTATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCCGTgg                                                                                >  2:3813792/1‑92 (MQ=39)
         tGTCGGCCCGTCTGGCTGCGGTAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGaaa                                                                            <  1:2597108/93‑1 (MQ=37)
           tCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAAc                                                                           <  2:3793522/93‑1 (MQ=255)
            cGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACg                                                                          <  1:606912/93‑1 (MQ=255)
              gACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGa                                                                        <  2:3417037/93‑1 (MQ=255)
              gACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGa                                                                        <  2:4510026/93‑1 (MQ=255)
              gACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGa                                                                        <  2:3417036/93‑1 (MQ=255)
              gACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGa                                                                        <  1:1597265/93‑1 (MQ=255)
              gACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGa                                                                        <  2:2174063/93‑1 (MQ=255)
               aCCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGat                                                                       <  2:2035482/93‑1 (MQ=255)
                 cGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatga                                                                     >  1:1423523/1‑93 (MQ=255)
                  gTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGGCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaa                                                                    <  2:2695871/93‑1 (MQ=255)
                  gTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaa                                                                    <  1:3159170/93‑1 (MQ=255)
                  gTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaa                                                                    >  2:691176/1‑93 (MQ=255)
                   tCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaat                                                                   <  1:1206128/93‑1 (MQ=255)
                     tGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaatga                                                                 >  1:210673/1‑93 (MQ=255)
                     tGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGatgaatga                                                                 >  2:4300300/1‑93 (MQ=255)
                       gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGaca                                                               <  1:3620303/93‑1 (MQ=255)
                       gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGaca                                                               >  2:2495911/1‑93 (MQ=255)
                        cTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGacac                                                              <  2:4138551/93‑1 (MQ=255)
                          gCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTc                                                            >  2:290672/1‑93 (MQ=255)
                           cGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTcc                                                           >  2:569439/1‑93 (MQ=255)
                           cGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTcc                                                           <  2:122820/93‑1 (MQ=255)
                           cGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTcc                                                           <  1:4421228/93‑1 (MQ=255)
                             gTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGc                                                         >  2:4011691/1‑93 (MQ=255)
                               aaaTCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGACACTCCGCca                                                       >  1:2007267/1‑93 (MQ=255)
                               aaaTCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGccg                                                       >  1:1564957/1‑92 (MQ=255)
                                aaTCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCcag                                                      >  2:2655985/1‑93 (MQ=255)
                                aaTCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGAGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCCCcac                                                      >  2:3997844/1‑92 (MQ=38)
                                 aTCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCcagc                                                     <  1:597634/93‑1 (MQ=255)
                                   cGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCcagcag                                                   >  2:2816606/1‑93 (MQ=255)
                                   cGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCcagcag                                                   >  1:4428802/1‑93 (MQ=255)
                                    gACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGa                                                  >  1:1933976/1‑93 (MQ=255)
                                     aCTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGaa                                                 <  1:379217/93‑1 (MQ=255)
                                      cTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTTATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAc                                                <  2:1555187/93‑1 (MQ=255)
                                      cTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAc                                                <  2:4419912/93‑1 (MQ=255)
                                      cTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAc                                                <  1:2403678/93‑1 (MQ=255)
                                       tttACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAcg                                               >  1:2936933/1‑93 (MQ=255)
                                        ttACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAcgc                                              <  2:3827374/93‑1 (MQ=255)
                                         tACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAAcgcg                                             >  1:3965667/1‑93 (MQ=255)
                                          aCTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcgg                                            >  1:2456448/1‑93 (MQ=255)
                                          aCTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcgg                                            <  1:3640637/93‑1 (MQ=255)
                                           cTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcggc                                           >  1:2691392/1‑93 (MQ=255)
                                           cTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcggc                                           <  2:4446121/93‑1 (MQ=255)
                                            tGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcggcg                                          >  1:684979/1‑93 (MQ=255)
                                            tGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACgcggcg                                          >  2:3169929/1‑93 (MQ=255)
                                             gcgcATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGt                                         <  2:3110335/93‑1 (MQ=255)
                                             gcgcATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGt                                         <  2:4160979/93‑1 (MQ=255)
                                               gcATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGCCGTTg                                       >  1:1653494/1‑93 (MQ=255)
                                                cATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTgg                                      >  2:242185/1‑93 (MQ=255)
                                                 aTGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGACCCTCCGCCAGCAGAACGCGGCGTTGGt                                     <  2:3756092/93‑1 (MQ=255)
                                                 aTGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGt                                     >  1:776532/1‑93 (MQ=255)
                                                   gATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTAt                                   >  2:3624819/1‑93 (MQ=255)
                                                   gATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTAt                                   <  2:2732421/93‑1 (MQ=255)
                                                     ttGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATgg                                 <  1:2173351/93‑1 (MQ=255)
                                                      tGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgt                                <  2:2712162/93‑1 (MQ=255)
                                                      tGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgt                                <  1:4231306/93‑1 (MQ=255)
                                                      tGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgt                                <  2:1238362/93‑1 (MQ=255)
                                                      tGCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgt                                <  1:900596/93‑1 (MQ=255)
                                                       gCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGCACGCGGCGTTGGTATGgtg                               >  1:1745377/1‑93 (MQ=255)
                                                       gCCGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgtg                               <  1:665837/93‑1 (MQ=255)
                                                        gcGGGCGTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGgtgt                              <  1:3491702/92‑1 (MQ=255)
                                                         ccggCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGtt                             <  2:3238599/91‑1 (MQ=255)
                                                         cGGGCTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGtt                             >  1:3332020/1‑93 (MQ=255)
                                                             cTTGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAg                         <  1:3815976/93‑1 (MQ=255)
                                                              ttGAGACGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGt                        >  2:3842721/1‑93 (MQ=255)
                                                                gagaCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCt                      <  1:2389523/93‑1 (MQ=255)
                                                                gagaCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCt                      >  1:1865882/1‑93 (MQ=255)
                                                                 agaCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCtt                     >  2:66117/1‑93 (MQ=255)
                                                                  gaCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTa                    >  2:507586/1‑93 (MQ=255)
                                                                  gaCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTa                    <  2:2036775/93‑1 (MQ=255)
                                                                   aCGATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTAc                   <  2:2523471/93‑1 (MQ=255)
                                                                     gATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTAcgc                 <  2:2094967/93‑1 (MQ=255)
                                                                     gATCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTAcgc                 >  1:3830341/1‑93 (MQ=255)
                                                                      aTCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTAcgcg                >  2:2895427/1‑93 (MQ=255)
                                                                       tCACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACgcgc               <  1:2761755/93‑1 (MQ=255)
                                                                        cACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATTAATGCCACTCCGCCAGCAGGACCCGGCGTTGGTATGGGGTTTCATTCTTACGCGct              >  1:3684909/1‑93 (MQ=255)
                                                                        cACCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGct              >  1:3684910/1‑93 (MQ=255)
                                                                         aCCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGctc             <  2:4512460/93‑1 (MQ=255)
                                                                         aCCAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGctc             >  1:2874723/1‑93 (MQ=255)
                                                                           cAGCGGCGACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGAATGGTGTTTCAGTCTTACGCGCTCTa           >  1:1971347/1‑93 (MQ=255)
                                                                              cggcgACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATcc        <  1:1118857/93‑1 (MQ=255)
                                                                               ggcgACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATccc       <  2:1546961/93‑1 (MQ=255)
                                                                                gcgACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATcccc      <  1:1795284/93‑1 (MQ=255)
                                                                                gcgACCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATcccc      <  2:929514/93‑1 (MQ=255)
                                                                                   aCCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCAcc   <  1:829438/93‑1 (MQ=255)
                                                                                   aCCTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCAcc   <  1:3329209/93‑1 (MQ=255)
                                                                                    ccTGTTCATCGGTGAGAAACGGATGAATGACACTCCGCCCGCAGAACGCGGCGGTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCACCt  >  2:2810196/1‑93 (MQ=255)
                                                                                    ccTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAGCGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCACCt  <  1:4232631/93‑1 (MQ=255)
                                                                                    ccTGCTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCACCt  >  1:3494621/1‑93 (MQ=255)
                                                                                        |                                                                                        
CGTGGTGTTTGTCGGACCGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGACACTCCGCCAGCAGAACGCGGCGTTGGTATGGTGTTTCAGTCTTACGCGCTCTATCCCCACCT  >  NC_000913/4244455‑4244631

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