Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007779 4,371,274 (A)5→3 pseudogene (179‑180/1304 nt) Y75_RS21600 ← anaerobic C4‑dicarboxylate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,371,2710A.97.1% 154.5 / ‑1.3 35pseudogene (183/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (33/1);  total (33/2)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_0077794,371,2720A.97.1% 156.4 / ‑3.3 35pseudogene (182/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (33/1);  total (33/2)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTAC  >  NC_007779/4371022‑4371508
                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                            
aTGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGa                                                                                                                                                                                                                                               >  1:150849/1‑250 (MQ=255)
 tGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGaa                                                                                                                                                                                                                                              >  1:409092/1‑250 (MQ=255)
  gCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGaaa                                                                                                                                                                                                                                             <  1:918659/250‑1 (MQ=255)
     gTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGaaatag                                                                                                                                                                                                                                          >  1:718946/1‑247 (MQ=255)
       gATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGaaatagcg                                                                                                                                                                                                                                        <  2:611228/250‑6 (MQ=255)
             aaTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTGCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                                <  1:438598/250‑1 (MQ=255)
                                 cAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                                                                                                                                                                                                            >  1:851145/1‑250 (MQ=255)
                                     aGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACa                                                                                                                                                                                                        >  1:386935/1‑250 (MQ=255)
                                         aGCTGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGa                                                                                                                                                                                                    >  1:820451/1‑250 (MQ=255)
                                          gCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGCTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGCATTTCGGGATACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGAAACCTGCATGGCAG‑‑AAATAGCGCCGATAACCGCCATGATAATGAAGCTGACAACGaa                                                                                                                                                                                                   >  2:104087/1‑250 (MQ=255)
                                           cGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                                                                                                                                                                                                  >  2:501517/1‑250 (MQ=255)
                                            ggACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGCGATGACCACGAACg                                                                                                                                                                                                 >  2:585013/1‑250 (MQ=255)
                                               cGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGa                                                                                                                                                                                              >  2:870462/1‑250 (MQ=255)
                                                          aaCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt                                                                                                                                                                                   >  1:225153/1‑250 (MQ=255)
                                                                ttGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAATCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACg                                                                                                                                                                             >  2:257376/1‑250 (MQ=255)
                                                                          acAACTTCAGCGAAAACTGGGGGGGTCAACGGAGAGAGGTTGCCAGAACAAGAAAAGATAGTCAGAAAAGGGGTCAGGATGGGGGCGAGGATGGTGTTTTTTGGCGGGTGCGGGCGCAGCAGTTTTTGTGTCTGAGAACCCAGATAGCCCAGC‑CCGCAAGCAACCGGCA‑GGGA‑AAAAATGGGGGAGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCgg                                                                                                                                                                   <  2:589208/249‑1 (MQ=255)
                                                                                 cAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAAc                                                                                                                                                            >  1:860279/1‑250 (MQ=255)
                                                                                           tGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAccccc                                                                                                                                                  >  2:782516/1‑250 (MQ=255)
                                                                                                 tgtCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCg                                                                                                                                            >  2:802950/1‑250 (MQ=255)
                                                                                                 tgtCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCg                                                                                                                                            >  2:1063461/1‑250 (MQ=255)
                                                                                                            agaTGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAAcc                                                                                                                                 >  2:46434/1‑250 (MQ=255)
                                                                                                                    ccAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTa                                                                                                                         >  1:172034/1‑250 (MQ=255)
                                                                                                                        tACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcc                                                                                                                     >  1:600176/1‑250 (MQ=255)
                                                                                                                                aaGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCACAACCAATACCTATTCCCCCCAATc                                                                                                             >  2:423521/1‑250 (MQ=255)
                                                                                                                                aaGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATc                                                                                                             >  1:17413/1‑250 (MQ=255)
                                                                                                                                               tAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGACGCCCAagaaga                                                                                              >  2:351104/1‑250 (MQ=255)
                                                                                                                                                  gTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCg                                                                                           >  1:1028807/1‑250 (MQ=255)
                                                                                                                                                     aCGGTCGGTGGGAGTATCTGGATTTATTTGGTGTCCCGGCGCAGCTGTTTTTTTGTCTGATGACCCAGATCCCCCACC‑CCGCAACCACCCTGCATGGC‑‑AAAAGTGGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCa                                                                                        <  2:600176/250‑1 (MQ=255)
                                                                                                                                                     aCGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCa                                                                                        >  2:725354/1‑250 (MQ=255)
                                                                                                                                                         tCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCaa                                                                                    >  1:147993/1‑250 (MQ=255)
                                                                                                                                                                               tttCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACt                                                              >  1:415445/1‑250 (MQ=255)
                                                                                                                                                                                ttCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTa                                                             >  1:932645/1‑250 (MQ=255)
                                                                                                                                                                                      ttACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCatat                                                       >  1:856411/1‑250 (MQ=255)
                                                                                                                                                                                                  tctttttCTGGCTGAGGAA‑CAAAAAGCCCAAACCCGCCAGAAACCGGCATGG‑A‑AAAAATGGCGGCAATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTg                                           <  2:652090/247‑1 (MQ=255)
                                                                                                                                                                                                     ttttCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGttt                                        >  1:247399/1‑250 (MQ=255)
                                                                                                                                                                                                                 tgaACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTg                            >  1:644497/1‑250 (MQ=255)
                                                                                                                                                                                                                        aaatcgcccggcccGCCAGCAACCTGCAGGGC‑‑AAAAATAGCGGCAATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTa                      <  1:4998/238‑1 (MQ=255)
                                                                                                                                                                                                                             cccaagcccGCCAGCAACCGGCGTGG‑A‑AAAAAGAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTg                 <  1:1022241/243‑1 (MQ=255)
                                                                                                                                                                                                                              ccAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt               >  1:399242/1‑250 (MQ=255)
                                                                                                                                                                                                                              ccAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt               >  2:470466/1‑250 (MQ=255)
                                                                                                                                                                                                                               cAGA‑CCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTa              >  1:734055/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cccGCCAGCACCCTGCAGGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttt           <  1:342752/249‑1 (MQ=255)
                                                                                                                                                                                                                                       ccAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTaa       >  2:1071397/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAg      >  1:273799/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTCCACCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATCTTAGATTGTTATTTTTAAGTTAc  >  2:305106/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                            
ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGA‑CCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTAC  >  NC_007779/4371022‑4371508

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