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.96.4% 249.9 / 0.8 56pseudogene (183/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (52/2);  total (52/4)
Fisher's exact test for biased strand distribution p-value = 3.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0077794,371,2720A.96.4% 251.3 / ‑0.6 56pseudogene (182/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (52/2);  total (52/4)
Fisher's exact test for biased strand distribution p-value = 3.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAA  >  NC_007779/4371023‑4371521
                                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                         
tGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaa                                                                                                                                                                                                                                                           >  1:640991/1‑250 (MQ=255)
   ggTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaata                                                                                                                                                                                                                                                        >  1:266476/1‑248 (MQ=255)
     tGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagc                                                                                                                                                                                                                                                      >  1:491979/1‑246 (MQ=255)
      gATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagcg                                                                                                                                                                                                                                                     >  1:537566/1‑245 (MQ=255)
             aTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATaa                                                                                                                                                                                                                                            <  1:535575/250‑1 (MQ=255)
                 gcgcGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGc                                                                                                                                                                                                                                        >  1:454505/1‑250 (MQ=255)
                 gcgcGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGc                                                                                                                                                                                                                                        <  1:235361/250‑1 (MQ=255)
                        tACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATa                                                                                                                                                                                                                                 >  1:11721/1‑250 (MQ=255)
                         aCCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATaa                                                                                                                                                                                                                                >  1:931348/1‑250 (MQ=255)
                                cAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGAt                                                                                                                                                                                                                         >  2:476666/1‑250 (MQ=255)
                                          gggAGGGAAAGTTTTACCTCTTTTTTCCTCAGCAACTTCAGCAAAACTTGGAGGTGCCGAAAGAGAATTTTGCCAAGTACAAGAAAAAATGGACAGAAAATGGGCCAGGATGGGGGGGGGGAGGGTGATTTATTCCGGTTACGGGCGCAGCACTTTTTCTGTCTAATAACCCAGATAGCCCACCCCGCCAGCACCCGGCAGGGAAAAAATGGGGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                                                                                                                                                                                                               <  2:488408/249‑1 (MQ=255)
                                          cGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAACGGCGATAACCGCCATGATAATGGAGATGACATCGAAc                                                                                                                                                                                                               >  2:209068/1‑250 (MQ=255)
                                                   aGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGTACCAGATAGTCCAGACCGCCAGCAACCTGCATGACAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGtga                                                                                                                                                                                                      >  1:842271/1‑248 (MQ=255)
                                                   aGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTa                                                                                                                                                                                                      >  1:754371/1‑250 (MQ=255)
                                                         aaCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGt                                                                                                                                                                                                >  1:43045/1‑250 (MQ=255)
                                                            gCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTa                                                                                                                                                                                             >  1:192229/1‑250 (MQ=255)
                                                              cTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAAc                                                                                                                                                                                           >  2:353727/1‑250 (MQ=255)
                                                                   tCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCaa                                                                                                                                                                                      >  1:783490/1‑250 (MQ=255)
                                                                       tCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGc                                                                                                                                                                                  >  2:303863/1‑250 (MQ=255)
                                                                        cGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCg                                                                                                                                                                                 >  1:46110/1‑250 (MQ=255)
                                                                           aaCTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCa                                                                                                                                                                              >  2:736347/1‑250 (MQ=255)
                                                                            aCTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCaa                                                                                                                                                                             >  1:24307/1‑250 (MQ=255)
                                                                               tCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGaa                                                                                                                                                                          >  2:729099/1‑250 (MQ=255)
                                                                                        aCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCAccc                                                                                                                                                                 >  1:696992/1‑250 (MQ=255)
                                                                                               gtgtCGGAAGAGAGATTTTGCGATTACCAAAAAAAAAAGTAAGAAAATGGTCCGGGAGGTGTGCGAGGATGTTGTTTTTTTTGGTGTTGCGGCGCAGCTGTTTTTGGGTGTGAAGCACCAATACCCCCAGACCCCCAGCAACCTGAAGGGC‑‑AAAAAGAGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATcc                                                                                                                                                          <  2:666137/250‑1 (MQ=255)
                                                                                                 gtCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGc                                                                                                                                                        >  1:695560/1‑250 (MQ=255)
                                                                                                      cAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCa                                                                                                                                                   >  2:594119/1‑250 (MQ=255)
                                                                                                       agagagATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCaa                                                                                                                                                  >  2:144120/1‑250 (MQ=255)
                                                                                                           agaTGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAAcc                                                                                                                                              >  1:111111/1‑250 (MQ=255)
                                                                                                                ttGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCACAACCAATAc                                                                                                                                         >  2:231967/1‑250 (MQ=255)
                                                                                                                           cAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTccccc                                                                                                                               >  2:404275/1‑250 (MQ=255)
                                                                                                                                aaGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATc                                                                                                                          >  1:872421/1‑250 (MQ=255)
                                                                                                                                       tCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGcc                                                                                                                   >  1:774699/1‑250 (MQ=255)
                                                                                                                                            aaaTAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaaga                                                                                                              >  2:767241/1‑250 (MQ=255)
                                                                                                                                                 ggTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGACGATc                                                                                                         >  1:482174/1‑250 (MQ=255)
                                                                                                                                                     aCGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCa                                                                                                     >  1:89233/1‑250 (MQ=255)
                                                                                                                                                        aTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCa                                                                                                  >  1:513401/1‑250 (MQ=255)
                                                                                                                                                                      cGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGt                                                                                    >  2:262021/1‑250 (MQ=255)
                                                                                                                                                                         gATGCATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTCTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCt                                                                                 >  2:164637/1‑250 (MQ=255)
                                                                                                                                                                            gTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACa                                                                              >  1:189018/1‑250 (MQ=255)
                                                                                                                                                                                  cGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGc                                                                        >  1:22756/1‑250 (MQ=255)
                                                                                                                                                                                       tACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATAtt                                                                   >  1:613957/1‑250 (MQ=255)
                                                                                                                                                                                                agcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGccttcc                                                           >  1:421994/1‑250 (MQ=255)
                                                                                                                                                                                                  cagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGccttcctt                                                         >  2:125130/1‑250 (MQ=255)
                                                                                                                                                                                                    gcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGt                                                       >  2:68155/1‑250 (MQ=255)
                                                                                                                                                                                                          ctgtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                 >  1:386831/1‑250 (MQ=255)
                                                                                                                                                                                                                      ccAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGatat                                     >  1:227649/1‑250 (MQ=255)
                                                                                                                                                                                                                          aTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGa                                 >  1:764001/1‑250 (MQ=255)
                                                                                                                                                                                                                              gTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTg                              >  2:183001/1‑250 (MQ=255)
                                                                                                                                                                                                                               tCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGt                             >  1:825034/1‑250 (MQ=255)
                                                                                                                                                                                                                                ccAGACCGCCAGCAACCTGCATGGCGG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt                            >  2:560023/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTa                           >  2:592137/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttt                        >  2:421967/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aCCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttt                        >  2:245229/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTa                     >  2:799227/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAg                   >  2:540367/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTa                >  2:593181/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTa                >  2:827189/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACt              >  2:310299/1‑250 (MQ=255)
                                                                                                                                                                                                                                                cTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGc            >  1:688569/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     aaggaaaaaaaTAGGGGGAAAAGCCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTca         <  2:727722/244‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaaaTGGGGGCGAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGaaa  <  1:209068/250‑1 (MQ=255)
                                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                         
TGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAA  >  NC_007779/4371023‑4371521

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