Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,326 G→A 17.2% A283V (GCC→GTC)  ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,3260GA17.2% 451.8 / 66.1 185A283V (GCC→GTC) ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (78/74);  new base A (16/16);  total (95/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 189 total aligned reads displayed.

CACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTC  >  NC_000913/3066227‑3066425
                                                                                                     |                                                                                                   
cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                                     >  1:4008215/1‑100 (MQ=255)
 aCGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGc                                                                                                    <  1:1727084/100‑1 (MQ=255)
 aCGTCTGCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTgg                                                                                                     <  2:3271963/99‑1 (MQ=255)
  cGTCTTCGACTTACATCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                                   <  1:4078922/100‑1 (MQ=21)
  cGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACg                                                                                                   <  1:1884120/100‑1 (MQ=21)
   gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCgg                                                                                                  >  1:294614/1‑100 (MQ=255)
   gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                                  <  1:1503368/100‑1 (MQ=21)
   gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACgg                                                                                                  <  1:2184041/100‑1 (MQ=21)
   gTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACGCGTCCCGCCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                                  <  2:1264173/100‑1 (MQ=21)
    tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                                 <  1:3976105/100‑1 (MQ=255)
    tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGt                                                                                                 >  1:3283264/1‑100 (MQ=255)
     cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                                <  1:2025460/100‑1 (MQ=255)
     cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                                <  2:1312905/100‑1 (MQ=255)
     cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                                >  1:2737514/1‑100 (MQ=255)
     cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                                <  2:1898727/100‑1 (MQ=255)
     cTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGta                                                                                                <  1:3080500/100‑1 (MQ=255)
       tCGACTTACACCGTAATTTAATTGCAACCGGCTCGGTGTTGGGTTAATGACTCGTTACGTCAAGAT‑CTTCAGCAATA‑TTGTGAAGACTCGTGTGGGtat                                                                                               >  2:1920790/1‑99 (MQ=255)
       tCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGGCTCGTCCCGGCACGATGCTTCAGCATAT‑TGCTGAAG‑ATCGTGACgg                                                                                                  <  2:511963/96‑1 (MQ=18)
         gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                            >  2:3185112/1‑100 (MQ=255)
         gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                            <  1:2382932/100‑1 (MQ=255)
         gACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                            <  2:2079165/100‑1 (MQ=255)
           cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAAc                                                                                            >  2:858951/1‑98 (MQ=255)
             tACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGATGAAGCATCGTGGCGGTATAACgatg                                                                                        >  1:1914640/1‑100 (MQ=255)
              acacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACgatga                                                                                       <  1:1396645/100‑1 (MQ=255)
                 ccGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                                    >  1:2866943/1‑100 (MQ=255)
                  cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGAttt                                                                                    >  2:3269239/1‑99 (MQ=255)
                  cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                   <  2:23334/100‑1 (MQ=255)
                   gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                  <  1:1027142/100‑1 (MQ=255)
                   gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                  <  2:1034246/100‑1 (MQ=255)
                   gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                  <  2:3286348/100‑1 (MQ=255)
                         taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGa                                                                            <  1:3575129/100‑1 (MQ=255)
                           attGCAACCGGCTCGATGCTGAGTTAATGGCTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                          <  2:2996491/100‑1 (MQ=255)
                              gCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                       >  1:3140840/1‑100 (MQ=255)
                               cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                                      >  1:637695/1‑100 (MQ=255)
                                aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                     <  1:2828108/100‑1 (MQ=255)
                                 aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat                                                                    <  1:3668662/100‑1 (MQ=255)
                                 aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGa                                                                     <  2:2559315/99‑1 (MQ=255)
                                  ccGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGata                                                                   <  1:1089485/100‑1 (MQ=255)
                                     gCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                >  1:1716766/1‑100 (MQ=255)
                                      cTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc                                                               <  1:2106768/100‑1 (MQ=255)
                                        cGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                             <  1:747045/100‑1 (MQ=255)
                                          aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                            <  2:902965/99‑1 (MQ=255)
                                          aTGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                           <  1:2597810/100‑1 (MQ=255)
                                             cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCatat                                                        <  1:3371864/100‑1 (MQ=255)
                                             cTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATGTCTGATTGAATAAAGATATTCCTGCatat                                                        >  2:1739484/1‑100 (MQ=255)
                                                   tAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                   >  1:478580/1‑99 (MQ=255)
                                                      atcgttataccGCCACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCAGc                                                >  2:446156/10‑99 (MQ=9)
                                                       gACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCt                                              >  2:2751220/1‑100 (MQ=255)
                                                            gTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                         >  1:1463637/1‑100 (MQ=255)
                                                             tCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                        <  1:333174/100‑1 (MQ=255)
                                                              cccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                       <  2:3696750/100‑1 (MQ=255)
                                                              cccGTCACGAT‑CTTCAGCAATA‑TGATGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATAGATCGGCCtgtg                                           >  2:19607/1‑96 (MQ=255)
                                                               ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGt                                       >  1:77580/1‑99 (MQ=255)
                                                               ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                      >  2:1024689/1‑100 (MQ=255)
                                                               ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                      >  1:2208635/1‑100 (MQ=255)
                                                               ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                      <  1:1021072/100‑1 (MQ=255)
                                                               ccGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTCCGGTg                                      >  2:3499516/1‑100 (MQ=255)
                                                                cgccACGATGCTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGc                                     >  2:4045511/4‑100 (MQ=21)
                                                                cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTg                                      <  2:2688319/99‑1 (MQ=255)
                                                                cGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGACATTCCTGCATATACATCGGCCTGTGCGGTGc                                     >  2:569619/1‑100 (MQ=255)
                                                                  tCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCAt                                   >  1:438568/1‑100 (MQ=255)
                                                                   cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                    <  2:4003109/98‑1 (MQ=255)
                                                                   cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                  >  1:2062196/1‑100 (MQ=255)
                                                                   cACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                  >  2:4097484/1‑100 (MQ=255)
                                                                    acgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                                 >  1:2298563/7‑100 (MQ=25)
                                                                     cgctgcTTCAGCATCT‑TGCTGAAG‑ATCGTGACGGTATAACCCTGATCTCTGATTGCCTACAGATATTCCTGCATATCCATCGGCCAGTGCGGCGCATCaa                                >  2:2234481/6‑100 (MQ=16)
                                                                     cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                >  1:1173122/6‑100 (MQ=25)
                                                                     cgatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                >  1:2636645/6‑100 (MQ=25)
                                                                     cGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                >  2:2115342/1‑100 (MQ=255)
                                                                      gatgcTTCAGC‑ATATTGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaa                               >  1:1637917/5‑100 (MQ=25)
                                                                       aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                              <  1:3370521/100‑1 (MQ=255)
                                                                       aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                              >  1:1136233/1‑100 (MQ=255)
                                                                       aT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTATGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                              >  2:287122/1‑100 (MQ=255)
                                                                           ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                           <  1:3911050/100‑1 (MQ=255)
                                                                           ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                           >  2:3952308/1‑100 (MQ=255)
                                                                           ttCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                           >  2:1988111/1‑100 (MQ=255)
                                                                            tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                          <  1:3175544/100‑1 (MQ=255)
                                                                            tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                          >  2:1121422/1‑100 (MQ=255)
                                                                            tCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                          >  1:944020/1‑100 (MQ=255)
                                                                            ctcagcatattGCTGAAG‑ATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                          <  1:3767176/90‑1 (MQ=38)
                                                                              cagcatattGCTGAAG‑ATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTc                        <  2:37359/92‑1 (MQ=255)
                                                                               gCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                       >  1:3483859/1‑100 (MQ=255)
                                                                                cAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAg                      <  1:372427/100‑1 (MQ=255)
                                                                                 aaTA‑TGATGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGc                     <  2:553688/100‑1 (MQ=255)
                                                                                    a‑tGCTGCAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCCGCGt                   <  2:138944/99‑1 (MQ=255)
                                                                                      tGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTg                 <  2:3561580/100‑1 (MQ=255)
                                                                                       gCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc                >  1:3349511/1‑100 (MQ=255)
                                                                                       gCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGc                >  1:3613408/1‑100 (MQ=255)
                                                                                        cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc               >  1:582666/1‑100 (MQ=255)
                                                                                        cTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGcc               >  2:2186158/1‑100 (MQ=255)
                                                                                         tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg              <  1:534538/100‑1 (MQ=255)
                                                                                         tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg              >  1:2770109/1‑100 (MQ=255)
                                                                                         tGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCg              <  2:1346017/100‑1 (MQ=255)
                                                                                          gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg             <  2:3842633/100‑1 (MQ=255)
                                                                                          gAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCgg             >  1:2477566/1‑100 (MQ=255)
                                                                                            aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGa            >  1:903651/1‑99 (MQ=255)
                                                                                            aGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt           <  1:514731/100‑1 (MQ=255)
                                                                                                 cGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAcc      <  1:1287791/100‑1 (MQ=255)
                                                                                                 cGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAcc      <  1:2039129/100‑1 (MQ=255)
                                                                                                    ggCGGTATAACGATGATTTCTGATTGAATAAAGATATTACTGCATATACATCGGCCTCTGAGGTGCATAAAACCATTACAGCGTTACCCGCTTCAcc      >  2:3644608/1‑97 (MQ=255)
                                                                                                     |                                                                                                   
CACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGAT‑CTTCAGCAATA‑TGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTC  >  NC_000913/3066227‑3066425

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