Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 888,315 C→T G243D (GGC→GAC)  b0847 ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913888,3150CT99.3% 498.9 / ‑5.4 141G243D (GGC→GAC) b0847putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (65/75);  total (65/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 153 total aligned reads displayed.

CAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888227‑888404
                                                                                        |                                                                                         
cAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGt                                                                                       >  1:2975814/1‑93 (MQ=255)
 aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCCCCGTATGGATTTGCCAGACTACCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      >  2:427171/1‑93 (MQ=255)
 aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      >  2:4227908/1‑93 (MQ=255)
 aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      >  1:3441431/1‑93 (MQ=255)
  aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGgtt                                                                                     >  2:188369/1‑93 (MQ=255)
   cGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTg                                                                                    >  1:3378498/1‑93 (MQ=255)
    gCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGt                                                                                   <  1:2685186/93‑1 (MQ=255)
    gCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGt                                                                                   <  1:2331860/93‑1 (MQ=255)
    gCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGt                                                                                   >  2:625496/1‑93 (MQ=255)
      tATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCg                                                                                 >  1:1796618/1‑93 (MQ=255)
         ttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGAtc                                                                              <  1:4423519/93‑2 (MQ=255)
          tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATaa                                                                             >  2:2963511/1‑93 (MQ=255)
          tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATaa                                                                             >  2:76899/1‑93 (MQ=255)
            aTCGCCCATTTGTAGCCCGGCATCACCGTCTGCATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTAGGTTTGTCGATAAAt                                                                           >  2:3389234/1‑93 (MQ=255)
            aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATCCGTTCAATGTAGCAGTCGGTTTGTCGATAAAt                                                                           <  1:2534292/93‑1 (MQ=255)
            aTCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAAt                                                                           <  2:3740897/93‑1 (MQ=255)
              ctcccATTTGTAGCACGGCATCCCCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAc                                                                         <  1:1134110/91‑1 (MQ=255)
              cGCCCATTTGTAGCGCGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAc                                                                         <  2:936320/93‑1 (MQ=255)
                cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       >  1:4650910/1‑93 (MQ=255)
                 ccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCa                                                                      <  1:1149403/93‑1 (MQ=255)
                  cATTTGTAGCACGTGATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  1:1742427/93‑1 (MQ=255)
                  cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  1:4650514/93‑1 (MQ=255)
                  cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     <  1:3684147/93‑1 (MQ=255)
                     ttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTc                                                                  >  1:1101667/1‑93 (MQ=255)
                      tGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCa                                                                 <  1:2895390/93‑1 (MQ=255)
                       gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                >  2:130297/1‑93 (MQ=255)
                       gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                >  1:2409351/1‑93 (MQ=255)
                        tAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACg                                                               <  1:1324934/93‑1 (MQ=255)
                         aGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGc                                                              >  2:4017343/1‑93 (MQ=255)
                          gCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCa                                                             >  1:4043140/1‑93 (MQ=255)
                            aCGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGa                                                           <  1:2234387/93‑1 (MQ=255)
                            aCGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGa                                                           >  2:3121019/1‑93 (MQ=255)
                            aCGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGa                                                           <  2:3780137/93‑1 (MQ=255)
                             cGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAt                                                          <  2:1093166/93‑1 (MQ=255)
                             cGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAt                                                          <  1:447323/93‑1 (MQ=255)
                             cGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAt                                                          <  1:2125194/93‑1 (MQ=255)
                               gCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttt                                                        <  2:2416803/93‑1 (MQ=255)
                                cATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAtttt                                                       <  2:197693/93‑1 (MQ=255)
                                cATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAtttt                                                       >  2:1603956/1‑93 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      <  1:3946063/93‑1 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  1:2230701/1‑93 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      <  1:3316508/93‑1 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  2:1871690/1‑93 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  1:2036695/1‑93 (MQ=255)
                                 aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTCCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  2:3362031/1‑93 (MQ=255)
                                  tCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTg                                                     >  2:441730/1‑93 (MQ=255)
                                    aCCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGcc                                                   <  2:3156987/93‑1 (MQ=255)
                                    aCCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGcc                                                   <  2:3883691/93‑1 (MQ=255)
                                      cGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGt                                                 <  1:164719/93‑1 (MQ=255)
                                      cGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGt                                                 >  1:1992945/1‑93 (MQ=255)
                                       gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTc                                                >  1:3974618/1‑93 (MQ=255)
                                       gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTc                                                >  1:2502539/1‑93 (MQ=255)
                                        tCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGTTTTGTCGATTAATACCCAGTTCACGCAGCTTTTTGCCGTCg                                               >  2:2810081/1‑93 (MQ=255)
                                        tCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCg                                               <  2:1000459/93‑1 (MQ=255)
                                        tCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCg                                               >  2:76982/1‑93 (MQ=255)
                                        tCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCg                                               >  2:3004258/1‑93 (MQ=255)
                                         cTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCgg                                              <  2:2748224/93‑1 (MQ=255)
                                          tGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGt                                             >  2:4107850/1‑93 (MQ=255)
                                            gATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTcc                                           <  2:388328/93‑1 (MQ=255)
                                            gATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTcc                                           >  2:2853851/1‑93 (MQ=255)
                                             aTTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCa                                          <  2:1945031/93‑1 (MQ=255)
                                              tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         <  1:2100806/93‑1 (MQ=255)
                                              tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         >  2:1695230/1‑93 (MQ=255)
                                              tttGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAg                                         >  2:633768/1‑93 (MQ=255)
                                                   cAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccac                                    <  1:3610688/93‑1 (MQ=255)
                                                      aaTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGATCACGCAGATTTTTGCAGTCGGTCCAGGCCACCAg                                 >  1:459504/1‑93 (MQ=255)
                                                        tCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCATGCAGATTTTTGCCGTCGGTCCAGGCCACCAGtt                               <  1:4569869/93‑1 (MQ=255)
                                                        tCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGtt                               <  2:3850448/93‑1 (MQ=255)
                                                        tCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGtt                               <  2:3653338/93‑1 (MQ=255)
                                                        tCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGtt                               >  2:2692609/1‑93 (MQ=255)
                                                        tCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGtt                               <  2:508186/93‑1 (MQ=255)
                                                          ccGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGCTTTTTGCCGTCGGTCCAGGCCACCAGTTcc                             >  1:1839621/1‑93 (MQ=255)
                                                          ccGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTcc                             >  2:4467040/1‑93 (MQ=255)
                                                           cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCATCGGTCCAGGCCACCAGTTCCg                            >  1:1052254/1‑93 (MQ=255)
                                                            gTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCgg                           <  1:3797482/93‑1 (MQ=255)
                                                             ttACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCggg                          <  2:2054864/93‑1 (MQ=255)
                                                             ttACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCggg                          >  2:1326646/1‑93 (MQ=255)
                                                              tACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGc                         <  1:4352769/93‑1 (MQ=255)
                                                                 gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      <  1:948136/93‑1 (MQ=255)
                                                                 gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      <  2:1748747/93‑1 (MQ=255)
                                                                  tCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAc                     <  2:1649984/93‑1 (MQ=255)
                                                                  tCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAc                     <  1:2821622/93‑1 (MQ=255)
                                                                   cGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAcg                    <  1:4527219/93‑1 (MQ=255)
                                                                    gAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAcgc                   >  1:1383467/1‑93 (MQ=255)
                                                                     aaTACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAcgcg                  >  1:1661205/1‑93 (MQ=255)
                                                                       tACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                <  1:493048/93‑1 (MQ=255)
                                                                       tACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                >  2:4418623/1‑93 (MQ=255)
                                                                         cGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAg              >  1:652700/1‑93 (MQ=255)
                                                                          gTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCCCCAGTTCCGGGCCGACGCCATAgg             >  2:1122026/1‑93 (MQ=255)
                                                                              aaTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         <  1:1331480/93‑1 (MQ=255)
                                                                                tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       >  2:3887210/1‑93 (MQ=255)
                                                                                tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       <  1:595920/93‑1 (MQ=255)
                                                                                tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       >  1:3821547/1‑93 (MQ=255)
                                                                                tGTAGCAGTCGGTTTGTCGATAAAGACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGa       >  2:4518102/1‑93 (MQ=255)
                                                                                  tAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATc     <  1:3725898/93‑1 (MQ=255)
                                                                                  tAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATc     <  2:2904982/93‑1 (MQ=255)
                                                                                   agcagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCa    <  2:4656423/93‑1 (MQ=255)
                                                                                   agcagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCa    >  1:4115647/1‑93 (MQ=255)
                                                                                    gcagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCAc   <  1:1428880/93‑1 (MQ=255)
                                                                                    gcagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCAc   <  1:3276326/93‑1 (MQ=255)
                                                                                        |                                                                                         
CAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888227‑888404

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