Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 363,821 T→C 17.4% D570G (GAC→GGC)  Y75_RS01770 ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779363,8210TC17.4% 127.7 / 2.2 47D570G (GAC→GGC) Y75_RS01770beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base T (34/4);  new base C (2/7);  total (36/11)
Fisher's exact test for biased strand distribution p-value = 1.57e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTGGCCGTGGTGTCGGCTC  >  NC_007779/363590‑364058
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
tGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGAc                                                                                                                                                                                                                                        >  1:457928/1‑239 (MQ=255)
      aaCAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAgcc                                                                                                                                                                                                                                  >  1:575089/1‑239 (MQ=255)
          gttgttgcggggTCATTTGGATGGTTGCCACGGAAAAAGGGACGAGAAAAACTGGCGGGTGGTTTTTGGCCTCCGACGGCGCGGGGTGGGGCGGGCGGGGGCAAAAACCAACACCTTTAAAAAAAACCGGGCGTTGGTTGGGCGTCCCCCAAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGccc                                                                                                                                                                                                                              <  2:88354/229‑1 (MQ=255)
                cGCTGCTCATTTGGATGTTTGCCCGGGAAAACGGAAACGGAAAAACCTGTTGTTGGTTTTTTCCTCCCGTAGGGGCGGGAGGGGGGGGGGGGCCGAAAAAACCCAGCCCGTTCATAAAAAACGGGCGATCGTTGGGGTTACCGCAAAATCCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTaaa                                                                                                                                                                                                                        <  2:858669/239‑1 (MQ=255)
                     gTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                                                                                                                                                                                   >  2:358387/1‑239 (MQ=255)
                     gTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACgggg                                                                                                                                                                                                                   >  2:1954076/1‑239 (MQ=255)
                     gTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACCAAGCCGCCCTGTAAACgggg                                                                                                                                                                                                                   >  2:1603109/1‑239 (MQ=255)
                         tttgggtGGTTCGCCGGGATAAACGAACGGAAAAAAGCTGCTGCTGGTGTGTTGCTCCCGGCAGGGGGGGAGGGGGTGTGCGGGCGGAAAACCCCAGACCGTTCAAACAGAACTGGCGTGCGTGCGGCTTACCACAAAAACCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATAc                                                                                                                                                                                                               <  2:948923/233‑1 (MQ=255)
                          ttCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGTGAt                                                                                                                                                                                                                 >  1:301176/1‑236 (MQ=255)
                                  tttGCCCGGATAAAAGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGa                                                                                                                                                                                                                                         >  1:672771/1‑204 (MQ=255)
                                    tGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGCCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGCAACg                                                                                                                                                                                                     >  2:1480990/1‑238 (MQ=255)
                                        cGGATAACGGGACCTGGAAAACCTGCTGCGGGTGTTTTGCTCCCGTCAGCGCGGGAGGCGGCGTGCGGTCGGAAAAGCCCAGCCCGTTCATCCAGACCTGGCGATCGTTCGGCGTATCGCAAAAATCCCCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGc                                                                                                                                                                                                <  1:1182160/239‑1 (MQ=255)
                                                             aCTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCa                                                                                                                                                                           >  1:119404/1‑239 (MQ=255)
                                                                  gctgctgGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACg                                                                                                                                                                                                     >  1:555027/1‑210 (MQ=255)
                                                                                                     ggcgTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGa                                                                                                                                      >  1:1921795/1‑238 (MQ=255)
                                                                                                        gTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTCTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTAcc                                                                                                                                                                >  1:1435590/1‑209 (MQ=255)
                                                                                                         tGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGc                                                                                                                                 >  1:1240385/1‑239 (MQ=255)
                                                                                                                            gaccGTTCATACAGAACTGGCGATCGTGCGGCGTATCGCCAAAATCACCGCAGTAACCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGCCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGTGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGt                                                                                                               >  1:266796/1‑239 (MQ=255)
                                                                                                                             accGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGa                                                                                                                                                                        >  1:525110/1‑181 (MQ=255)
                                                                                                                             accGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAggga                                                                                                              >  1:908021/1‑237 (MQ=255)
                                                                                                                                  tCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCa                                                                                                                                   >  1:581465/1‑213 (MQ=255)
                                                                                                                                  tCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGGGGGCGTATTCGCAAAGGATCAGCGGGCGCGCCTCTCCCGGCAGCGac                                                                                                         >  1:1035795/1‑238 (MQ=255)
                                                                                                                                       cAGAACTGGCGATCGTTCGGCGCATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCa                                                                                                    >  2:945418/1‑239 (MQ=255)
                                                                                                                                          aaCTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACggg                                                                                                                                                                                                                    >  1:1249006/1‑124 (MQ=255)
                                                                                                                                            cTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCGTATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCAttttt                                                                                               >  1:433109/1‑239 (MQ=255)
                                                                                                                                               gCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGa                                                                                            >  2:296806/1‑239 (MQ=255)
                                                                                                                                                      ttCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAgg                                                                                                                >  1:11465/1‑212 (MQ=255)
                                                                                                                                                         ggCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCAtt                                                                                   >  1:2058008/1‑238 (MQ=255)
                                                                                                                                                         ggCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGAAACTGCCGAACCGCCTCCCAGAATTAAGCGAAACCCCCAGGCCTGTTACCCATCGCGGGGCCGTATTCGCAAAGGACCAGCGGGCGCGTCTCCCCCGGTAGCGAAAGCCATTTTTTGATGGACCAttg                                                                                  >  2:2022717/1‑238 (MQ=255)
                                                                                                                                                         ggCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAACCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCAtt                                                                                   >  1:1664608/1‑238 (MQ=255)
                                                                                                                                                             tgcatataaaaCCCCCGCCCTCAGCCGGCCGGGGTTGGCCGTTTCCATCTTATTAAACCAGCGACTACCCCCCCAAGCCCCATAGAACGCCCCCCTGAAAGCGGGGATCCTGACAAAGCCCCGGCCAGTATTAAGCAAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGc                                                                              <  2:119404/232‑1 (MQ=255)
                                                                                                                                                                     aaaTCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCgg                                                                      >  1:2041688/1‑239 (MQ=255)
                                                                                                                                                                                        aCCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCg                                                                                >  1:21869/1‑210 (MQ=255)
                                                                                                                                                                                         ccACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATcc                                                   >  2:2208032/1‑238 (MQ=255)
                                                                                                                                                                                                 tGCCGTTTTCACCATTTTTAATCAGCGCCTGATCCCCCCTGCCCAGGCCGAGGCCGCCCTGAAACCGGGGATACGGACGAACCCCCTGCCATTTTTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCAcgcgcgcg                                          <  2:301176/239‑1 (MQ=255)
                                                                                                                                                                                                          catcttATTTAACCGGCGACTGCCCCCCCCAGCCCCAGACGAGCCCCCCCTGTAACGGGGGATACTGACAAACCCCCTGCCGGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCggg                                 <  2:1240385/239‑1 (MQ=255)
                                                                                                                                                                                                          aatcatATTTAACCGGGGACTGACCCCCCCAGCCCCAGACGAAGCCCCCCTGTAAACGGGGTGACTGACGAAACCCCTGCCAGTATTTAGCAAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCggg                                 <  1:2067181/238‑1 (MQ=255)
                                                                                                                                                                                                            tcatATTTAATCAGCGACTGACCCACCCAGCCCCAGACGAACCCGCCCGGTAAACGGGGATCCTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCa                               <  1:1940849/239‑1 (MQ=255)
                                                                                                                                                                                                            caatttttAATCAGCGCCTTCCCCCCCCGGCCCCAGCCGAAGCCGCCCTGTAAACGGGGATCCTGCCGAACCCCCTGCCAGTATTTGGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCa                               <  2:1664608/234‑1 (MQ=255)
                                                                                                                                                                                                                     accaccGACGGACCCCCCCAGCCCCGGAGGAAGCCGCCCTGAAAACGGGGATCCTGCCAAACCGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATc                      <  2:880462/234‑1 (MQ=255)
                                                                                                                                                                                                                     aTCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATc                      >  2:338314/1‑239 (MQ=255)
                                                                                                                                                                                                                          cGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTgg                 >  2:704394/1‑239 (MQ=255)
                                                                                                                                                                                                                          cGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTg                  >  1:476446/1‑238 (MQ=255)
                                                                                                                                                                                                                          cGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTg                  >  2:1408612/1‑238 (MQ=255)
                                                                                                                                                                                                                          cGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAAAGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTg                  >  1:535395/1‑238 (MQ=255)
                                                                                                                                                                                                                            aCTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGCCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTGGc                >  1:882530/1‑238 (MQ=255)
                                                                                                                                                                                                                                       cAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTGGCCGTGGTGTCg      >  1:1974910/1‑237 (MQ=255)
                                                                                                                                                                                                                                          tCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCAt                                                     >  1:1485606/1‑187 (MQ=255)
                                                                                                                                                                                                                                          tCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACa                                      >  1:1493136/1‑202 (MQ=255)
                                                                                                                                                                                                                                          tCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCa                               <  1:2165207/209‑1 (MQ=255)
                                                                                                                                                                                                                                          tCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCa                               >  2:2165207/1‑209 (MQ=255)
                                                                                                                                                                                                                                          tCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGCCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTGGCCGTGGTGTCGGCTc  >  1:645112/1‑238 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                             
TGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCGCTGGATGCGGCGTGCGGTCGGCAAAGACCAGACCGTTCATACAGAACTGGCGATCGTTCGGCGTATCGCCAAAATCACCGCCGTAAGCCGACCACGGGTTGCCGTTTTCATCATATTTAATCAGCGACTGATCCACCCAGTCCCAGACGAAGCCGCCCTGTAAACGGGGATACTGACGAAACGCCTGCCAGTATTTAGCGAAACCGCCAAGACTGTTACCCATCGCGTGGGCGTATTCGCAAAGGATCAGCGGGCGCGTCTCTCCAGGTAGCGAAAGCCATTTTTTGATGGACCATTTCGGCACAGCCGGGAAGGGCTGGTCTTCATCCACGCGCGCGTACATCGGGCAAATAATATCGGTGGCCGTGGTGTCGGCTC  >  NC_007779/363590‑364058

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