Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,220,610 A→G E61G (GAG→GGG)  b1170 → putative part of putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,220,6100AG98.5% 399.5 / ‑2.5 130E61G (GAG→GGG) b1170putative part of putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (66/61);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 136 total aligned reads displayed.

CTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGAGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGCT  >  NC_000913/1220520‑1220699
                                                                                          |                                                                                         
cTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGc                                                                                         >  2:719292/1‑93 (MQ=255)
cTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGGGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGc                                                                                         >  1:1173293/1‑93 (MQ=255)
 tAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCg                                                                                        <  2:4008266/93‑1 (MQ=255)
  aaCCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGt                                                                                       <  1:2805791/93‑1 (MQ=255)
  aaCCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCCACATCGATCCTGTAACTGAACAGGGGCGt                                                                                       <  2:529519/93‑1 (MQ=255)
    ccTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTc                                                                                     >  2:2110555/1‑93 (MQ=255)
     cTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCa                                                                                    >  2:4293613/1‑93 (MQ=255)
     cTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCa                                                                                    >  2:2870408/1‑93 (MQ=255)
      tGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCaa                                                                                   >  1:3391725/1‑93 (MQ=255)
      tGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCaa                                                                                   <  2:4446322/93‑1 (MQ=255)
      tGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCaa                                                                                   <  2:569166/93‑1 (MQ=255)
      tGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCaa                                                                                   >  1:196640/1‑93 (MQ=255)
       ggCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAg                                                                                  >  2:1458903/1‑93 (MQ=255)
        gCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGc                                                                                 <  1:4020523/93‑1 (MQ=255)
        gCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGc                                                                                 <  1:1106102/93‑1 (MQ=255)
         cAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGcc                                                                                <  1:4346670/93‑1 (MQ=255)
           gctgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGa                                                                              >  2:3738210/1‑93 (MQ=255)
           gctgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGa                                                                              <  2:1976277/93‑1 (MQ=255)
            ctgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGAc                                                                             <  1:4543346/93‑1 (MQ=255)
            ctgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGAc                                                                             >  1:3398427/1‑93 (MQ=255)
             tgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACt                                                                            >  1:3911374/1‑93 (MQ=255)
             tgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACt                                                                            >  2:4311361/1‑93 (MQ=255)
             tgctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTAAAGCCGACt                                                                            >  2:2179241/1‑93 (MQ=255)
              gctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACtt                                                                           <  2:4029362/93‑1 (MQ=255)
               ctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACttt                                                                          <  1:4264663/93‑1 (MQ=255)
               ctAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACttt                                                                          <  2:505991/93‑1 (MQ=255)
                tAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTg                                                                         <  2:599695/93‑1 (MQ=255)
                tAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTg                                                                         >  1:1347842/1‑93 (MQ=255)
                tAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTg                                                                         >  1:2537477/1‑93 (MQ=255)
                   cacTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCt                                                                      >  1:698590/1‑93 (MQ=255)
                   cacTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCt                                                                      >  1:2311169/1‑93 (MQ=255)
                    acTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTa                                                                     >  2:3339330/1‑93 (MQ=255)
                    acTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTa                                                                     >  2:803006/1‑93 (MQ=255)
                    acTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTa                                                                     >  2:3614963/1‑93 (MQ=255)
                      tctTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACg                                                                   >  1:3956600/1‑93 (MQ=255)
                       ctTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGt                                                                  <  1:108000/93‑1 (MQ=255)
                        tttttGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTc                                                                 >  2:2338467/1‑93 (MQ=255)
                         ttttGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCa                                                                >  2:616263/1‑93 (MQ=255)
                          tttGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCaa                                                               <  1:2100773/93‑1 (MQ=255)
                            tGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAAt                                                             <  2:3713185/93‑1 (MQ=255)
                             gTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAAtt                                                            >  1:4114434/1‑93 (MQ=255)
                             gTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAAtt                                                            <  1:3904859/93‑1 (MQ=255)
                             gTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGGTCAAGCCGACTTTGGCTACGTCAAAtt                                                            >  1:1860311/1‑93 (MQ=255)
                              ttATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCACGCCGACTTTGGCTACGTCAAATTg                                                           >  1:1550132/1‑93 (MQ=255)
                               tATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTgg                                                          <  1:3282776/93‑1 (MQ=255)
                                aTGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTggg                                                         >  1:3310755/1‑92 (MQ=255)
                                aTGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGc                                                         <  1:745985/93‑1 (MQ=255)
                                aTGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGc                                                         <  2:2674268/93‑1 (MQ=255)
                                  gCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCgg                                                       <  1:615342/93‑1 (MQ=255)
                                     tCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCa                                                    <  2:2613846/93‑1 (MQ=255)
                                       tGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATa                                                  >  2:3304637/1‑93 (MQ=255)
                                        gAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATaa                                                 >  1:3688234/1‑93 (MQ=255)
                                        gAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATaa                                                 <  1:137509/93‑1 (MQ=255)
                                        gAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATaa                                                 >  2:1615028/1‑93 (MQ=255)
                                           cGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGc                                              >  1:1094158/1‑93 (MQ=255)
                                            gACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATggc                                             >  2:1228034/1‑91 (MQ=255)
                                            gACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGcc                                             >  1:117764/1‑93 (MQ=255)
                                              ccGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTg                                           >  1:3411400/1‑93 (MQ=255)
                                               cGTGCGGGTGAACCCCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTgg                                          <  2:2656767/93‑1 (MQ=255)
                                               cGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTgg                                          >  2:660391/1‑93 (MQ=255)
                                                gTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGACTGGc                                         >  1:2231116/1‑93 (MQ=255)
                                                gTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGAGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGc                                         <  2:2363536/93‑1 (MQ=255)
                                                   cGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTg                                      <  1:926781/93‑1 (MQ=255)
                                                     ggTGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGAc                                    >  2:4426703/1‑93 (MQ=255)
                                                       tGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAg                                  >  2:899434/1‑93 (MQ=255)
                                                       tGAAACGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAg                                  <  2:211701/93‑1 (MQ=255)
                                                         aaaCGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCa                                <  1:1542/93‑1 (MQ=255)
                                                         aaaCGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCa                                >  1:234471/1‑93 (MQ=255)
                                                          aaCGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCaa                               >  2:412397/1‑93 (MQ=255)
                                                           aCGCGCTACATCGATCCTGTAACTGAACAGGTGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAAc                              <  1:1344667/93‑1 (MQ=255)
                                                           aCGCGCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAAc                              <  1:3098167/93‑1 (MQ=255)
                                                            cgcgCTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCCGGCGTGACAGCAACg                             >  2:664301/1‑93 (MQ=255)
                                                             gcgcTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACgg                            >  2:1828454/1‑93 (MQ=255)
                                                              cgcTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGa                           <  2:1298668/93‑1 (MQ=255)
                                                                cTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACa                         <  2:2864491/93‑1 (MQ=255)
                                                                cTACATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACa                         <  1:3667294/93‑1 (MQ=255)
                                                                 tACATCGATCCTGTAACTGAACAGGGGCGTCCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAg                        <  2:760790/93‑1 (MQ=255)
                                                                  aCATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACTGt                       >  1:430791/1‑93 (MQ=255)
                                                                  aCATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGt                       <  1:434557/93‑1 (MQ=255)
                                                                  aCATCGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGt                       <  1:3112302/93‑1 (MQ=255)
                                                                      cGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTgag                   >  1:4054515/1‑93 (MQ=255)
                                                                      cGATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTgag                   <  2:4175188/93‑1 (MQ=255)
                                                                       gATCCTGTAACTGAACTGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTgaga                  >  2:3030777/1‑93 (MQ=255)
                                                                       gATCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTgaga                  <  1:3675987/93‑1 (MQ=255)
                                                                        aTCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGaa                 >  2:878443/1‑93 (MQ=255)
                                                                         tCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGGCAGTTGAGAAc                <  1:100304/93‑1 (MQ=255)
                                                                         tCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAAc                <  1:740446/93‑1 (MQ=255)
                                                                         tCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAAc                <  2:3129897/93‑1 (MQ=255)
                                                                         tCCTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAAc                <  1:740447/93‑1 (MQ=255)
                                                                          ccTGTAACTGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGGCAGCAATGGACAGTTGAGaagg               >  1:2695074/1‑91 (MQ=255)
                                                                                  tGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCa       <  2:1230634/93‑1 (MQ=255)
                                                                                  tGAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCa       <  2:1230635/93‑1 (MQ=255)
                                                                                  ngAACAGGGGCGTTCAAGCCGNCTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCa       <  1:3852106/92‑1 (MQ=255)
                                                                                   gAACAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGCCCTCGCAt      <  2:811084/93‑1 (MQ=255)
                                                                                    aaCAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATc     <  2:3152276/93‑1 (MQ=255)
                                                                                     aCAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCg    >  1:790602/1‑93 (MQ=255)
                                                                                      cAGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGc   >  1:2546792/1‑93 (MQ=255)
                                                                                      cAGGAGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGc   <  2:3069552/93‑1 (MQ=255)
                                                                                      cAGGAGCGTTCAAGCCAACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGc   >  2:4508696/1‑93 (MQ=255)
                                                                                       aGGGGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGCt  >  2:1899255/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTAACCTGGCAGCTGCTAACACTCTTTTTGTTATGCGTCTGAACGACCGTGCGGGTGAAACGCGCTACATCGATCCTGTAACTGAACAGGAGCGTTCAAGCCGACTTTGGCTACGTCAAATTGGCGGGCATAATGCCTGGCGTGACAGCAACGGACAGTTGAGAACGACCTCGCATCGCT  >  NC_000913/1220520‑1220699

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