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.93.9% 202.4 / 7.7 49pseudogene (183/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (43/3);  total (43/6)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0077794,371,2720A.93.9% 206.4 / 3.7 49pseudogene (182/1304 nt)Y75_RS21600anaerobic C4‑dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (43/3);  total (43/6)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTA  >  NC_007779/4371022‑4371502
                                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                      
aTGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGa                                                                                                                                                                                                                                         <  1:358006/250‑1 (MQ=255)
  gCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaa                                                                                                                                                                                                                                       >  1:163302/1‑250 (MQ=255)
       gATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagcg                                                                                                                                                                                                                                  <  1:163745/250‑6 (MQ=255)
            cAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                           >  1:328579/1‑250 (MQ=255)
             aaTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                          <  2:471630/250‑1 (MQ=255)
             aaTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                          >  2:610863/1‑250 (MQ=255)
              aTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATaa                                                                                                                                                                                                                         <  2:615754/250‑1 (MQ=255)
                            cACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATg                                                                                                                                                                                                           >  1:488816/1‑250 (MQ=255)
                              cTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGga                                                                                                                                                                                                         >  1:731510/1‑250 (MQ=255)
                              cTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCAAGGATAATGga                                                                                                                                                                                                         >  2:250100/1‑250 (MQ=255)
                               tGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGgag                                                                                                                                                                                                        >  1:382767/1‑250 (MQ=255)
                                    tAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGAc                                                                                                                                                                                                   >  1:719505/1‑250 (MQ=255)
                                          gggggcggggaagTTTAACGCTTTTTTT‑CCTTCGAAGCTCAGGGGATAAGG‑GGAGGTGTCGCCA‑GAGAGATGTTGCCAGAACCAAAAAAGATAGGAAGAAAAAGGTTCAGGAGGGGGGGGGGGTTGGTGTTTTTTTGGGGTTTGCGGCGCGCCAGTTTTTCTGTTGAAAGAACAAGATAGCCCAGCCCGCAAGCAACGTGAATGGC‑‑AAAAAGAGGGGCAAAAACCGCCATGATAATGGAGATGACATCGa                                                                                                                                                                                              <  2:496589/238‑1 (MQ=255)
                                                   caGGTTTAACCC‑CTTTTTCCCTCCGCAACTAAAG‑AAAAAGT‑GGGAGGGTGGGCAGGAGAGGTGTGGCAAGTACAAAAAAATATAGTAAGAAAAAGGTTCAGGATGGGTGGGAGGATGGTAATGTTTTTGGGTTTCCGGCGCAGAAGTTTTTTTGTTTGAAAACCCAAAAAGCCCAGCCCGCCAGCAACCTGCATGG‑A‑AAAAATAGCGGAGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGt                                                                                                                                                                                     <  2:41903/249‑1 (MQ=255)
                                                       tttAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCg                                                                                                                                                                                >  1:154095/1‑250 (MQ=255)
                                                                  tGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGc                                                                                                                                                                      >  1:1026976/1‑250 (MQ=255)
                                                                   gTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GGGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGcc                                                                                                                                                                     >  2:1000334/1‑250 (MQ=255)
                                                                   gTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATTGAACGGGATGTTACCGGGTTTAACGcc                                                                                                                                                                     >  1:466535/1‑250 (MQ=255)
                                                                        cTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCCGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATACTGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCACTa                                                                                                                                                                 >  1:466401/1‑250 (MQ=255)
                                                                          ccGCAACTGCACCGA‑AACTGGGCGGGGACGCCG‑GAGGGATGTCGCCAGTCCAAACAAAGATAGACAGAAAATGGGTCACGCTGGGCGCGGGTCTCTTTATTTATTTCGTTTTCCGCCCCCGCGTCTTTTCTGTCTGATAACCCATATAGCCCACCCCCCCACCACCCTGCATGGAAG‑‑AAATGGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGc                                                                                                                                                               <  2:103234/249‑1 (MQ=255)
                                                                            gCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAATGGGATGTTACCGGGTTTAACGCCAATAGCgg                                                                                                                                                             >  2:781400/1‑250 (MQ=255)
                                                                            gCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCgg                                                                                                                                                             >  1:996572/1‑250 (MQ=255)
                                                                                    aGCGACAACG‑GGGGGTGCCGCAG‑GAGGGATTTGGCAAGAACCAGAAAAGATAGTGAGAAAATAGGTCACGATGGGGGGGAGGATGGTGATGTATTTGGGGTTACGGCGCAGCATTTTTTCTGTCTGATAAACCAGATAGCCCAGCCCGCCAGCAACCTGCATGG‑A‑AAAAAGGGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAAcc                                                                                                                                                     <  2:382767/250‑1 (MQ=255)
                                                                                                  aGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATc                                                                                                                                        >  2:744243/1‑250 (MQ=255)
                                                                                                     gtCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGc                                                                                                                                     >  1:780072/1‑250 (MQ=255)
                                                                                                          cA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCa                                                                                                                                >  2:623757/1‑250 (MQ=255)
                                                                                                          cA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCa                                                                                                                                >  2:724335/1‑250 (MQ=255)
                                                                                                                        ccAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCACAACCAATACCTa                                                                                                                   >  2:76298/1‑250 (MQ=255)
                                                                                                                        ccAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTa                                                                                                                   >  1:884794/1‑250 (MQ=255)
                                                                                                                         cAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTAt                                                                                                                  >  2:429468/1‑250 (MQ=255)
                                                                                                                               cAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTccccc                                                                                                            >  1:528695/1‑250 (MQ=255)
                                                                                                                                aGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcccccc                                                                                                           >  1:244628/1‑250 (MQ=255)
                                                                                                                                aGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTcccccc                                                                                                           >  1:368650/1‑250 (MQ=255)
                                                                                                                                           tCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTACCCCCAATCTgccgcc                                                                                                >  1:19950/1‑250 (MQ=255)
                                                                                                                                             aGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCa                                                                                              >  1:356372/1‑250 (MQ=255)
                                                                                                                                               gAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaag                                                                                            >  2:598283/1‑250 (MQ=255)
                                                                                                                                                        cACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGca                                                                                   >  2:78936/1‑249 (MQ=255)
                                                                                                                                                         aCGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGTGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCa                                                                                  >  1:191313/1‑250 (MQ=255)
                                                                                                                                                          cGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAg                                                                                 >  1:429482/1‑250 (MQ=255)
                                                                                                                                                                  gCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTa                                                                         >  1:736414/1‑250 (MQ=255)
                                                                                                                                                                             gATGTATTTCGGGTTACGGCGCAGCAGCCTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCt                                                              >  1:273649/1‑250 (MQ=255)
                                                                                                                                                                                   tttCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATATGGCCCCCAAGAAGACCGCCAGCAAAACAATGATAAGATCTACactt                                                        >  2:532266/1‑247 (MQ=255)
                                                                                                                                                                                    ttCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTa                                                       >  2:604271/1‑250 (MQ=255)
                                                                                                                                                                                        ggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCat                                                   >  1:848568/1‑250 (MQ=255)
                                                                                                                                                                                           tACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGACGATCGCCAGCAACACTATGATGAGTTCTACAACCAGCATAtt                                                >  2:367448/1‑250 (MQ=255)
                                                                                                                                                                                              ggCGCAGCAGCTTTTCTGTCTGATGAACCAGATATTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAATTATGATGAGTTCTACAACTAGCATATTAGc                                             >  1:889451/1‑250 (MQ=255)
                                                                                                                                                                                                               gtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATT‑CCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACa                           >  2:406080/1‑250 (MQ=255)
                                                                                                                                                                                                                                gTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTg           >  2:644945/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt         >  2:925437/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTa        >  1:202482/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTa        >  2:116516/1‑250 (MQ=255)
                                                                                                                                                                                                                                        caccggCAACCTGCATGGCAG‑‑AAATGGCGGAGAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttttt   <  1:610830/245‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGTTGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTa  >  1:315075/1‑250 (MQ=255)
                                                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                      
ATGCGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGC‑CTTGTT‑CCTTCGCAACTTCAGCGATAACT‑GGCAGTGTCGCCA‑GAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTA  >  NC_007779/4371022‑4371502

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: ‑