Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 990,324 C→T 17.9% S152N (AGC→AAC)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779990,3240CT17.9% 145.7 / 9.9 75S152N (AGC→AAC) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (19/40);  new base T (5/8);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCA  >  NC_007779/990088‑990562
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ggCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGCGAAGTCGACGAATTTTCtt                                                                                                                                                                                                                                                  >  1:1240822/1‑239 (MQ=255)
   tGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTcc                                                                                                                                                                                                                                               <  2:1161684/239‑1 (MQ=255)
      tttGCTGATGTGCGTGGAACCGTTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGGCGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATGCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGATTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGc                                                                                                                                                                                                                                           <  2:555733/239‑1 (MQ=255)
       ttGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGATAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCg                                                                                                                                                                                                                                          <  2:546321/239‑1 (MQ=255)
            gATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCg                                                                                                                                                                                                                                     <  2:1252126/239‑1 (MQ=255)
                 gCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGc                                                                                                                                                                                                                                <  1:87023/239‑1 (MQ=255)
                  ccTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAAAGCCGAGAAAACGAATAAATTTGTTTTCCACGTGTCGgggggt                                                                                                                                                                                                                               >  1:1351703/1‑237 (MQ=255)
                           attctggtgCCTGGGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACg                                                                                                                                                                                                                    <  2:825836/233‑1 (MQ=255)
                           aTTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGGGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAACGGCGGGAAGACGTCTAATTTTTTTTCCAGCGGGTCGGCGGGGGGCGCAACg                                                                                                                                                                                                                    >  1:441559/1‑239 (MQ=255)
                            ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCCGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTGAACACCCGGAAATAGACCAATATTCTTTCCCTAGTGTGGAGGGCGTTCGCaaac                                                                                                                                                                                                                    >  2:304979/1‑236 (MQ=255)
                            ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCCCGGAGGCGAATAATTTGGTTTCCAGCGTTTTGGGGGGGTGGGCAACGt                                                                                                                                                                                                                   >  1:2098918/1‑239 (MQ=255)
                            ttCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCCAGCAGAGCGCCAAATTGCTTTCCCGCGCGTTGGGGGGCTGCGCaaag                                                                                                                                                                                                                    >  1:1089889/1‑236 (MQ=255)
                                 gtgCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGc                                                                                                                                                                                                              <  1:952121/238‑1 (MQ=255)
                                 gtgCCTGTGTTCCCATCGGCGTGAGGGAAAGACGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGc                                                                                                                                                                                                              <  2:2030530/238‑1 (MQ=255)
                                      tgtgTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCGAGAAGTCGGCGAATTTGTTTTTCACCGTGTCGGGGGCTTGCGCGACGTCCGCCTCCgg                                                                                                                                                                                                        >  1:1222456/1‑239 (MQ=255)
                                         gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                      <  1:1059553/238‑1 (MQ=255)
                                         gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                      <  1:1387386/238‑1 (MQ=255)
                                           tCCCCTCGGCGTGAAGGAAAGCGGGCGCGCCGGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGTCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGag                                                                                                                                                                                                    <  2:1227458/238‑1 (MQ=255)
                                           tCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGag                                                                                                                                                                                                    <  1:948646/238‑1 (MQ=255)
                                               aTCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCGCGAAAGCGACTAATTTTTTTTTCAGCGTGGCGGGGGCGTGCGCCACGTCGGCCTGCGGGGagggata                                                                                                                                                                                               >  1:514735/1‑239 (MQ=255)
                                                tCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatag                                                                                                                                                                                              <  1:967367/239‑1 (MQ=255)
                                                  ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCTAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  1:1867622/239‑1 (MQ=255)
                                                  ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  2:733671/239‑1 (MQ=255)
                                                  ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  1:156272/239‑1 (MQ=255)
                                                  ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCGGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  1:126229/239‑1 (MQ=255)
                                                  ggCGTGAAGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  1:1120748/239‑1 (MQ=255)
                                                   gCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGGCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTTTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                      <  1:1916812/228‑1 (MQ=255)
                                                   gCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTGTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                            <  2:1232274/238‑1 (MQ=255)
                                                    cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGCGAAGGAGACTCATTTGCTTTTCAGACGGGTGGGGGGGTGCGCCACGGCGGGCTGCGGGGggcggtagagat                                                                                                                                                                                          >  1:723279/1‑239 (MQ=255)
                                                     gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACACACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                         <  1:966238/239‑1 (MQ=255)
                                                     gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                         <  1:1048705/239‑1 (MQ=255)
                                                     gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                         <  1:1392710/239‑1 (MQ=255)
                                                     gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                         <  1:127355/239‑1 (MQ=255)
                                                       gAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTa                                                                                                                                                                                       <  1:168327/239‑1 (MQ=255)
                                                                    cgcgCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCCACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGAt                                                                                                                                                                          >  1:1992816/1‑239 (MQ=255)
                                                                         cAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTa                                                                                                                                                                                       <  1:1743788/221‑1 (MQ=255)
                                                                          aGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCCAAAATATGCATAATTTTCCTTTCAGCGGGTCGAGGGCCGTCGCAACCTCGGCCTTGGGGTAGCGTGAGCGATGTACCATTTCACTGACaaact                                                                                                                                                                     >  2:195094/1‑235 (MQ=255)
                                                                                               acaAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGAGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCACAGAAGTCGACTAAATTGCTTTCCAGCGTGTCGAGGGCGTGCGCAACGTCGGCCTGCGGGGGGCGATAGCGATGTACCACTTCTCTGATAAACGCCCGCAAAGGAACCTCCCagag                                                                                                                                               >  1:160297/1‑239 (MQ=255)
                                                                                                 aaaCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAAAGCCGAGAAGTGGACTAATTTGGTTTTCAGCGTGGCGAGGGGTTGCGCGAACGTGGCCCGCGGTGAGGGAGAGAGCAGTACCACTTCTCTGAAAAACCCCAGCCAAGGAACCTTCCCGAGGa                                                                                                                                             >  2:1387799/1‑239 (MQ=255)
                                                                                                         aTTCCTGTGGCGGCCGCTAAACGGTGGTTTCTTGGCCTCCGGGGAAAAACCGGCGCGGGGGGCAAAAACCCGTCGGATGGAGGCGCAAAAAATAAAGCCCGGCGGTTACCGCCGAGAAGGCGATTAATTGGCTTTCCGGGGGGTGGGGGGCTTGCGCAGCGTCGGCCTGGGGTGGGCGATAGAGGTGTCCCTTTCCACGGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                      <  2:1222456/238‑1 (MQ=255)
                                                                                                             cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTTCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAAACGCCCGCAAAAGaa                                                                                                                                                         >  1:1317638/1‑215 (MQ=255)
                                                                                                               gTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCCAGAAGTCGACTAATTTGCTTTTCAGCGGGTCGAGGGCTTTCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTGATAACCGCCAGCCAAAGAACTTCCCCGAGGATGACTTTCCGccc                                                                                                                                >  1:1353043/1‑237 (MQ=255)
                                                                                                                        gCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCtt                                                                                                                      >  1:88441/1‑239 (MQ=255)
                                                                                                                            aaCGATGGTTTCTTGTGCTCCGGGAAAAAACGGGGGCGGGGGCCAAAAACCAGCGGAAGGGAGGGGCGCATGTCAGACACGGGCGTTTACCGCGGAGAGGGCGGCTAATTGGTTTCCCGGTGTGGGGGGGGTTGCGGCAGCGTCGGCCGGCGGGGAGGGAGAGCGTGCTCCGATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAAc                                                                                                                  <  2:277131/239‑1 (MQ=255)
                                                                                                                                 tGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGGCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                         >  1:1764203/1‑195 (MQ=255)
                                                                                                                                   gTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAGGACCGGGGGTTAACGCCGAGAAGTCGCCTAATTTGCTTTCCAGCGTGTGGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCTTTTCACTGATAACCGCCAGCAAAgg                                                                                                                                                           <  1:385115/191‑1 (MQ=255)
                                                                                                                                    tttCTTGTACTTCGCGAGAAAAACGGCGACGGGTGTCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGCGACCAATTTGATTTTCAGCGTGTCGGGGGGTTGCGCAACCGCGGCCTGCGGGGGGGGAGAGCGATGTACCATTTTACTGATAACCGCCCGCAAAAGGACATTCCCGAGGAGGGAATCTCCCCACGGGCCGCCTTAACAAatcgcc                                                                                                          >  2:543117/1‑235 (MQ=255)
                                                                                                                                      tCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCCAGAAGTCGACTCATTTGTTTTCCAGCGTGTCGAGGGGATGCCCAACATCGCCCCCCGGCGGGCGATCGCGATGGATCGTTTCTCTTATAACCGCCCGCCAACGAAGTTCCCAGAGGGTGAATTCACACCCCTCACCGCCTTAACCAATATcac                                                                                                         >  2:1777142/1‑236 (MQ=255)
                                                                                                                                             cTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCTTCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCCAGCAATCGACTAAATTTCTTTCCACCGTGTCGAGGGGCTGCGCCACGTCGGCCTGCGGGGAGCGCTAGCGATGTACCATTTCTCTCATAACCGCCAGCAAAGGGACATCCCCGAGGGTGACTTTACGCCCCAGGGCCCGTAAACGAAAATcccacctgcc                                                                                                  >  2:554552/1‑230 (MQ=255)
                                                                                                                                              ttCGCGAGAAAACCGGCGACGGGTGCAAAAATCCTTCAGAGGGAAGCGCGAAATACCAAGCCCGTCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGTTTGCGAAACGTCGGCCTGCGGGGAGAGATAGCGATGTACCTTTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                                                                <  2:1353043/239‑1 (MQ=255)
                                                                                                                                                cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCCAGAAGTCGAATAATTTTCTTTTCCGCGCGGCGAGGGGCTGCGCCACGTCGGCCTGCGGGGGGGGGTAGCGGTGGTACATTTCACTGGTAACCGCCCGCAAAGGGACTTCCCCGAGGATGACTTCACCCCACGGGCCGCCTAAACAAAAATCCACCCTGCCATTa                                                                                              >  2:156904/1‑239 (MQ=255)
                                                                                                                                                cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCCAGAAATCGAATAATTTTCTTTCCCACCTGTCGAGGGGGTGCCCAACGTCGGCCTGCCGGGAGCGATAGCGATGTAACATTTCACTGATAACCGCCCGCAAAGGAACTCCCCCGAGGGTGACTTCACGCCCCGGGCCCCCTACACGAAAATACAACCTGcgcttt                                                                                              >  2:731884/1‑234 (MQ=255)
                                                                                                                                                     gaAAAACGGCGAGGGGTGCCAAAGCCCATCAAAGGGAAGCGCTGAACAACAAGCCGGGGGTTTCTCGCCGGGAAGCCGCATATTTTGTTTGCCGGAGGGAGGGGGGTTTGCCACTTTCCGGCCGGGGGGGGGCAATGCGGGTGAACCTTTTCATTGAAATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCg                                                                                           <  2:1539586/237‑1 (MQ=255)
                                                                                                                                                                 cGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGCCCGGCGGTTACCGCCGAGAAGTCGACTATTTTGCTTCCCAGCGTGTCGAGGGCTTGCGCACCGTCGGCCGGCGGTGAGCGATAGCGATGTCCCATTTCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAc                                                                             <  1:156904/239‑1 (MQ=255)
                                                                                                                                                                       ccAAAATCCATGAGGTGGAACGGCGACATATCAGGCCCGGCGGTTCAGGCCGAGAAGTCGTCTAATTTGCTTTCGCGCGTGTCGAGGCCTGGCGCAACGGCGGCCTGGGGGGAGCGATAGAGATGTACGATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTc                                                                       <  2:1850461/239‑1 (MQ=255)
                                                                                                                                                                             tCCCTCAGAGGGAAGCGCGAAATATCAAGACCGGCGGTTAACGCGGAGAAGTCGCCTAATTTGTTTTCCAGCGTGTGGAGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGtt                                                                 <  1:731884/239‑1 (MQ=255)
                                                                                                                                                                                          ggcgcGAATTATCAAGCCCGGCGGTTAACGCGGAGAGGTCGACTAATTTGTTTTCCAGCGTGTCGAGGGTTGGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCTTTTCACGGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                      <  2:155283/156‑1 (MQ=255)
                                                                                                                                                                                               gataaacgaaGACCGGCGTTTACGGCGGAAGAGCCCATTATTTTGTTTCCCGGCGTGTCGAGGGTTGGCGAAACGTCGGCCGGCGGTAGGCTATGGCGATGACCCATTTCACAGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                               <  2:723279/231‑1 (MQ=255)
                                                                                                                                                                                               gaaatatCAAGCCCGGGGGTGAACGCCAAGAGGTCGACTATTTTGTTTCGCGGCGTGTGGAGGGCTTGCGCAACGTCGGCCGGCGGGGAGCGATAGCGATGTCCCATTCCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                               <  2:1587940/236‑1 (MQ=255)
                                                                                                                                                                                               gACATACCAGGACGGGCGGTTACGGCGGAGAAGCCGATTAATTTGCTTTCCAGGGTGTCGAGGGTTTGCGCAACGTCGGCCTGCGGTGGGGGATGGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                               <  1:543117/239‑1 (MQ=255)
                                                                                                                                                                                                 ataatCGAGACCGGCGGGTAACGGCGAGAGTGCGACTATTTGGTTTTCCAGCGTGTGGAGTTCTTGCGAGACGTCGGCCTGGGGTGAGAGATAGCGATGCACCATTTCTCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAAcc                                             <  2:160297/236‑1 (MQ=255)
                                                                                                                                                                                                    atCAAGACCGGCGGTTAACGCCGAGAAGTCGGCTAATTTTCTTTCCAGCGTGTCGAGGGCTTGCGCCACGTCGGCCTGCGGGGAGCGATAGCGATGTAACCTTTCACTGATAACCCCCAGCCAAAGAACTTCCCCGAGGGTGACTTTACGcccccgg                                                                                                                            >  1:155283/1‑152 (MQ=255)
                                                                                                                                                                                                    atCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGTTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTa                                                                                              <  1:1869731/187‑1 (MQ=255)
                                                                                                                                                                                                              ggcggTAAACGCCGAGAAGTCGACTATTGTGCTTTCCGGGGTGTCGCGGGCTTGCGCAACGTGGCCCGGCGGGGAGCGATAGCGATGTACCATTTCACTGAAACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                 <  1:1656320/239‑1 (MQ=255)
                                                                                                                                                                                                               ggggtaaccGCCGAGAGGGCGCTTATTTGGCTTCCCAGCGTGTCGAGGGCTTGCGAAACGTCGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                 <  1:506449/230‑1 (MQ=255)
                                                                                                                                                                                                               gcgtTTAGCGCCGAGATGTGACCTATTTGCGTTTCCAGGGTGTCGGGGGCTTGCGCAACGTCGGCATGGGGGGGGCGGTAGCGATGTACCATTTCCCTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                 <  2:1240822/238‑1 (MQ=255)
                                                                                                                                                                                                               gcggTTAACGCGGAGAAGTCGACTATTTTGCTTCCCGGGGTGTCGAGGGTTTGCGCAACGTGGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                 <  1:1759609/238‑1 (MQ=255)
                                                                                                                                                                                                                    gaccGCCGAGATCTGGACTAATTGGCTTCCCAGCGTGTCGGGGGCTTGCCCAACGTGGGCCTGCGGAGAGCAATAGCGATGCACCTTTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                           <  1:1732736/236‑1 (MQ=255)
                                                                                                                                                                                                                    gaaCGCCGAAAAGTGGACTATTTTTCTTTCCAGCGTGTGGAGGGTTTGCGCAACGTCGGCCGGCGGTGAGCGATAGCGGTGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCt                           <  1:304979/238‑1 (MQ=255)
                                                                                                                                                                                                                          cGAGATGTCGTGTATTTTGCTCCCCAGCGTGTGGCGGGCTTGCGAAACGCCGGCCGGGGGTGAGCATTAGCGATGTACCTTTTCATTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGGACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaag                      <  2:514735/238‑1 (MQ=255)
                                                                                                                                                                                                                            agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                   >  1:37923/1‑239 (MQ=255)
                                                                                                                                                                                                                                 tCGATTAATTTGTTTTCCGCGGTGTCGAGGGTTGGCGAAACCTCGGCCTGCGGTGAGGAATTGCGATGTCCCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggca               <  1:1387799/238‑1 (MQ=255)
                                                                                                                                                                                                                                 tCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag              >  1:302642/1‑239 (MQ=255)
                                                                                                                                                                                                                                  cGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgccagg             >  2:2055140/1‑239 (MQ=255)
                                                                                                                                                                                                                                            tGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAg     >  1:62272/1‑238 (MQ=255)
                                                                                                                                                                                                                                              cTTTCCAGCGTTTGGAGGGCTTGCGCAACGGCTGCCTGGGGTGGGCGATAGCGATGTCCCTTTTCCGTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa  <  1:1620775/239‑1 (MQ=255)
                                                                                                                                                                                                                                              cTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc   >  1:219619/1‑238 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
GGCTGATTTGCTGATGTGCCTGGAACCATTCGTGTGCCTGTGTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCA  >  NC_007779/990088‑990562

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