Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,439,588 G→C 16.0% intergenic (‑79/+28) yhdN ← / ← rplQ DUF1992 family protein/50S ribosomal subunit protein L17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,439,5880GC16.0% 288.2 / 45.4 113intergenic (‑79/+28)yhdN/rplQDUF1992 family protein/50S ribosomal subunit protein L17
Reads supporting (aligned to +/- strand):  ref base G (45/50);  new base C (15/3);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 8.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 127 total aligned reads displayed.

CACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCG  >  NC_000913/3439490‑3439685
                                                                                                  |                                                                                                 
cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                  >  1:827706/1‑100 (MQ=255)
cACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                  <  2:817649/100‑1 (MQ=255)
cACTGGTCAAGTAACCACATAAATACTCCAAAGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgg                                                                                                  <  1:332505/100‑1 (MQ=255)
 aCTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCggg                                                                                                 >  1:1501155/1‑100 (MQ=255)
   tGGTCAAGTAACAACATAAATACTGCAGGGGATGAGCAAAAAAAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCgggg                                                                                                <  1:2172334/99‑1 (MQ=255)
    ggTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCg                                                                                              >  2:1870166/1‑100 (MQ=255)
      tCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAGAAACCCGCCGGGGCggg                                                                                            <  1:744189/100‑1 (MQ=255)
       cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                            >  1:854392/1‑99 (MQ=255)
       cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCggg                                                                                            >  2:520242/1‑99 (MQ=255)
       cAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGt                                                                                           <  2:1852761/100‑1 (MQ=255)
         aGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttt                                                                                         >  2:57442/1‑100 (MQ=255)
            aaCCACATAAATCCTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCCCGGCGGGtttttt                                                                                      <  1:1772845/100‑1 (MQ=255)
            aaCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGtttttt                                                                                      <  2:2122575/100‑1 (MQ=255)
             aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTGTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGttttttt                                                                                     <  1:1813360/100‑1 (MQ=255)
             aCCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGAATAAAAAAATCCGCAGGGGCGGGtctttt                                                                                      >  1:1275179/1‑99 (MQ=255)
              ccACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTa                                                                                    <  2:878483/100‑1 (MQ=255)
               cacaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAATCCGCCCCGGCGGGTTTTTTTat                                                                                   <  1:1506437/100‑2 (MQ=39)
                acaTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACg                                                                                  >  1:1794017/1‑100 (MQ=255)
                 caTAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGt                                                                                 >  1:2477342/1‑100 (MQ=255)
                   tAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTg                                                                               <  2:609479/100‑1 (MQ=255)
                           ccAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAACAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAt                                                                       <  2:1158386/100‑1 (MQ=255)
                            cAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGAtt                                                                      >  2:1713213/1‑100 (MQ=255)
                                  aTGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                <  1:1561557/100‑1 (MQ=255)
                                  aTGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTg                                                                <  1:1561763/100‑1 (MQ=255)
                                   tGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                               >  2:2221480/1‑100 (MQ=255)
                                     naCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                             >  1:1944798/2‑100 (MQ=255)
                                     aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcag                                                             <  2:1959768/100‑1 (MQ=255)
                                     aaCAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgc                                                               <  1:1619272/98‑1 (MQ=255)
                                      aCAAAAAGGGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagc                                                            <  1:2023297/100‑1 (MQ=255)
                                       cAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTgcagca                                                           >  2:678208/1‑100 (MQ=255)
                                        aaaaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGcagcag                                                          >  1:2187547/1‑100 (MQ=255)
                                          aaaGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCt                                                        >  2:2238696/1‑100 (MQ=255)
                                            agagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTGCAGAGTAATCTgaagcaacgta                                                      <  2:170239/100‑6 (MQ=38)
                                             gagTACAGCTATTGTAGCTAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCATAAGATTACTCTGCAGCAGCTTCt                                                     >  1:580286/1‑100 (MQ=255)
                                             gagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCt                                                     <  2:291925/100‑1 (MQ=255)
                                             gagTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTGCAGAGTAATCTgaagcaacgtaa                                                     <  1:1633417/100‑7 (MQ=37)
                                              agTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTg                                                    <  2:1481803/100‑1 (MQ=255)
                                               gTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGc                                                   <  2:230187/100‑1 (MQ=255)
                                                tACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCt                                                  <  2:2450161/100‑1 (MQ=255)
                                                 aCAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCtt                                                 <  1:49825/100‑1 (MQ=255)
                                                  cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                >  1:970791/1‑100 (MQ=255)
                                                  cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                >  1:2389959/1‑100 (MQ=255)
                                                  cAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCttt                                                >  2:1687665/1‑100 (MQ=255)
                                                     cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                             <  1:1102998/100‑1 (MQ=255)
                                                     cTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTtctc                                             <  1:744811/100‑1 (MQ=255)
                                                       aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTctct                                            <  1:614982/99‑1 (MQ=255)
                                                       aTTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTg                                           >  2:490578/1‑100 (MQ=255)
                                                          gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCTGGGGGGGGGATGTTTACGTTTCTTCGGATTACCCGCAAGAAGCTTCCGCTTTCTCTGa                                          >  1:145651/1‑98 (MQ=255)
                                                          gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                        <  2:102944/100‑1 (MQ=255)
                                                          gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                        >  2:741335/1‑100 (MQ=255)
                                                          gTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAAc                                        >  1:1221098/1‑100 (MQ=255)
                                                             gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                     <  1:480259/100‑1 (MQ=255)
                                                             gATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGCGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGAt                                     <  1:2132190/100‑1 (MQ=255)
                                                               tAAGTGGGGATACTACGGGTATAAAAAACCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCa                                   <  1:1380740/100‑1 (MQ=255)
                                                               tAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCCGCAGCAGCCTCTGCTTTTTCTTAACGATc                                    >  1:2532488/1‑99 (MQ=255)
                                                                 aGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCTCAAc                                 >  2:1031168/1‑100 (MQ=255)
                                                                  gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                >  2:917275/1‑100 (MQ=255)
                                                                  gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                >  2:1533641/1‑100 (MQ=255)
                                                                  gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAcc                                <  1:283676/100‑1 (MQ=255)
                                                                  gTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCACATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGCTCAAcc                                >  1:912041/1‑100 (MQ=255)
                                                                    ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                              >  2:90707/1‑100 (MQ=255)
                                                                    ggggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGGTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAg                              >  2:600462/1‑100 (MQ=255)
                                                                     gggATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGc                             <  2:2021161/100‑1 (MQ=255)
                                                                      ggATACTACGTGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCt                            <  1:597777/100‑1 (MQ=255)
                                                                        aTACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCg                          <  2:49169/100‑1 (MQ=255)
                                                                           ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                       >  2:363219/1‑100 (MQ=255)
                                                                           ctacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATg                       <  2:1630913/100‑1 (MQ=255)
                                                                            tacGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGt                      <  1:1737660/100‑1 (MQ=255)
                                                                             acGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTa                     <  2:1980637/100‑1 (MQ=255)
                                                                                ggTATAAAAAGACCCGCCGGGGCGGGGTGTTTTACGTTGCTTCAGATTGCTCTGCAGCAGCTTCTGCTTTCTCTGAACGGTCAACCAGCTCGATGTAAGc                  <  1:69144/100‑1 (MQ=255)
                                                                                ggTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGc                  <  2:1271825/100‑1 (MQ=255)
                                                                                 gTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                 <  2:1221102/100‑1 (MQ=255)
                                                                                 acgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                 >  1:2028644/4‑100 (MQ=38)
                                                                                 acgtaAAAAAACCCGCCCCGGCGGGTGTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGcc                 >  1:1832483/4‑100 (MQ=38)
                                                                                  cgtaAAAAAACCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGNAAGCCa                >  2:2104225/3‑100 (MQ=255)
                                                                                    taAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATc              >  2:291535/1‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg             >  2:1774221/1‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg             >  1:1086764/1‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCg             >  2:1384776/1‑100 (MQ=255)
                                                                                     aaaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCGTTCTCTGAACGGTCAACCAGCTCGATGTAAGCCATCg             >  1:857865/1‑100 (MQ=38)
                                                                                      aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            >  1:1859052/1‑100 (MQ=255)
                                                                                      aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            >  2:1579470/1‑100 (MQ=255)
                                                                                      aaaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            >  2:2562963/1‑100 (MQ=255)
                                                                                      aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAAGCAGCTCGATGTAAGCCAtag             >  1:649215/1‑97 (MQ=38)
                                                                                      aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            >  1:2557182/1‑100 (MQ=255)
                                                                                      aaaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            >  2:85159/1‑100 (MQ=255)
                                                                                      aaaaaaCCAGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCgg            <  1:1747849/100‑1 (MQ=255)
                                                                                       aaaaaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc           <  2:523932/100‑1 (MQ=255)
                                                                                       aaaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgc           >  2:2180041/1‑100 (MQ=255)
                                                                                        aaaaCCCGCCCCGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGgcg          >  2:255992/1‑100 (MQ=255)
                                                                                          aaCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg        <  2:847045/100‑1 (MQ=255)
                                                                                           aCCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt       >  2:174786/1‑100 (MQ=255)
                                                                                            cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGcgcg        <  1:826793/98‑1 (MQ=255)
                                                                                            cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGtt      >  2:1402830/1‑100 (MQ=255)
                                                                                            cccGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGt       <  2:1274389/99‑1 (MQ=255)
                                                                                              cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGt    >  2:504129/1‑100 (MQ=255)
                                                                                              cgccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGt    >  2:379139/1‑100 (MQ=255)
                                                                                               gccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTc   >  2:2020120/1‑100 (MQ=255)
                                                                                                ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg  <  2:1222821/100‑1 (MQ=255)
                                                                                                ccgGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCg  <  1:2082965/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
CACTGGTCAAGTAACCACATAAATACTCCAGGGGATGAACAAAAAGAGTACAGCTATTGTAGATAAGTGGGGATACTACGGGTATAAAAAAACCCGCCGGGGCGGGTTTTTTTACGTTGCTTCAGATTACTCTGCAGCAGCTTCTGCTTTCTCTGAACGATCAACCAGCTCGATGTAAGCCATCGGCGCGTTGTCG  >  NC_000913/3439490‑3439685

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: