Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,856 A→G 12.4% pseudogene (1465/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,8560AG12.4% 113.3 / 8.3 48pseudogene (1465/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base G (2/4);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGT  >  NC_000913/1466635‑1467083
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
ttaaggggtaggaCGCCACGGCAAACAACAACGGAAAAGCCACCTTTGACGGCAAGGATTCCAACGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACg                                                                                                                                                                                                                               <  2:526509/217‑1 (MQ=14)
            aCGCCACGGCAAACAACAACGGCAAAACCACCGGTGACGGCAAGGATTCCACCGGTATGGAAATTGCTGGTAATGACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCCCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCa                                                                                                                                                                                                                   <  1:659883/228‑1 (MQ=255)
                            caacGGCAAAACCAGCGTTGGCGGCAAGTATTCCACCGGTACGGAAATTGATGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAgg                                                                                                                                                                                                   <  1:580078/228‑1 (MQ=255)
                                    acaccccggTTCACGGAAAGTATTCCATCTGTACAGAACATGCTGGCAATAACGGGAAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAACGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                                                                                                                                                                           <  1:255561/220‑1 (MQ=17)
                                       ccccctttGACGGTAAGATTTCCAACGGTACCGAAATAAACGGAAATAACGGGAAAGTGATTCAGGACGGCGACCTGGATGTCAGCGGCGGCGGTAACGGTATTGATAACACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTa                                                                                                                                                                                        <  2:513205/222‑1 (MQ=11)
                                                       ggATTCCCCCGGGCCGGAAATTGGCGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                                                                                                                        <  1:548064/228‑1 (MQ=255)
                                                               cctgTACGGAAATGGGGGGCAAATACGGGAAGGTTGCTCGGGACGGCGGTCTGGATGTCAGCGGCAGAGGTCACGGTATTGATATCCCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt                                                                                                                                                                 <  1:8464/224‑1 (MQ=17)
                                                                 ggTACGGAAATTGCTGGCAATAACGGGAGGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAg                                                                                                                                                              <  1:263398/228‑1 (MQ=255)
                                                                  gTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATAACCAATGGcggcaccggccccccg                                                                                                                                                              >  1:448876/1‑221 (MQ=255)
                                                                   tACGGAAATTGGTGGCAATAACGGGAAGGTGATTCAGGACGGGTATCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt                                                                                                                                                                 <  2:448876/223‑1 (MQ=255)
                                                                   tACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGGATCCAGATTGACGGCGACCAGGCCATGTCCAATAACGAAGGCGAGAGACCTATAACCAATGGcggcaccgggccacaaaa                                                                                                                                                            >  2:632025/1‑219 (MQ=17)
                                                                                    attaCGGGAAGGTGATTCAGGTCGGCGTGTTGGAGGTACGTGGCGGCGTTACCGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCa                                                                                                                                           <  1:281808/225‑1 (MQ=17)
                                                                                    aTAACGGAAAAGTTTGTCAGGAAGGGGATCTGGATGGCAGCGGCGGCGGGCACGGTATTGATCTCACCGGCGACTGCGCAGCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCa                                                                                                                                           <  2:652820/228‑1 (MQ=14)
                                                                                          ggAAGGTGATTCAGGTCGGCGAGCTGTAGGTACGCGGAGGAGATCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgcg                                                                                                                                      <  1:392683/227‑1 (MQ=17)
                                                                                                ggATTCAGGGCGGGGTTTTGGGTGTCATGGGGCGCGGTCACGGTTTTGAAACCACCGGCCACAGCGCAACGGGGGATAACAAGGGCTCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCAGCCAGGCCATCGTCAATAAGGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACa                                                                                                                               <  2:594216/227‑1 (MQ=17)
                                                                                                     cAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCACCGTCCATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCCCTCAGAACAACGGTAAAGACCCGACCCCGAATACCAGTGGa                                                                                                                          >  1:586194/1‑228 (MQ=17)
                                                                                                          cGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCCTCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGAGAGAGCAATCTCACCAATGGCCGCACCCGCACTCCGATCAACGGGAACGACGCCCCCGCGAATAAAAGAGGAAAAAc                                                                                                                     >  2:336107/1‑228 (MQ=17)
                                                                                                            gggATCTGGAGTTCAGCGGGGGGGGTCACGGTTTTGATATCACCGGCGACCGGGAAACGGTGGAAAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAcc                                                                                                                    <  2:11891/225‑1 (MQ=18)
                                                                                                            gCGATCTGGATGTCAGGGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAc                                                                                                                     <  1:422055/226‑1 (MQ=255)
                                                                                                                    gATGGCAGCGGCGGCGGTCACGGTATTGATATCCCCGGCGACAGCGCAACGGTGGATAACAAGGGCCCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGAt                                                                                                           <  2:242124/228‑1 (MQ=255)
                                                                                                                          aGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGAAACGACGCCACCCCGAGTAAAAGTGGAAAAACCCCTGTTGATGGaaaa                                                                                                     >  1:569154/1‑228 (MQ=14)
                                                                                                                               cggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAACACTGTTGGTGTAAAAGATTc                                                                                                >  1:360486/1‑228 (MQ=255)
                                                                                                                                                atatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGGACCAAAATcg                                                                               >  1:404897/1‑228 (MQ=255)
                                                                                                                                                     aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAcc                                                                                                                                                                                                                         >  1:171459/1‑85 (MQ=33)
                                                                                                                                                     aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGCCCATCGTCAATAACGTAGGCGAGAGCACTATCACCAATGGCGGCACCGGCCCTCAGATCAACGGTAACGCCCCCACACCGATTAACAGGGCACAAACACCTGTTTTAGGAACAGATTACACGGGTACCAAAATCACggga                                                                          >  2:180149/1‑227 (MQ=18)
                                                                                                                                                      ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCa                                                                                                                                                                                                                        <  2:171459/85‑1 (MQ=33)
                                                                                                                                                       cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGCGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCAACGGCACTCAGACCGACGGCAACGACGCCACCCCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAAATTCCCCGGGACCAAAAACGCGGGCaaa                                                                       >  2:500480/1‑227 (MQ=21)
                                                                                                                                                       cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGTGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGCTGGAAAAGATTCCACAGGTACCAAAATCGCGGGCaa                                                                        >  1:446017/1‑228 (MQ=255)
                                                                                                                                                       cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAACAAGGAAAAACACATGTTGATGGAAAAGATTCCACGGGGACCAAAAtaacggaaaa                                                                        >  2:337568/1‑219 (MQ=21)
                                                                                                                                                             cAGCGCAACGGATGATAACAAGGGCACCATGACCGTCGCCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAACAGGCCATCGTAAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCgg                                                                  <  1:604214/228‑1 (MQ=255)
                                                                                                                                                                cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGCATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCAt                                                               >  1:386143/1‑228 (MQ=255)
                                                                                                                                                                cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCCCGGGTACCAAAATCGCGGGCAATATCGGCAt                                                               >  1:70285/1‑228 (MQ=255)
                                                                                                                                                                cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGAGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAAACCCGGGAAAAATCGGCAt                                                               >  2:598775/1‑228 (MQ=255)
                                                                                                                                                                      ggTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATCGGTAGCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAATATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTaa                                                          <  1:336107/227‑1 (MQ=255)
                                                                                                                                                                        tGGATAACAAGGGAACCCGGACCGGCCCCGATCCGGAGTCCAGCGGTATCCAGATTGACGGCGACCAGGCCAGCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGACACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                        <  2:586194/227‑1 (MQ=255)
                                                                                                                                                                         ggtaaaCAAGGGCCCCCTGGCCGGCCCCGAATCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCt                                                      <  2:569154/224‑1 (MQ=18)
                                                                                                                                                                            tAACAAGGGCACCATGGCCCTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGGCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCt                                                      <  1:500480/225‑1 (MQ=255)
                                                                                                                                                                             aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACCGTTGATGGAAAAGATTCCACGGGGACCAAAATCGCGGGCAATAACGGCATTGAAAAACTGGaa                                                  >  2:18161/1‑227 (MQ=255)
                                                                                                                                                                               cAAGGGCCCCCAGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAAGATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATgg                                                <  1:484789/228‑1 (MQ=255)
                                                                                                                                                                                  gggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCCCGGGCCATATCGGCATTTGAAATCTGGAGGGTAg                                             >  1:588025/1‑228 (MQ=255)
                                                                                                                                                                                  gggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACAAAAATCGCGGGCAAAATCGGCATTTCAACTTTGGATGGTa                                              >  1:405464/1‑227 (MQ=255)
                                                                                                                                                                                    gccccATGACCGTCACCGATCCGGAGTCCATCGGTCACCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCCCCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGc                                            <  1:455061/224‑1 (MQ=255)
                                                                                                                                                                                       ccATGGCCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCGCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTg                                         <  1:180418/227‑1 (MQ=255)
                                                                                                                                                                                                tCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAATATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGGGAATAACAGGGGGAAAACCAATGTTAATGGAAAAGGTCCACCGGGAACCAAAATTGCGGGCAAAATCGGCATTGTAAAGATGATAGGAACCCTGACTGTTACAg                               >  2:610127/1‑228 (MQ=21)
                                                                                                                                                                                                  aCCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCTGCATTGGAAATCTGGATGGTAGCCTGACTGTTAAAGGc                             >  2:77458/1‑228 (MQ=255)
                                                                                                                                                                                                   ccGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGGCAATAACGAAGGCGAGAGCACTATCACCAAAGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCGCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCggg                                                                           <  1:204449/181‑1 (MQ=255)
                                                                                                                                                                                                           gAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGAAGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAACATGGATGGTAGCCTGACTGTTACAGGCGGGGGGCa                     >  1:404319/1‑227 (MQ=255)
                                                                                                                                                                                                                cATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGCCAATATCGGCATTGTAAATCTGGAGGTGAGCCTGACTGTTCCAGGCGGGGCGCCTGgggt               >  1:344330/1‑228 (MQ=37)
                                                                                                                                                                                                                             aTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGGGTTTAGAACATTggg  >  1:76904/1‑227 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
TTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGT  >  NC_000913/1466635‑1467083

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