Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,587,936 A→G R30R (CGT→CGC b1505 ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,9360AG99.7% 1020.5 / ‑5.1 318R30R (CGT→CGCb1505putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (148/169);  total (149/169)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 331 total aligned reads displayed.

TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025
                                                                                          |                                                                                         
tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                         >  1:4979113/1‑93 (MQ=255)
tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                         >  2:4473263/1‑93 (MQ=255)
  aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCt                                                                                       <  1:6091326/93‑1 (MQ=255)
   cGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCtt                                                                                      <  1:9030932/93‑1 (MQ=255)
    gtgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTa                                                                                     >  1:1655252/1‑93 (MQ=255)
    gtgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTa                                                                                     <  1:2111899/93‑1 (MQ=255)
     tgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTAc                                                                                    >  1:572995/1‑93 (MQ=255)
      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   >  1:5921176/1‑93 (MQ=255)
      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   >  1:2513994/1‑93 (MQ=255)
       ttCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTg                                                                                  <  2:9904569/93‑1 (MQ=255)
       ttCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCGGACCGCGCTTACTg                                                                                  <  1:5429638/93‑1 (MQ=255)
        tCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGt                                                                                 >  1:5184222/1‑93 (MQ=255)
        tCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGt                                                                                 <  1:6025899/93‑1 (MQ=255)
         ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                <  2:3975002/93‑1 (MQ=255)
          cAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAg                                                                               <  2:9832055/93‑1 (MQ=255)
          cAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAg                                                                               <  1:76016/93‑1 (MQ=255)
            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                             >  1:9833749/1‑93 (MQ=255)
            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                             <  1:8247102/93‑1 (MQ=255)
             tAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAAc                                                                            <  1:4426004/93‑1 (MQ=255)
              aaCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACa                                                                           <  1:8169981/93‑1 (MQ=255)
              aaCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACa                                                                           <  2:4806071/93‑1 (MQ=255)
                ttttGGCGACTGGCTCTGAAAAATGTCGTAGCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  2:4109906/92‑1 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGNCCGCGCTTACTGTAGAACAgt                                                                         >  1:7807490/1‑92 (MQ=255)
                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACCGa                                                                         >  1:4044892/1‑93 (MQ=255)
                 tttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAt                                                                        >  1:1501016/1‑93 (MQ=255)
                  ttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAtt                                                                       <  2:5611770/93‑1 (MQ=255)
                  ttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAtt                                                                       <  2:6148327/93‑1 (MQ=255)
                    ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                     >  2:118070/1‑93 (MQ=255)
                     gCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGCGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaa                                                                    >  2:5783214/1‑93 (MQ=255)
                     gCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTCCTGTAGAACAGATTaaa                                                                    <  2:10041610/93‑1 (MQ=255)
                      cGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaaa                                                                   <  1:502660/93‑1 (MQ=255)
                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  2:6823251/93‑1 (MQ=255)
                        aCTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaata                                                                 <  2:3667485/93‑1 (MQ=255)
                          tGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataat                                                               <  2:5477311/93‑1 (MQ=255)
                          tGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataat                                                               >  2:9641875/1‑93 (MQ=255)
                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                              >  1:324811/1‑93 (MQ=255)
                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                              <  2:7352667/93‑1 (MQ=255)
                             cACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCaa                                                            <  1:901325/93‑1 (MQ=255)
                             cACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCaa                                                            >  2:3416309/1‑93 (MQ=255)
                             cACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCaa                                                            >  1:2887634/1‑93 (MQ=255)
                               cTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATa                                                          <  1:632066/93‑1 (MQ=255)
                               cTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATa                                                          >  2:2869212/1‑93 (MQ=255)
                               cTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATa                                                          >  2:1199858/1‑93 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        <  2:4407468/93‑1 (MQ=255)
                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        >  1:699454/1‑93 (MQ=255)
                                    aaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtg                                                     <  1:2790210/93‑1 (MQ=255)
                                    aaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTGATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtg                                                     <  2:19842/93‑1 (MQ=255)
                                      aTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                   <  1:3629917/93‑1 (MQ=255)
                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  2:2101571/93‑1 (MQ=255)
                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 >  2:7513177/1‑93 (MQ=255)
                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 >  1:3004674/1‑93 (MQ=255)
                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  2:7388458/93‑1 (MQ=255)
                                         tCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCt                                                <  1:896314/93‑1 (MQ=255)
                                             aCCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTc                                            <  2:2746330/93‑1 (MQ=255)
                                             aCCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTc                                            >  1:4212231/1‑93 (MQ=255)
                                              ccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCa                                           <  1:5247913/93‑1 (MQ=255)
                                              ccGTAATTTCCAAGATGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTACAATAATCAATCAATTGTGTCTGCTCa                                           >  1:5793727/1‑93 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         <  2:8183458/93‑1 (MQ=255)
                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         <  2:7573176/93‑1 (MQ=255)
                                                 tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                        >  1:10023366/1‑93 (MQ=255)
                                                 tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                        >  1:5669100/1‑93 (MQ=255)
                                                 tAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTg                                        >  1:6442368/1‑93 (MQ=255)
                                                     ttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtgt                                    <  1:1384189/93‑1 (MQ=255)
                                                     ttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtgt                                    <  1:4650720/93‑1 (MQ=255)
                                                      tCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtc                                   >  1:8346736/1‑93 (MQ=255)
                                                      tCCAAGCTGCTGAGAGATGCTTATTTGTCCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtc                                   <  2:890001/93‑1 (MQ=255)
                                                       ccAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtct                                  <  2:1891266/93‑1 (MQ=255)
                                                        cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                 >  2:7942607/1‑93 (MQ=255)
                                                        cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                 >  1:8847045/1‑93 (MQ=255)
                                                         aaGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTcc                                >  1:7882433/1‑93 (MQ=255)
                                                         aaGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTcc                                <  1:744104/93‑1 (MQ=255)
                                                           gctgctGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCgg                              >  1:9673592/1‑93 (MQ=255)
                                                           gctgctGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCgg                              >  1:1721628/1‑93 (MQ=255)
                                                           gctgctGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCgg                              <  1:9461526/93‑1 (MQ=255)
                                                            ctgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGt                             <  2:6240869/93‑1 (MQ=255)
                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            >  2:952576/1‑93 (MQ=255)
                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            >  2:7136817/1‑93 (MQ=255)
                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            <  1:9304657/93‑1 (MQ=255)
                                                              gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTag                           >  2:8048326/1‑93 (MQ=255)
                                                              gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGGGTCTCCGGTag                           <  2:432616/93‑1 (MQ=255)
                                                               ctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagg                          <  2:9942336/93‑1 (MQ=255)
                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         >  2:3499124/1‑93 (MQ=255)
                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         <  1:7211264/93‑1 (MQ=255)
                                                                 gagagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTaggag                        >  1:1181676/1‑93 (MQ=255)
                                                                 gagagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTaggag                        <  2:558901/93‑1 (MQ=255)
                                                                  agagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTaggagg                       >  2:3537792/1‑93 (MQ=255)
                                                                    agaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCTGGTAGGAGGCt                     <  1:3846127/93‑1 (MQ=255)
                                                                     gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAGATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                    >  2:6258421/1‑93 (MQ=255)
                                                                     gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                    <  1:4571390/93‑1 (MQ=255)
                                                                       tGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGa                  >  2:6409840/1‑93 (MQ=255)
                                                                       tGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGa                  <  1:7278638/93‑1 (MQ=255)
                                                                         cTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACg                <  2:1732385/93‑1 (MQ=255)
                                                                           tATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgt              <  1:8056127/93‑1 (MQ=255)
                                                                           tATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTATCTCCGGTAGGAGGCTTGACGgt              >  2:5522987/1‑93 (MQ=255)
                                                                            aTTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgtg             <  1:6000284/93‑1 (MQ=255)
                                                                              ttGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTa           >  2:7785284/1‑93 (MQ=255)
                                                                                gTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAc         <  2:1258089/93‑1 (MQ=255)
                                                                                gTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAc         <  1:909772/93‑1 (MQ=255)
                                                                                gTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAc         <  1:2624286/93‑1 (MQ=255)
                                                                                    cTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACt     >  1:2210559/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025

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