Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,609,260 G→A H207Y (CAT→TAT)  b1522 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,609,2600GA86.8% 691.1 / 76.6 250H207Y (CAT→TAT) b1522hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/17);  new base A (97/119);  total (113/137)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-02

Only 100 of 256 total aligned reads displayed.

CGACAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGA  >  NC_000913/1609170‑1609349
                                                                                          |                                                                                         
cGACAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATAcc                                                                                         <  2:6229256/93‑1 (MQ=255)
 gACAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCg                                                                                        <  2:429147/93‑1 (MQ=255)
  aCAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATGCCgg                                                                                       >  2:1795811/1‑93 (MQ=255)
  aCAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCgg                                                                                       >  2:1388379/1‑93 (MQ=255)
  aCAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCgg                                                                                       <  1:3734435/93‑1 (MQ=255)
   cAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCTCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCggg                                                                                      <  2:1189811/93‑1 (MQ=255)
     ggATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCa                                                                                    >  1:6147361/1‑93 (MQ=255)
          cACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAacac                                                                               >  1:5890257/1‑93 (MQ=255)
          cACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAacac                                                                               >  2:2693207/1‑93 (MQ=255)
          cACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAacac                                                                               <  2:1425510/93‑1 (MQ=255)
           aCCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACg                                                                              >  2:5713802/1‑93 (MQ=255)
            ccGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGc                                                                             <  1:3340798/93‑1 (MQ=255)
             cGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGTGCAAACACGCt                                                                            <  1:4222464/93‑1 (MQ=255)
             cGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCt                                                                            <  1:1037033/93‑1 (MQ=255)
               aCGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTa                                                                          >  1:5811757/1‑93 (MQ=255)
                 ggCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATGCCGGGCaaa                                                                                  <  2:238704/83‑1 (MQ=255)
                  gCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAaca                                                                       >  2:2575946/1‑93 (MQ=255)
                  gCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAaca                                                                       >  1:5030745/1‑93 (MQ=255)
                   cACTGCAACAGCAAATTCTTCGCCGCGCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAacac                                                                      >  1:4237168/1‑93 (MQ=255)
                   cACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAacac                                                                      <  1:6208162/93‑1 (MQ=255)
                     cTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACtt                                                                    >  1:2358180/1‑93 (MQ=255)
                      tGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACttt                                                                   <  2:3624752/93‑1 (MQ=255)
                         aaCAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATc                                                                <  1:2044015/93‑1 (MQ=255)
                             gCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCca                                                            <  2:2264561/93‑1 (MQ=255)
                             gCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTAAATATACCGGGCAAACACGCTTAACACTTTATCGCca                                                            <  2:1491309/93‑1 (MQ=255)
                               aaaTTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCcaca                                                          <  1:4819297/93‑1 (MQ=255)
                                aaTTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACAt                                                         <  2:6561071/93‑1 (MQ=255)
                                 aTTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACAtt                                                        <  1:5388747/93‑1 (MQ=255)
                                  ttcttcGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTc                                                       >  1:3965944/1‑93 (MQ=255)
                                   tcttcGCCGCCCATTCGCGCCCCCAGCCCCTTATCACCGTCAATCTTCTGACTCTACCGGGCCAACCCGCTTCACACTTTATCGCCACATTCa                                                      <  2:6648472/93‑1 (MQ=255)
                                     ttcGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATg                                                    >  1:5828263/1‑93 (MQ=255)
                                     ttcGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATg                                                    <  1:2968855/93‑1 (MQ=255)
                                        gccgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTcc                                                 >  2:1507251/1‑93 (MQ=255)
                                        gccgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTcc                                                 <  1:1580403/93‑1 (MQ=255)
                                        gccgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTcc                                                 >  1:721492/1‑93 (MQ=255)
                                        cccgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTcc                                                 <  2:1652187/92‑1 (MQ=255)
                                          cgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCAt                                               <  2:836/93‑1 (MQ=255)
                                          cgccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCAt                                               <  2:5933080/93‑1 (MQ=255)
                                           gccCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATa                                              <  1:438118/93‑1 (MQ=255)
                                            cccATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAg                                             <  2:4115215/93‑1 (MQ=255)
                                            cccATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAg                                             >  2:2326375/1‑93 (MQ=255)
                                            cccATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGACACATTCATGTCCATAg                                             >  2:1293377/1‑93 (MQ=255)
                                               aTTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTa                                          >  2:910372/1‑93 (MQ=255)
                                               aTTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTa                                          >  1:3638171/1‑93 (MQ=255)
                                                ttCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTAt                                         >  1:6351530/1‑93 (MQ=255)
                                                ttCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTAt                                         >  1:1154444/1‑93 (MQ=255)
                                                ttCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTAt                                         >  2:4897806/1‑93 (MQ=255)
                                                ttCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTAt                                         >  2:2821306/1‑93 (MQ=255)
                                                  cgcgCCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATct                                       >  1:5615159/1‑92 (MQ=255)
                                                   gcgcCACCAGCCCCTTATCTCCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGt                                      >  1:6465796/1‑93 (MQ=255)
                                                   gcgcCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGt                                      <  2:2638046/93‑1 (MQ=255)
                                                   gcgcCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGt                                      <  1:3035455/93‑1 (MQ=255)
                                                     gcCACCAGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTa                                    >  1:5086917/1‑93 (MQ=255)
                                                       caccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAAt                                  <  2:3081564/93‑1 (MQ=255)
                                                       caccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAAt                                  <  2:649175/93‑1 (MQ=255)
                                                        cccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATg                                 <  2:1620168/92‑1 (MQ=255)
                                                        cccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATg                                 <  2:1111695/92‑1 (MQ=255)
                                                        accaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATg                                 <  1:5239996/93‑1 (MQ=255)
                                                         ccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGc                                >  1:4448530/1‑93 (MQ=255)
                                                         ccaGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGc                                <  1:328211/93‑1 (MQ=255)
                                                           aGCCCCTTATCACCGACAATCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCtt                              >  1:283454/1‑93 (MQ=255)
                                                           aGCCCCTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCtt                              >  2:5577865/1‑93 (MQ=255)
                                                             ccccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCtttt                            <  1:1738367/93‑1 (MQ=255)
                                                             ccccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCtttt                            >  1:2555759/1‑93 (MQ=255)
                                                             ccccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCtttt                            <  2:5198280/93‑1 (MQ=255)
                                                              cccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTg                           <  2:2989698/93‑1 (MQ=255)
                                                              cccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTg                           >  2:3679164/1‑93 (MQ=255)
                                                               ccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGa                          >  1:5454350/1‑93 (MQ=255)
                                                               ccTTATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGa                          <  2:5821989/93‑1 (MQ=255)
                                                                 ttATCACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGaaa                        <  1:2124731/93‑1 (MQ=255)
                                                                     cACCGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGt                    <  1:5680247/93‑1 (MQ=255)
                                                                       ccGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCg                  >  2:4437559/1‑93 (MQ=255)
                                                                        cGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGa                 >  1:2640171/1‑93 (MQ=255)
                                                                        cGACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGa                 <  2:4061868/93‑1 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                <  2:318341/93‑1 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                >  2:3506909/1‑93 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                >  2:1978362/1‑93 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                >  2:6279718/1‑93 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                <  1:2390102/93‑1 (MQ=255)
                                                                         gACAATCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGat                >  2:4589373/1‑93 (MQ=255)
                                                                           cAATCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGatat              <  2:3019516/93‑1 (MQ=255)
                                                                            aaTCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATc             >  2:4873498/1‑93 (MQ=255)
                                                                            aaTCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATc             >  2:380283/1‑93 (MQ=255)
                                                                             aTCTTCTGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCa            >  1:135703/1‑93 (MQ=255)
                                                                              tcttctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCaa           >  1:1106690/1‑93 (MQ=255)
                                                                               cttctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAag          <  1:4913900/93‑1 (MQ=255)
                                                                               cttctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAag          >  1:715869/1‑93 (MQ=255)
                                                                                ttctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagc         <  2:2923609/93‑1 (MQ=255)
                                                                                ttctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagc         >  1:1296791/1‑93 (MQ=255)
                                                                                 tctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagca        >  1:2165626/1‑93 (MQ=255)
                                                                                 tctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGAGATCAagca        <  1:6491571/93‑1 (MQ=255)
                                                                                  ctGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagcag       <  1:5969218/93‑1 (MQ=255)
                                                                                  ctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagcag       <  2:1845813/93‑1 (MQ=255)
                                                                                  ctGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagcag       >  2:5025610/1‑93 (MQ=255)
                                                                                   tGAATATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAagcagc      >  1:4114184/1‑93 (MQ=255)
                                                                                    taaTATACCGGGCAAACACGCTTAACGCTTTATCGCCACATTCATGTCAATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAgcagca     <  1:3597252/92‑1 (MQ=255)
                                                                                     aaTATACCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCAt    <  2:332854/93‑1 (MQ=255)
                                                                                       tataCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGa  >  2:5269463/1‑93 (MQ=255)
                                                                                       tatGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGa  <  2:6112347/93‑1 (MQ=255)
                                                                                       tatGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGa  <  1:1329577/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGACAGGATTCACCGACGGCACTGCAACAGCAAATTCTTCGCCGCCCATTCGCGCCACCAGCCCCTTATCACCGACAATCTTCTGAATATGCCGGGCAAACACGCTTAACACTTTATCGCCACATTCATGTCCATAGTTATCGTTAATGCTTTTGAAATAGTCGATATCAAGCAGCATGA  >  NC_000913/1609170‑1609349

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