Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 311,355 T→A 5.3% intergenic (‑19/‑757) ecpR ← / → ykgL DNA‑binding transcriptional dual regulator MatA/uncharacterized protein YkgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913311,3550TA 5.3% 540.0 / 12.1 166intergenic (‑19/‑757)ecpR/ykgLDNA‑binding transcriptional dual regulator MatA/uncharacterized protein YkgL
Reads supporting (aligned to +/- strand):  ref base T (72/85);  new base A (5/4);  total (77/89)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 175 total aligned reads displayed.

ACAATCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATATGTGAAC  >  NC_000913/311206‑311504
                                                                                                                                                     |                                                                                                                                                     
acaATCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAt                                                                                                                                                       >  1:938872/1‑150 (MQ=255)
 caATCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAtt                                                                                                                                                      >  1:1217009/1‑150 (MQ=255)
    tCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAttttt                                                                                                                                                   <  1:1378439/150‑1 (MQ=255)
    tCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAttttt                                                                                                                                                   <  1:629550/150‑1 (MQ=255)
     cAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTa                                                                                                                                                  >  1:15499/1‑150 (MQ=255)
      aGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTAc                                                                                                                                                 >  1:195222/1‑150 (MQ=255)
      aGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTAc                                                                                                                                                 >  2:105008/1‑150 (MQ=255)
      aGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTAc                                                                                                                                                 <  2:136298/150‑1 (MQ=255)
       gTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACa                                                                                                                                                >  1:998107/1‑150 (MQ=255)
          cAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGt                                                                                                                                             >  1:1006567/1‑150 (MQ=255)
           aGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAAAATTACAGGtt                                                                                                                                            <  2:85402/150‑2 (MQ=255)
             aTAGTCCTTTTTAGAAGTATGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCGTACTGTAATCACTTTGCCATGTCACTACTTTCCAATCCTGTAATTTTTACAGGTGCa                                                                                                                                          <  2:1370708/150‑1 (MQ=255)
                gTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTc                                                                                                                                       <  1:26283/150‑1 (MQ=255)
                  ccTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAg                                                                                                                                     >  2:1381894/1‑150 (MQ=255)
                  ccTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAg                                                                                                                                     >  1:1407007/1‑150 (MQ=255)
                   cTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTAATCAGAACGGTTTTGACATTCCATATCGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGt                                                                                                                                    >  1:509118/1‑150 (MQ=255)
                      tataGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTc                                                                                                                                 >  1:288619/1‑150 (MQ=255)
                       ataGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTcc                                                                                                                                <  1:823007/150‑1 (MQ=255)
                        taGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCt                                                                                                                               >  2:958041/1‑150 (MQ=255)
                        taGAAGTAGGCGTCATGAGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCt                                                                                                                               >  2:1203896/1‑150 (MQ=255)
                         aGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCtt                                                                                                                              <  1:1053191/150‑1 (MQ=255)
                          gAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTa                                                                                                                             >  1:806705/1‑150 (MQ=255)
                          gAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTa                                                                                                                             >  1:708597/1‑150 (MQ=255)
                          gAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTa                                                                                                                             >  1:402944/1‑150 (MQ=255)
                           aaGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAg                                                                                                                            <  1:159325/150‑1 (MQ=255)
                            aGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGt                                                                                                                           >  1:18977/1‑150 (MQ=255)
                            aGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGt                                                                                                                           <  1:1354381/150‑1 (MQ=255)
                                   gTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCgg                                                                                                                    <  2:1401281/150‑1 (MQ=255)
                                   gTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCgg                                                                                                                    <  1:853851/150‑1 (MQ=255)
                                        ggggAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTa                                                                                                               <  1:415655/150‑1 (MQ=255)
                                           gagaCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAAt                                                                                                            <  1:1203896/150‑1 (MQ=255)
                                             gaCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATaa                                                                                                          <  2:195222/150‑1 (MQ=255)
                                              aCCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAg                                                                                                         <  1:1096305/150‑1 (MQ=255)
                                                cAGATGTATTTATCAGAAAGGTTTTGACATTCCATATGGTTTTTAACATCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGAt                                                                                                       >  2:616847/1‑150 (MQ=255)
                                                 ataTGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCACTACTGTAATCACTTTTCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCTGAAGTAAATAAGATa                                                                                                      <  2:1126236/148‑1 (MQ=255)
                                                   aTGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCGATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAATAAACCGGAAGTAAATAAGATACg                                                                                                    <  2:1415113/150‑1 (MQ=255)
                                                   aTGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACg                                                                                                    >  1:1254473/1‑150 (MQ=255)
                                                    tGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTAAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgt                                                                                                   <  2:939687/150‑1 (MQ=255)
                                                      tATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtag                                                                                                 <  2:882293/150‑1 (MQ=255)
                                                       aTTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACgtagt                                                                                                >  1:1128604/1‑150 (MQ=255)
                                                         ttATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCTGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTAc                                                                                              >  2:602385/1‑150 (MQ=255)
                                                          tATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGt                                                                                                                           >  2:489340/1‑120 (MQ=255)
                                                             cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACatat                                                                                          >  2:207248/1‑150 (MQ=255)
                                                             cAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAAAAATTACAGGtttgg                                                                                                                                         <  1:246107/103‑5 (MQ=255)
                                                               gAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCa                                                                                        >  1:522833/1‑150 (MQ=255)
                                                                   ggTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATaa                                                                                    >  2:237199/1‑150 (MQ=255)
                                                                    gTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATaaa                                                                                   >  2:767591/1‑150 (MQ=255)
                                                                           cATTCCATATGGTTTTTAACCTCATAGTCCCTTCTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTaaa                                                                            <  1:168335/150‑1 (MQ=255)
                                                                            aTTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAAc                                                                           <  2:1144480/150‑1 (MQ=255)
                                                                                cATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACgaga                                                                       <  2:138872/150‑1 (MQ=255)
                                                                                    tGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTaa                                                                   >  2:1184873/1‑150 (MQ=255)
                                                                                     ggTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAg                                                                                                                 >  1:1309994/1‑103 (MQ=255)
                                                                                     ggTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAg                                                                                                                 <  2:1309994/103‑1 (MQ=255)
                                                                                           tAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAat                                                            >  1:976885/1‑150 (MQ=255)
                                                                                           tAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAat                                                            >  2:59349/1‑150 (MQ=255)
                                                                                           tAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAat                                                            <  1:176738/150‑1 (MQ=255)
                                                                                             aCCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTAAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAatta                                                          <  2:382571/150‑1 (MQ=255)
                                                                                                tCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTAttt                                                       >  1:1248100/1‑150 (MQ=255)
                                                                                                tCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTAttt                                                       <  2:749734/150‑1 (MQ=255)
                                                                                                 cATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTAtttt                                                      >  2:330889/1‑150 (MQ=255)
                                                                                                  aTAGTCGCTACTGTAATCTCTTTTCGATGTCTCTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTTAATCGGAATTAAATAACATACTTAGTACATTTCATTTTAGGTTTATGGATATTAAGTTTTATTAttttt                                                     <  2:442381/150‑2 (MQ=255)
                                                                                                  aTAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTc                                                     <  2:864042/150‑1 (MQ=255)
                                                                                                    aGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGa                                                   <  2:402944/150‑1 (MQ=255)
                                                                                                      tCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGAtt                                                 >  1:1033624/1‑150 (MQ=255)
                                                                                                        ccTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTgg                                               <  2:422611/150‑1 (MQ=255)
                                                                                                        ccTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTgg                                               >  1:661489/1‑150 (MQ=255)
                                                                                                         cTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTggg                                              <  1:59349/150‑1 (MQ=255)
                                                                                                         cTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTggg                                              <  2:660457/150‑1 (MQ=255)
                                                                                                              gTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTaa                                         >  2:700501/1‑150 (MQ=255)
                                                                                                                  tCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAAttt                                     <  2:15499/150‑1 (MQ=255)
                                                                                                                   cACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTa                                    >  2:219222/1‑150 (MQ=255)
                                                                                                                   cACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTa                                    <  2:806705/150‑1 (MQ=255)
                                                                                                                    aCTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTAt                                   >  2:1256708/1‑150 (MQ=255)
                                                                                                                     cTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATg                                  >  1:431958/1‑150 (MQ=255)
                                                                                                                     cTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATg                                  <  2:1248100/150‑1 (MQ=255)
                                                                                                                         gCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttt                              <  2:69373/150‑1 (MQ=255)
                                                                                                                           cATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttttt                            <  1:85841/150‑1 (MQ=255)
                                                                                                                           cATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGtttttt                            >  2:690977/1‑150 (MQ=255)
                                                                                                                            aTGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGttttttt                           <  1:200445/150‑1 (MQ=255)
                                                                                                                              gTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTc                         <  2:288619/150‑1 (MQ=255)
                                                                                                                              gTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTc                         <  1:466009/150‑1 (MQ=255)
                                                                                                                               tCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTcc                        >  1:565153/1‑150 (MQ=255)
                                                                                                                               tCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTcc                        <  2:250025/150‑1 (MQ=255)
                                                                                                                               tCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTcc                        <  2:186778/150‑1 (MQ=255)
                                                                                                                                cACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTccc                       <  2:1418074/150‑1 (MQ=255)
                                                                                                                                cACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTccc                       <  2:526529/150‑1 (MQ=255)
                                                                                                                                       ttCCAAACCTGTAATTTTTACAGGTGCATTCTGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgt                <  1:602385/150‑1 (MQ=255)
                                                                                                                                       ttCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgt                >  1:1352414/1‑150 (MQ=255)
                                                                                                                                       gaatgcaCCTGTAAAAATTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgt                >  1:692332/7‑150 (MQ=255)
                                                                                                                                        tCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgtg               >  2:1151385/1‑150 (MQ=255)
                                                                                                                                         atgcaCCTGTAAAAATTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAgtgt              >  1:750642/5‑150 (MQ=255)
                                                                                                                                          cAAACCTGTAATTTTTACATGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAACTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTc             <  1:964727/150‑1 (MQ=255)
                                                                                                                                            caCCTGTAAAAATTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCaa           >  2:893830/2‑150 (MQ=255)
                                                                                                                                            aaCCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCaa           >  2:182019/1‑150 (MQ=255)
                                                                                                                                             aCCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAat          >  2:1431186/1‑150 (MQ=255)
                                                                                                                                                tGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATAtg       <  1:1287066/150‑1 (MQ=255)
                                                                                                                                                tGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATAtg       >  1:732594/1‑150 (MQ=255)
                                                                                                                                                 gTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATAtgt      >  2:311996/1‑150 (MQ=255)
                                                                                                                                                   aaTTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGCTTAAATTTATGTTTTTTTTCCAGATGAGTGTCAATATGTGaa   >  1:719542/1‑150 (MQ=255)
                                                                                                                                                     tttttACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATATGTGAAc  >  1:1023129/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAATCAGTTCAGATAGTCCTTTATAGAAGTAGGCGTCATGGGGAGACCAGATGTATTTATCAGAACGGTTTTGACATTCCATATGGTTTTTAACCTCATAGTCCCTACTGTAATCACTTTGCCATGTCACTACTTTCCAAACCTGTAATTTTTACAGGTGCATTCAGTTTCCTTAGTAAACCGGAAGTAAATAAGATACGTAGTACATATCAATAAAGGTTAAACGAGATTAAGTTTAATTATTTTCGATTGGGCTTAAATTTATGTTTTTTTTCCCGATGAGTGTCAATATGTGAAC  >  NC_000913/311206‑311504

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: