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,3150CT97.0% 447.7 / ‑0.2 134G243D (GGC→GAC) b0847putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (65/64);  total (67/67)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 141 total aligned reads displayed.

CCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888226‑888404
                                                                                         |                                                                                         
ccAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCgg                                                                                        <  2:3897759/93‑1 (MQ=255)
ccAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCgg                                                                                        >  1:1568083/1‑93 (MQ=255)
 cAACGCTATTTCATCGCCCATTTGTAGCCCGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGt                                                                                       <  1:1075635/93‑1 (MQ=255)
 cAACGCTATTTCATCGCCCATTTGTAGCACGGCATCCCCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGt                                                                                       <  1:1075634/93‑1 (MQ=255)
 cAACGCTATTTCATCGCCCATCTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGt                                                                                       >  1:1440037/1‑93 (MQ=255)
  aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      >  2:929079/1‑93 (MQ=255)
  aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      <  2:2573978/93‑1 (MQ=255)
  aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      >  1:489068/1‑93 (MQ=255)
  aaCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGtt                                                                                      <  1:3373358/93‑1 (MQ=255)
   aCGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGttt                                                                                     <  2:1757002/93‑1 (MQ=255)
     gCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGt                                                                                   >  2:2810207/1‑93 (MQ=255)
      cTATTTCATCGCCCATTTGTAGCACGGCCTCACCGTCTGGATTTGCCGGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTc                                                                                  <  2:4581075/93‑1 (MQ=255)
      cTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTc                                                                                  <  2:97398/93‑1 (MQ=255)
       tATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCg                                                                                 >  2:529486/1‑93 (MQ=255)
       tATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCg                                                                                 <  1:4597910/93‑1 (MQ=255)
       tATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCg                                                                                 >  1:3070624/1‑93 (MQ=255)
        ttttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGa                                                                                <  1:3787892/92‑1 (MQ=255)
        aTTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGa                                                                                >  2:2829945/1‑93 (MQ=255)
         tttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGAt                                                                               <  2:857472/93‑1 (MQ=255)
         tttCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGAt                                                                               >  1:4190297/1‑93 (MQ=255)
           tCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATaa                                                                             <  1:212177/93‑1 (MQ=255)
                gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAcc                                                                        >  1:1725730/1‑93 (MQ=255)
                gCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAcc                                                                        >  2:3041642/1‑93 (MQ=255)
                 cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       >  1:3054540/1‑93 (MQ=255)
                 cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       >  2:2636325/1‑93 (MQ=255)
                 cccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATAccc                                                                       <  2:2055040/93‑1 (MQ=255)
                  ccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCa                                                                      <  2:2080006/93‑1 (MQ=255)
                  ccATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCa                                                                      <  2:2079894/93‑1 (MQ=255)
                   cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     >  1:1820810/1‑93 (MQ=255)
                   cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     >  2:881971/1‑93 (MQ=255)
                   cATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAg                                                                     >  2:3627303/1‑93 (MQ=255)
                     tttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGtt                                                                   >  2:1338167/1‑93 (MQ=255)
                     tttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGtt                                                                   <  2:1606946/93‑1 (MQ=255)
                     tttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGtt                                                                   >  2:2371409/1‑93 (MQ=255)
                     tttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGtt                                                                   >  2:1790357/1‑93 (MQ=255)
                     tttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGtt                                                                   <  1:2641668/93‑1 (MQ=255)
                      ttGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTc                                                                  <  1:3046620/93‑1 (MQ=255)
                       tGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCa                                                                 <  1:79107/93‑1 (MQ=255)
                        gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                >  2:461877/1‑93 (MQ=255)
                        gTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCAc                                                                <  2:3337461/93‑1 (MQ=255)
                           gCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCa                                                             <  2:3847277/93‑1 (MQ=255)
                           gCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTAACGCa                                                             >  2:4528188/1‑93 (MQ=255)
                              cGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAt                                                          <  2:3044535/93‑1 (MQ=255)
                              cGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAt                                                          >  2:269065/1‑93 (MQ=255)
                               ggCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAtt                                                         <  2:2022642/93‑1 (MQ=255)
                               ggCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAtt                                                         <  2:1409415/93‑1 (MQ=255)
                                gCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttt                                                        <  2:3641095/93‑1 (MQ=255)
                                gCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttt                                                        >  1:1001587/1‑93 (MQ=255)
                                 cATCACCGTCGGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAtttt                                                       <  2:3880595/93‑1 (MQ=255)
                                  aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  1:2774165/1‑93 (MQ=255)
                                  aTCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGAttttt                                                      >  2:4274534/1‑93 (MQ=255)
                                   tCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTg                                                     >  1:3111239/1‑93 (MQ=255)
                                    cACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGc                                                    <  2:3126944/93‑1 (MQ=255)
                                    cACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGc                                                    <  2:1029802/93‑1 (MQ=255)
                                    cACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGc                                                    <  2:3973149/93‑1 (MQ=255)
                                       cGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGt                                                 >  2:3566518/1‑93 (MQ=255)
                                        gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTc                                                >  2:1139089/1‑93 (MQ=255)
                                        gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTc                                                >  2:2868982/1‑93 (MQ=255)
                                        gTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACACAGTTCACGCAGATTTTTGCCGTc                                                <  2:3959137/93‑1 (MQ=255)
                                                  gCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGcc                                      <  1:2482374/93‑1 (MQ=255)
                                                    cAGAATCCCGTTACGTCGGATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccac                                    >  1:2736787/1‑93 (MQ=255)
                                                    cAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccac                                    >  2:396771/1‑93 (MQ=255)
                                                     aGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacc                                   <  1:3795556/93‑1 (MQ=255)
                                                      gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCTGTCCAGGccacca                                  >  1:3155133/1‑93 (MQ=255)
                                                      gAATCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGccacca                                  <  2:3019407/93‑1 (MQ=255)
                                                       aaTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAg                                 <  2:4490631/93‑1 (MQ=255)
                                                       aaTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAg                                 >  2:3208206/1‑93 (MQ=255)
                                                       aaTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAg                                 >  1:3325046/1‑93 (MQ=255)
                                                        aTCCCGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGt                                >  2:113389/1‑93 (MQ=255)
                                                           ccGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTcc                             <  2:805341/93‑1 (MQ=255)
                                                           ccGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTcc                             <  2:805342/93‑1 (MQ=255)
                                                            cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCGGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCg                            <  2:3822123/93‑1 (MQ=255)
                                                            cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCg                            <  1:2345443/93‑1 (MQ=255)
                                                            cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCg                            <  2:1596564/93‑1 (MQ=255)
                                                            cGTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCg                            <  1:844156/93‑1 (MQ=255)
                                                             gTTACGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCgg                           >  2:1262890/1‑93 (MQ=255)
                                                                 cGTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCg                       >  1:3013666/1‑93 (MQ=255)
                                                                  gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      <  1:1311713/93‑1 (MQ=255)
                                                                  gTCGAATACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGa                      >  1:4311481/1‑93 (MQ=255)
                                                                      aaTACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGAcgcg                  >  1:490774/1‑93 (MQ=255)
                                                                       aTACGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGa                 >  2:4553624/1‑93 (MQ=255)
                                                                        tACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGAt                <  1:3155573/93‑1 (MQ=255)
                                                                         aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  1:124287/93‑1 (MQ=255)
                                                                         aCGTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCAGTCCAGGCCACCAGTTCCGGGCCGACGCGATa               <  1:4194368/93‑1 (MQ=255)
                                                                           gTTCAATGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAgg             <  2:426950/93‑1 (MQ=255)
                                                                            ttCAATGTAGCAGTCGGTTTGTCGATAAAGACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGgc            <  1:4601949/93‑1 (MQ=255)
                                                                               aaTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         <  2:2058399/93‑1 (MQ=255)
                                                                               aaTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         <  2:2333903/93‑1 (MQ=255)
                                                                               aaTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         >  1:3750261/1‑93 (MQ=255)
                                                                               aaTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGcgcg         >  1:1231364/1‑93 (MQ=255)
                                                                                aTGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCgg        <  2:2267624/93‑1 (MQ=255)
                                                                                 tGTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCggc       >  1:1334502/1‑92 (MQ=255)
                                                                                  gTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGGCGCGATAGGCGCGGAt      >  1:2571389/1‑93 (MQ=255)
                                                                                  gTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGAt      >  1:2571390/1‑93 (MQ=255)
                                                                                  gTAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGAt      <  1:3659117/93‑1 (MQ=255)
                                                                                   tAGCAGTCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATc     <  1:3567044/93‑1 (MQ=255)
                                                                                    agcagtCGGTTTGTCGATAAATAACCAGTTCACGCAGATTTTTGACGTCGGTCCAGGCTACCAGTTCCGGGCCGACGAGATAGGCGCGGATCa    >  2:1170796/1‑93 (MQ=255)
                                                                                     gcagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCAc   >  2:1428767/1‑93 (MQ=255)
                                                                                      cagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCAcc  >  1:3317261/1‑93 (MQ=255)
                                                                                      cagtCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCAcc  <  1:1311388/93‑1 (MQ=255)
                                                                                         |                                                                                         
CCAACGCTATTTCATCGCCCATTTGTAGCACGGCATCACCGTCTGGATTTGCCAGAATCCCGTTACGTCGAATACGTTCAATGTAGCAGCCGGTTTGTCGATAAATACCCAGTTCACGCAGATTTTTGCCGTCGGTCCAGGCCACCAGTTCCGGGCCGACGCGATAGGCGCGGATCACC  >  NC_000913/888226‑888404

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