Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,138 G→A 11.5% pseudogene (747/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,1380GA11.5% 270.2 / 16.2 108pseudogene (747/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (26/66);  new base A (5/7);  total (32/76)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

GCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCC  >  NC_000913/1465998‑1466285
                                                                                                                                            |                                                                                                                                                   
gcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                                                                                                                                            <  1:343944/150‑1 (MQ=255)
         gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc                                                                                                                                   >  1:754395/1‑150 (MQ=255)
         gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc                                                                                                                                   <  2:205281/150‑1 (MQ=255)
          cAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATACCGGGAAAGTGATTCAGGACGGCg                                                                                                                                  >  1:495532/1‑150 (MQ=255)
             ataaGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATc                                                                                                                               <  1:504659/148‑1 (MQ=255)
              gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt                                                                                                                              <  1:1645133/150‑1 (MQ=255)
               aaGTGGCAATATAATCACAACGTTGAGCTTGTTATTACGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATTACGGGAAAGTGATTCAGGACGGCGATCTg                                                                                                                             <  1:196110/150‑1 (MQ=255)
                aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGAAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                            >  1:538626/1‑150 (MQ=255)
                   ggCAATATAATCACAACGGTGAGCTTGTTATTACTGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATg                                                                                                                         <  2:1890197/150‑1 (MQ=18)
                          taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTTAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                  >  1:74838/1‑150 (MQ=255)
                          taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                  >  2:475045/1‑150 (MQ=255)
                           aaTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggc                                                                                                                 <  2:740681/150‑1 (MQ=255)
                             tCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                                                                                                                           <  2:1763046/138‑1 (MQ=255)
                             tCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGa                                                                                                                           >  1:1763046/1‑138 (MQ=255)
                              cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc                                                                                                              >  2:1741449/1‑150 (MQ=255)
                               acaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                                                                                                             <  1:1892815/150‑1 (MQ=255)
                                caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                            <  1:398465/150‑1 (MQ=255)
                                caAAGGTGAGCTTTTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACTGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                            <  2:1176621/150‑1 (MQ=255)
                                aaaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATCGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGACATCATCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                            <  1:89509/149‑1 (MQ=255)
                                 aaCGGTGAGCTTGTTAATACGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGCTGACGGTAAAGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCTATCTGGATGTCAGCGGCGGCGGt                                                                                                           <  1:590522/150‑1 (MQ=18)
                                     gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                       <  1:1560182/150‑1 (MQ=255)
                                      tGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                      <  1:357348/150‑1 (MQ=255)
                                         gCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                   >  2:1333778/1‑150 (MQ=255)
                                           ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                 <  1:475045/150‑1 (MQ=255)
                                           ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                 <  1:944890/150‑1 (MQ=255)
                                           ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                 <  1:859182/150‑1 (MQ=255)
                                           ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                 <  1:1612738/150‑1 (MQ=255)
                                             gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                               >  1:1315851/1‑150 (MQ=255)
                                              ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                              <  1:1180766/150‑1 (MQ=255)
                                                attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                            <  2:1267993/150‑1 (MQ=255)
                                                  taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAc                                                                                          >  1:1559804/1‑150 (MQ=255)
                                                    cGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCg                                                                                        >  1:1599725/1‑150 (MQ=255)
                                                     gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                       >  2:1346844/1‑150 (MQ=255)
                                                     gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                       <  2:1709217/150‑1 (MQ=255)
                                                     gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                       <  2:337767/150‑1 (MQ=255)
                                                         taTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGAc                                                                                   <  2:1184949/149‑1 (MQ=255)
                                                         gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGAc                                                                                   <  1:1436366/150‑1 (MQ=255)
                                                         gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGAc                                                                                   <  2:789336/150‑1 (MQ=255)
                                                         gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                                                                                                                       >  1:965911/1‑114 (MQ=255)
                                                         gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                                                                                                                       <  2:965911/114‑1 (MQ=255)
                                                             aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAACAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAgcg                                                                               >  1:1735035/1‑150 (MQ=255)
                                                              tGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgcgc                                                                              <  1:563830/150‑1 (MQ=255)
                                                                  acaGTAAACAACAATGGCAAAACCACCGTTGACGGGAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                           <  1:747697/85‑1 (MQ=11)
                                                                  acaGTAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTc                                                                                                                                           >  2:747697/1‑85 (MQ=11)
                                                                    aGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                                                           <  1:207431/131‑1 (MQ=12)
                                                                    aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCa                                                                                           >  2:207431/1‑131 (MQ=255)
                                                                     gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGTAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCTGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTg                                                                       <  1:498816/150‑1 (MQ=255)
                                                                      tAAACACCAATCGCAAAACCACCGTTGACGGTAAAGACTCCACCGGTACGGAAAGCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                      <  1:1715780/150‑1 (MQ=14)
                                                                       aaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGa                                                                     <  2:528651/150‑1 (MQ=255)
                                                                        aacaacAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                    >  2:1485759/1‑150 (MQ=12)
                                                                        aacaacAATGGCAAAACCACCGTAGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                    <  1:395981/150‑1 (MQ=255)
                                                                               aTGGCAAAACCTCCGTTGACGTTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                             <  2:1715484/150‑1 (MQ=255)
                                                                               aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                             <  1:897701/150‑1 (MQ=255)
                                                                               aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAgg                                                             >  2:1333472/1‑150 (MQ=25)
                                                                                   cAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAc                                                         >  2:629311/1‑150 (MQ=255)
                                                                                    aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                        <  2:316332/150‑1 (MQ=255)
                                                                                     aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATCACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCa                                                       <  1:1183056/150‑1 (MQ=255)
                                                                                     aaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCa                                                       <  1:1465285/150‑1 (MQ=255)
                                                                                     aaaCCACCGTTGACGGTAAAGAATCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCa                                                       <  1:117002/150‑1 (MQ=255)
                                                                                         caccGTTGACGGTAAAGCTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAc                                                   <  1:1738720/150‑1 (MQ=255)
                                                                                           ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtta                                                                                                                                              <  2:664140/57‑3 (MQ=255)
                                                                                           ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGtta                                                                                                                                              >  1:664140/1‑55 (MQ=255)
                                                                                            cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGt                                                <  1:1760178/150‑1 (MQ=255)
                                                                                             gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTc                                               <  2:1592877/150‑1 (MQ=255)
                                                                                              ttGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACA‑‑GCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                            <  1:1599675/150‑1 (MQ=255)
                                                                                                gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                            <  1:209918/150‑1 (MQ=255)
                                                                                                  cGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACAGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGa                                          <  1:849566/150‑1 (MQ=255)
                                                                                                   ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                         <  1:1509999/150‑1 (MQ=255)
                                                                                                     tAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                       >  1:408622/1‑150 (MQ=255)
                                                                                                      aaaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgca                                                                             >  2:804648/1‑110 (MQ=21)
                                                                                                       aaGATTCCACCGGTACGGAAATCAACGGTACTAACGGGAAAGTGATGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCTATCCgg                                     <  1:769026/150‑1 (MQ=255)
                                                                                                       aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                     <  2:191608/150‑1 (MQ=40)
                                                                                                         gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                   <  2:217263/150‑1 (MQ=39)
                                                                                                         gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATGCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                   <  2:1021669/150‑1 (MQ=255)
                                                                                                          aTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                  <  2:1315851/150‑1 (MQ=39)
                                                                                                          aTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                  <  2:615620/150‑1 (MQ=39)
                                                                                                             ccACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCa                               <  2:953158/150‑1 (MQ=39)
                                                                                                             ccACCGGTACGGAAATCAACGGAAATAACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAAAAAGGGCACCATGACCGTCACCGATCCGGAGTCCa                               <  2:320023/150‑1 (MQ=255)
                                                                                                              cACCGGTACTGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCGGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAt                              <  1:334960/150‑1 (MQ=255)
                                                                                                               aCCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATg                             <  1:1087382/150‑1 (MQ=40)
                                                                                                                ccGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                               <  2:1650260/83‑1 (MQ=18)
                                                                                                                ccGGTACGGAAATCAACGGTAATAACGGCAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                               >  1:1650260/1‑83 (MQ=18)
                                                                                                                    tACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTAt                        >  1:448110/1‑150 (MQ=255)
                                                                                                                     aCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATc                       >  2:1025160/1‑150 (MQ=255)
                                                                                                                      cGGAACTCAACGCTAATAACGGGAAAGTGATTCAGGACGGCGATCGTGATGTCAGCGGCGGCGGTCACGGTAGAGATATCACCGGTGACAGCGCGACGGTGGATAACAATGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATcc                      <  2:283014/150‑1 (MQ=14)
                                                                                                                          aaTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGAt                  <  2:829625/150‑1 (MQ=255)
                                                                                                                            tCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                <  1:1333778/150‑1 (MQ=255)
                                                                                                                             cAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                    >  1:1544305/1‑145 (MQ=21)
                                                                                                                             cAACGGTAATAACGGGAAAGTCATTCAGGACGGCGATCTGGATGTCAGCGTCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGa               <  1:304436/150‑1 (MQ=255)
                                                                                                                               aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACg             <  2:1636108/150‑1 (MQ=255)
                                                                                                                                cGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATTCAGATCGACgg            >  2:1727596/1‑150 (MQ=255)
                                                                                                                                cGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACgg            <  1:91189/150‑1 (MQ=255)
                                                                                                                                  gTAATACCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTg          <  1:366715/150‑1 (MQ=255)
                                                                                                                                  gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTg          <  2:1706588/150‑1 (MQ=255)
                                                                                                                                      taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                      <  1:464569/54‑1 (MQ=35)
                                                                                                                                      taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg                                                                                                      >  2:464569/1‑54 (MQ=35)
                                                                                                                                      taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACaa      >  2:203305/1‑150 (MQ=255)
                                                                                                                                       aaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAg     >  1:191953/1‑150 (MQ=255)
                                                                                                                                         cGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACTGGATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGc   <  1:314397/150‑1 (MQ=255)
                                                                                                                                          gggAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGcc  <  1:761402/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCC  >  NC_000913/1465998‑1466285

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