Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,667,825 A→C W8G (TGG→GGG)  hcaT ← putative 3‑phenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,667,8250AC77.2% 140.5 / 60.7 79W8G (TGG→GGG) hcaTputative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base C (36/25);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

GGCGACAGCAAAGAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGGCGATGAAATCGCGTAATG  >  NC_000913/2667589‑2668061
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
ggCGACAGCAAAGAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCgc                                                                                                                                                                                                                                            >  2:242220/1‑239 (MQ=255)
    aCAGCAAAGAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGt                                                                                                                                                                                                                                        >  1:202952/1‑239 (MQ=255)
                       gcagTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCCCCAGAAAGGGAGAAAAATGCCGTAACTAAAAAAGGATGTGAAATAGCCGAGCGCCAACCCGCGCGAGGATTGCAAAACCATGAg                                                                                                                                                                                                                      >  1:1105730/1‑238 (MQ=255)
                         aggggcAGCACGCGCAAGGCGGAAATCAGGGGGGAAGGGTCCCTGACGCGGGGCCCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGGCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCa                                                                                                                                                                                                                          <  2:1134078/227‑1 (MQ=255)
                                     cgcAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCagga                                                                                                                                                                                                       >  2:1038196/1‑239 (MQ=255)
                                     cgcAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCagga                                                                                                                                                                                                       >  1:271981/1‑239 (MQ=255)
                                         aGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCg                                                                                                                                                                                                   <  2:202952/239‑1 (MQ=255)
                                         aGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGGCAGGAGGCg                                                                                                                                                                                                   <  1:199427/239‑1 (MQ=255)
                                                aaTCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCt                                                                                                                                                                                            <  2:588249/239‑1 (MQ=255)
                                                                                      gCAAACGCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCACAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGTGGCAAGTTGAAAAATAAGa                                                                                                                                                      >  2:579092/1‑239 (MQ=255)
                                                                                        aaaCGCCCGAGGAAACGGGCAACCAGCCCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAGCCCGCGCGTGGATTGCAAAACCATGAGTCCCGCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGa                                                                                                                                                      <  1:238430/237‑1 (MQ=255)
                                                                                          acgcccGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAatat                                                                                                                                                  >  2:278493/4‑239 (MQ=255)
                                                                                                        gggCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTc                                                                                                                                     >  2:444450/1‑238 (MQ=255)
                                                                                                           cAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCg                                                                                                                                  <  1:689923/238‑1 (MQ=255)
                                                                                                            aaCCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGAt                                                                                                                                <  1:1038196/239‑1 (MQ=255)
                                                                                                             aCCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATg                                                                                                                               <  2:1148820/239‑1 (MQ=255)
                                                                                                              ccAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTccc                                                                                                                                   >  2:139450/1‑234 (MQ=255)
                                                                                                                   ccTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATg                                                                                                                               <  2:1083748/233‑1 (MQ=255)
                                                                                                                        ccccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGc                                                                                                                     <  1:242220/238‑1 (MQ=255)
                                                                                                                        ccccAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGGGCCTGCTTCCCGATGTCTGTAACGc                                                                                                                     <  1:61213/238‑1 (MQ=255)
                                                                                                                                    gCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGc                                                                                                         >  1:557551/1‑238 (MQ=255)
                                                                                                                                    gCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGt                                                                                                                                                                                                            >  2:305123/1‑139 (MQ=255)
                                                                                                                                    gCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGt                                                                                                                                                                                                            <  1:305123/139‑1 (MQ=255)
                                                                                                                                     ccGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTcaca                                                                                                            <  1:139450/234‑1 (MQ=255)
                                                                                                                                     ccGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAga                                                                                                         >  2:593086/1‑236 (MQ=255)
                                                                                                                                     ccGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAg                                                                                                          <  1:444450/236‑1 (MQ=255)
                                                                                                                                       gATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTa                                                                                                     <  1:716095/239‑1 (MQ=255)
                                                                                                                                        aTGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAACATTCGAGCCTGCTTCCCGATGTCTGTAACGCAAGTTCACAAGCACTaa                                                                                                    >  1:565726/1‑239 (MQ=255)
                                                                                                                                        aTGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTaa                                                                                                    >  2:904773/1‑239 (MQ=255)
                                                                                                                                                  gCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTgaga                                                                                          >  2:447817/1‑239 (MQ=255)
                                                                                                                                                      tAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATcc                                                                                       >  2:922545/1‑238 (MQ=255)
                                                                                                                                                             aTCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCa                                                                               >  1:958422/1‑239 (MQ=255)
                                                                                                                                                                   ttAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGc                                                                         >  1:1179946/1‑239 (MQ=255)
                                                                                                                                                                          aGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaata                                                                  >  1:298222/1‑239 (MQ=255)
                                                                                                                                                                            aCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGt                                                                                                                                                                                              >  2:989655/1‑113 (MQ=255)
                                                                                                                                                                            aCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGt                                                                                                                                                                                              <  1:989655/113‑1 (MQ=255)
                                                                                                                                                                            aCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaataat                                                                >  2:864185/1‑239 (MQ=255)
                                                                                                                                                                            aCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaataat                                                                >  1:9248/1‑239 (MQ=255)
                                                                                                                                                                            aCGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaataat                                                                >  2:1159080/1‑239 (MQ=255)
                                                                                                                                                                             cGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaata                                                                  <  2:958422/236‑1 (MQ=255)
                                                                                                                                                                             cGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAaata                                                                  <  2:1105730/236‑1 (MQ=255)
                                                                                                                                                                                   aGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAAGGCCAGTTCACAAGCACTAAAAAACTGAGATCAGAGACGCATGTTGCAAAAATAATACGTAAt                                                         >  2:1045936/1‑239 (MQ=255)
                                                                                                                                                                                      aaagTAAAAAAATGGCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTa                                                       <  1:435362/235‑1 (MQ=255)
                                                                                                                                                                                        ggTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATg                                                     >  2:622446/1‑238 (MQ=255)
                                                                                                                                                                                        ggTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATg                                                     >  1:645478/1‑238 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGt                                                                                                                          <  1:1043388/167‑1 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGt                                                                                                                          >  2:1043388/1‑167 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTa                                                       <  1:411977/234‑1 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTa                                                       >  2:411977/1‑234 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTAt                                                      >  2:1072234/1‑235 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTAt                                                      >  1:1173586/1‑235 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTAt                                                      >  2:1115635/1‑235 (MQ=255)
                                                                                                                                                                                          tAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAACATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTa                                                       <  2:565726/234‑1 (MQ=255)
                                                                                                                                                                                          aaGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACAATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAACAATAAGAATATTCGAGCCTGCTTCCCGATGTCCGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGAGACGCATGTCGCAACAATAAGGCGTAATTAt                                                      >  2:250499/2‑235 (MQ=255)
                                                                                                                                                                                           aGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCg                                                  <  1:278493/238‑1 (MQ=255)
                                                                                                                                                                                                aaTGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGc                                            <  2:271981/239‑1 (MQ=255)
                                                                                                                                                                                                   gCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCtt                                          >  1:957473/1‑238 (MQ=255)
                                                                                                                                                                                                   gCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCtt                                          >  2:111300/1‑238 (MQ=255)
                                                                                                                                                                                                   gCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCt                                                                                                                                                                                            <  1:1052537/92‑1 (MQ=255)
                                                                                                                                                                                                   gCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCt                                                                                                                                                                                            >  2:1052537/1‑92 (MQ=255)
                                                                                                                                                                                                   gCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCtt                                          >  2:385304/1‑238 (MQ=255)
                                                                                                                                                                                                       tAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATTCCGTTACGCTTGCCaa                                     >  1:924020/1‑239 (MQ=255)
                                                                                                                                                                                                        aaCTAAAAAAGTATGTGAAATAGCCGGGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCaaa                                    <  1:460969/239‑1 (MQ=255)
                                                                                                                                                                                                        aaCTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCaaa                                    <  1:1072234/239‑1 (MQ=255)
                                                                                                                                                                                                        aaCTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCaaa                                    <  2:645478/239‑1 (MQ=255)
                                                                                                                                                                                                        aaCTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCaaa                                    <  1:286289/239‑1 (MQ=255)
                                                                                                                                                                                                             aaaaaGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCagg                                                                                                                                                                                                        <  1:200258/70‑1 (MQ=255)
                                                                                                                                                                                                             aaaaaGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCagg                                                                                                                                                                                                        >  2:200258/1‑70 (MQ=255)
                                                                                                                                                                                                             aaaaaGTATGTGAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTc                               <  2:1179946/239‑1 (MQ=255)
                                                                                                                                                                                                                     tgtgAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAg                                                                                                          >  1:96199/1‑156 (MQ=255)
                                                                                                                                                                                                                     tgtgAAATAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAg                                                                                                          <  2:96199/156‑1 (MQ=255)
                                                                                                                                                                                                                     tgtgAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGt                                                                                                                                                                                              <  1:1056581/72‑1 (MQ=255)
                                                                                                                                                                                                                     tgtgAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGt                                                                                                                                                                                              >  2:1056581/1‑72 (MQ=255)
                                                                                                                                                                                                                            tAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATg                                                               >  2:1050178/1‑192 (MQ=255)
                                                                                                                                                                                                                            tAGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATg                                                               <  1:1050178/192‑1 (MQ=255)
                                                                                                                                                                                                                             aGCCGAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGTCGATg               >  2:826844/1‑239 (MQ=255)
                                                                                                                                                                                                                             aGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGGCGATg               >  2:1007128/1‑239 (MQ=255)
                                                                                                                                                                                                                                 gAGCGCCAACCCGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGGCGATGaaa            >  2:234550/1‑238 (MQ=255)
                                                                                                                                                                                                                                          ccAGCGCTTGGATTGCAAAACCATGATTCCCTCCCGTCAGGAAGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGGTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGGCGATGAAATCGCGTAATg  <  2:261559/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                            
GGCGACAGCAAAGAGAAGTGTCAGCAGTGCCAGCACGCGCAAGGCGGAAATCAGGCGGGAAGGATCGCTGACGCGGGGCGCGATGAGCAAACTCCCGAGGAAACGGGCAACCAGACCTGCCCCCAATAACAGGCCGATGGTTTCTGGCGTTAAACCAATCCCTTTAAGCCAGACGCTCCAGAAAGGTAGAAAAATGCCGTAACTAAAAAAGTATGTGAAATAGCCGAGCGCCAACCAGCGCGTGGATTGCAAAACCATGAGTCCCTCCCGTCAGGAGGCGTTAGTCTGGGGTTAATCGAATGTTGGTGCAAGTTGAAAAATAAGAATATTCGAGCCTGCTTCCCGATGTCTGTAACGCCAGTTCACAAGCACTAAAAAACTGAGATCCGCGACGCATGTTGCAAAAATAATGCGTAATTATGCCGTTACGCTTGCCAAACGTTCCTGCACAATGGCGATGAAATCGCGTAATG  >  NC_000913/2667589‑2668061

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