Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 514,434 C→T G142G (GGG→GGA ybbJ ← NfeD‑like family protein YbbJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,4340CT92.3% 98.7 / ‑3.2 40G142G (GGG→GGAybbJNfeD‑like family protein YbbJ
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (18/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGAT  >  NC_000913/514193‑514676
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
taaTCTGTTTAACGACCCGCAGCGGCACAGCGCCTACGTCAAACGGCGCACGCGGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGATCCGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCTACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGTTTATTCCTTCTAt                                                                                                                                                                                                                                            <  2:195558/249‑1 (MQ=255)
       tttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATg                                                                                                                                                                                                                                     >  1:298117/1‑250 (MQ=255)
        ttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGa                                                                                                                                                                                                                                    <  1:320763/250‑1 (MQ=255)
              aCCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAAc                                                                                                                                                                                                                             >  1:217569/1‑251 (MQ=255)
                 cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATg                                                                                                                                                                                                                          >  1:41680/1‑251 (MQ=255)
                 cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATg                                                                                                                                                                                                                          <  1:312729/251‑1 (MQ=255)
                          acaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCCCTGGCGAATGCCTGCCCTGGCGATGACAGAGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGTGGAGGAGGATTAAGACGAGACGGCCCGGAGATGGCGGGTTCATCCTTCTACCGCAATGACTTCTACATGTGTACCTGc                                                                                                                                                                                                                 >  2:85221/1‑251 (MQ=255)
                                 cGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg                                                                                                                                                                                                             >  1:27908/1‑248 (MQ=255)
                                  gACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCatcg                                                                                                                                                                                                             >  1:23423/1‑244 (MQ=255)
                                                      ggAGAAGGTGGCTGGGATCGAACGACACCTTGAGGAGTTGCACTCCAAGCGCTCCCAGTGGATGTCCCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTTTCGAAAAACTCTCCGGCTGCTGTCATCACCGGGCAGGGTGTGTAAGGCGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTATCTGCACCGAGATCCTCGCTGTCGCTGACAGGc                                                                                                                                                                                     <  2:111013/251‑1 (MQ=255)
                                                            gggggCGGAGATCGAACGACACATTGAGGAGCTGCAAACCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGaa                                                                                                                                                                               <  2:27908/248‑1 (MQ=255)
                                                                  ggAGATCGAACGACCCATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCGCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCa                                                                                                                                                                         <  2:27206/251‑1 (MQ=255)
                                                                   gagaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACGGTCAc                                                                                                                                                                        >  1:371473/1‑251 (MQ=255)
                                                                    agaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGCCCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGAGCTGACAGCGCCGCCTGACCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGAGATGTAGGGATATTCCTTCTAACGCAATGACTTCAACATGCGAACCTGCACCGAGATCCTCGCTGGAGCTGCCCGGCCATGTACTGtaccc                                                                                                                                                                       >  2:23047/1‑247 (MQ=255)
                                                                     gaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCa                                                                                                                                                                                                                                        >  1:304945/1‑185 (MQ=255)
                                                                     gaTCGAACGACACATTGAGGAGCTGCAACCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCa                                                                                                                                                                                                                                        <  2:304945/185‑1 (MQ=255)
                                                                                acaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGACTGCATGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGAGCTGACAGGCCATGAACTGTCACCGACGCGCATATg                                                                                                                                                           >  1:232/1‑251 (MQ=255)
                                                                                           cTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCACCATCGGGCAGGGAGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAATATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGt                                                                                                                                                >  1:155791/1‑251 (MQ=255)
                                                                                                     tGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGTCCAGCgg                                                                                                                                      >  2:135247/1‑251 (MQ=255)
                                                                                                     ggcgggacaggatgcgggcctgggCGAATGCCTGCCCTGGCGAGGACAGCGCCGACTGCCCGATTAGCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGGGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCGGACAGGCCATGAACTGTCAACGACGCGCATATGACCGCGCCCGTTGACCAGCgg                                                                                                                                      <  2:68723/229‑1 (MQ=255)
                                                                                                       cgcgACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGAAGACAGCACCGACCGCCAGCATATCGAAAATCTCTCAGGCTGCTGTCAACATCGGCCAGGGTGATTAAGACGCGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCCACATCCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGCACTGTCAACGACGCGCATCTGACCGCACCCGTTGACCAGCggtg                                                                                                                                    >  1:363038/1‑249 (MQ=255)
                                                                                                         cgACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGAt                                                                                                                                  <  1:5842/251‑1 (MQ=255)
                                                                                                          gACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGCTGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCCTCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCCACATGCGTACCTGCACCGAGCTCCTCACTGCCGCTGACAGGCCATGAACTGTCACCCACGCGCATAGGACCGCGCCCCTTGACCAGCGGAGAtt                                                                                                                                 >  2:378564/1‑251 (MQ=255)
                                                                                                               gctgctGGCACTGGCGAATGCCCGCCCTGGCGATGACAGCGCCGACTGCCCGATTCGCGAAAATCTCGCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGc                                                                                                                            <  2:59486/251‑1 (MQ=255)
                                                                                                                           ggCGAATGCCTGCCCTGGCGATGACAGCGCCGCCTGCCCGATTAACGAAAATCTCTCCGGCTGCGTCCATCATCGGCAAGGGTGATTAAGACGAGTCGGCTCGGATATGTAGGGTTATGCCTTCTATCGCTATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTc                                                                                                                  <  2:152596/249‑1 (MQ=255)
                                                                                                                                aTGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAAt                                                                                                           <  1:338289/251‑1 (MQ=255)
                                                                                                                                   ccggccatggggATGACAGCGGCGACTGACCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGGGATTAAGACGAGACGGCTCGGATGTGTAGGGTTATTCCTTCTATCGCAATGACTACAACATGCGTCCCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCATAGGAGATTCCAGCACAAATCGTCGGCAAATCAg                                                                                                        <  2:303116/240‑1 (MQ=255)
                                                                                                                                               atgacatcgccgccgtcgcggtgaTCGAACATCTCTCCGGCTTAGGTCATGATCGGGCAGTGTGATTAAGACGAGACTGCTCGGATATGTAGAGGTATTCCTGCTATCGCAATGACTTCAACTTGCCTACCTGCGCCGAGATCCTCGCTGGCGCTGACAGGCCTGGAACTGTCACCGACCCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCTCAAAACGTCGGCCAATCAGCTGCTGcca                                                                                               <  2:31242/225‑2 (MQ=255)
                                                                                                                                                 gACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATTCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGCCGGCCAATCAGCTGCTGCCCGCGCt                                                                                          <  1:139191/251‑1 (MQ=255)
                                                                                                                                                                  gtttATCGAAAATCCTTTCGGCTGTTGTCATTATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATTACTTCAACAGGCGTAACTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGTCCGCGCTGGTTTAAATGACTATCg                                                                         <  2:280160/249‑1 (MQ=255)
                                                                                                                                                                     tATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTTCTGCCCGCGCTGGTTTAAATGACTATCGCt                                                                       <  1:159112/250‑1 (MQ=255)
                                                                                                                                                                            aaTCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCtttt                                                               <  1:116692/251‑1 (MQ=255)
                                                                                                                                                                                           tgTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgcc                                                >  1:66448/1‑251 (MQ=255)
                                                                                                                                                                                                           gggTGATTAAGACGAGCCGGCTCAGATATGTAGGGATATTCCTTCTATCGCAATGACTTCAAAATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCCTGCACTGTCACCGACGCGCATATGGCCGCGCCCGTTGACCAGCGCAGATACCAGCACCAAACGTCGTACAACCACCTGAGTCCCGGGCTAGGTTAAATCACTCACGATGTGCTTTGGTTCGCCAAACCGCCGCGACCACCATCTcc                                 >  2:374571/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCATCACTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCAGCGACAACAATTTccacaac                            >  1:163094/1‑247 (MQ=255)
                                                                                                                                                                                                                  tAAGCCGAGACGGCTCGGCTATGTAGGGTTATTCCTTCTATCGCACTGCCTTCAACATGCGTCCCTGCCCCGAGATCCTCGCTGGCGCTGACAGGCCATGAACGATCACCGCCGCGCATATGACCACGCCCCTACTCGAGCGGAGCTTCCAGCACCACCCGTCGGCCCATCAGCTCCCGCCCGCGCTGGGTTAAATGCCAATCACTTAGCTTTTGTTCGCGCACCCACCCCCACCACCATTTCccccaccg                          >  2:72846/1‑249 (MQ=255)
                                                                                                                                                                                                                         aggcgGCGCGGATTTGTAGGTTTATTCCTTTTACCGAAATGACTCCAACAGGCGACCCTCCCCCGGGATCCCCGGGGGCGCGGCCAGGCCATAAACGGCCCCCGCCGCGAATATGACCGCGCCCGTTGACCAGCGGAGATCCCAGCAAAAAACGCCGGCAAATCAGCTGCTCCCCGCGCTGTTTTAAATGACTATCGCTGTGCTTTTGTTCGCCCACCCGCCGCGACAACCATTTCCACCACAGCCAggcg                   <  1:135247/248‑1 (MQ=255)
                                                                                                                                                                                                                            cgGCGCGGATATGTAGGGTTATCCCTTCTATCGCAAGGCCTCCAACACGCGTCCCTGCCCCGAGACCCTCGCGGGCGCTGACAGGCCATGAACTGCCACCGACGCGCATATGACCGCCCCCGTGGACCAGGGGAGATCCCAGCAAAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGTTTTAAATGCCTACCGCTGTGTTTTTGTCCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGgcggcg                <  1:195118/251‑1 (MQ=255)
                                                                                                                                                                                                                              gctaggATATGTAGGGGTATTCCTTATATAGCAATGACGTCAACATGCGTACCTGCACCGAGATCCTCGATGGCGCTGACAGGCCATGAACTGTCACCGACTCGCATATGACCGCGCCCGTTGACCAGAGGAGATTCCAGCACAAAACGGCGGCCAATCATCTGCTGACCGCGCTGGTTTAAATGACTATCGCTGAGCTTTGGATCGCGCACCCGCCGAGAACACCATTTCCACCACAGCCAggcacgc                >  2:124071/5‑245 (MQ=255)
                                                                                                                                                                                                                                          aGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGAt  >  1:94095/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
GAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGAT  >  NC_000913/514193‑514676

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