Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,002,969 A→G 29.1% intergenic (‑196/‑36) pitA ← / → Y75_RS19800 inorganic phosphate transporter PitA/NAD(P)/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,002,9690AG29.1% 64.9 / 28.8 41intergenic (‑196/‑36)pitA/Y75_RS19800inorganic phosphate transporter PitA/NAD(P)/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (8/21);  new base G (5/7);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCCCGGTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTACCAACCTTTATGTCGAACCAGGCGCTTATCTGAGC  >  NC_007779/4002731‑4003196
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
aCAGCCCGGTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTa                                                                                                                                                                                                                                     <  2:831374/239‑1 (MQ=255)
   gCCCGGTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCaa                                                                                                                                                                                                                                  >  2:1818836/1‑239 (MQ=255)
    cccGGTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAAt                                                                                                                                                                                                                                 >  2:1707141/1‑239 (MQ=255)
        gTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCGTTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCAc                                                                                                                                                                                                                             >  1:1681206/1‑239 (MQ=255)
         tATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCAc                                                                                                                                                                                                                             >  2:339441/1‑238 (MQ=255)
                aTCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttt                                                                                                                                                                                                                     <  2:1810375/239‑1 (MQ=255)
                  ccAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttt                                                                                                                                                                                                                   <  1:1257859/239‑1 (MQ=255)
                                 aaTGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCa                                                                                                                                                                                                    <  2:604849/239‑1 (MQ=255)
                                         aTAGGACTTTACGGCGTTTTGAGGACATGAACGCGGCGCCTTATCAGTTACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGTATTTCACTCAAAAAATAATTATCATATATTATTCATGCAAATACTTTCTCTTAGTAATATTGATGCTGCTGCTACCACTGAGGAATCTTCGCAATTCACGCCCGTTTCTTCTAAGAGGGGCGCAACCTGGa                                                                                                                                                                                             <  2:503823/238‑1 (MQ=255)
                                              aCGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGtt                                                                                                                                                                                       >  2:78252/1‑239 (MQ=255)
                                                gTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTg                                                                                                                                                                                     <  1:1818836/239‑1 (MQ=255)
                                                  tACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGa                                                                                                                                                                                    <  1:1566871/238‑1 (MQ=255)
                                                  tACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGAt                                                                                                                                                                                   >  2:1527796/1‑239 (MQ=255)
                                                  tACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGAt                                                                                                                                                                                   <  1:782353/239‑1 (MQ=255)
                                                         ttttGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGc                                                                                                                                                                                 <  2:405771/234‑1 (MQ=255)
                                                         ttttGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGc                                                                                                                                                                                 >  1:405771/1‑234 (MQ=255)
                                                                   aggAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATACTTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGgcgt                                                                                                                                                                  >  2:667932/3‑238 (MQ=255)
                                                                   aTGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGgc                                                                                                                                                                    >  1:163189/1‑237 (MQ=255)
                                                                           ggcgCATTATCAGTGACTTTCGCGTCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGACTTCACTCAAAAAATAATTATCTGATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTgcgg                                                                                                                                                          <  2:1471945/239‑1 (MQ=255)
                                                                           ggcgCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTgcgg                                                                                                                                                          <  2:1510336/239‑1 (MQ=255)
                                                                             cgcATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTgcggcg                                                                                                                                                        >  1:1525854/1‑239 (MQ=255)
                                                                              gcATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGcggcgg                                                                                                                                                       <  2:567031/239‑1 (MQ=255)
                                                                                aTTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGt                                                                                                                                                     <  1:187087/239‑1 (MQ=255)
                                                                                     cAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTAtgtt                                                                                                                                                <  1:1707141/239‑1 (MQ=255)
                                                                                          aCTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCACTGGTGCGGCggggttggttttgt                                                                                                                                           >  1:362383/1‑228 (MQ=255)
                                                                                                    cgcggGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGaa                                                                                                                                                                                            <  1:1858057/180‑1 (MQ=255)
                                                                                                    cgcggGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGaa                                                                                                                                                                                            >  2:1858057/1‑180 (MQ=255)
                                                                                                                   atatAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACg                                                                                                                  <  2:1525854/239‑1 (MQ=255)
                                                                                                                    tataGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAgg                                                                                                                              <  1:904154/226‑1 (MQ=255)
                                                                                                                                         ccccTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAAt                                                                                            <  2:1803201/239‑1 (MQ=255)
                                                                                                                                                                 aataatTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCaaa                                                                       <  1:1047549/236‑1 (MQ=255)
                                                                                                                                                                                     cAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTgcgg                                                                                                                                                          <  2:1848082/133‑1 (MQ=255)
                                                                                                                                                                                     cAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTgcgg                                                                                                                                                          >  1:1848082/1‑133 (MQ=255)
                                                                                                                                                                                             tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCaa                                          <  1:210229/237‑1 (MQ=255)
                                                                                                                                                                                                  ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTAcc                                   <  1:596055/239‑1 (MQ=255)
                                                                                                                                                                                                  ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTAcc                                   <  1:1856086/239‑1 (MQ=255)
                                                                                                                                                                                                  ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTAcc                                   <  1:1293335/239‑1 (MQ=255)
                                                                                                                                                                                                             aagtGATGCTGGTGCGACCACTGAGGAATTTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTACCAACCTTTATGt                        <  2:362383/236‑1 (MQ=255)
                                                                                                                                                                                                                      tGGTGCGACCACTGAGGAATCTTTGCAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTACCAACCTTTATGTCGAACCAgg               <  1:339441/239‑1 (MQ=255)
                                                                                                                                                                                                                        gTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTc                                                                                                                            <  1:1365317/128‑1 (MQ=255)
                                                                                                                                                                                                                                   gAGGAATCTTTACACTTCACGGCCGTTTTTTCTAAGAGGAGCGCAACCTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTCTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTACCAACCTGTATGTCGAACCAGGCGCTTATCTGAGc  <  2:965903/239‑1 (MQ=255)
                                                                                                                                                                                                                                           tttACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTctg                                                                                                        <  1:166132/129‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                   
ACAGCCCGGTATGCAAATCCAGGCCAGCAAACAAATGTAGCATAGGACGTTACGCCATTTTGAGGACATGAACGCGGCGCATTATCAGTGACTTTCGCGGCGCGGGCAAAGTGAAATATAGACTTTTTTTGATTTGCCCCCTGTATGGATTTCACTCAAAAAATAATTATCTTATATAATTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAGGTTTGATGCCATTATTATAGGCGCTGGTGCGGCGGGTATGTTCTGTTCTGCGCTGGCAGGTCAGGCAGGACGCCGGGTTCTGCTGATCGATAATGGTAAAAAACCAGGGCGCAAAATCCTTATGTCTGGCGGTGGGCGCTGCAACTTTACCAACCTTTATGTCGAACCAGGCGCTTATCTGAGC  >  NC_007779/4002731‑4003196

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