Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,856 A→G 11.7% 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,8560AG11.7% 121.8 / 8.1 53pseudogene (1465/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base G (2/4);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

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