Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,144 G→T 8.8% pseudogene (753/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,1440GT 8.8% 338.5 / 11.6 130pseudogene (753/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (45/72);  new base T (7/6);  total (52/78)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT  >  NC_000913/1465995‑1466290
                                                                                                                                                     |                                                                                                                                                  
aGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACTGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAtgtg                                                                                                                                                    <  2:1579660/150‑1 (MQ=255)
 gcgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACACCAATGGCATAACCACAGTTGACGGTAACGATTCCCCAGGTACGGAAATCAAACGTATAAACGGGAATGTGa                                                                                                                                                   >  1:342930/1‑150 (MQ=255)
  cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAt                                                                                                                                                  <  1:737305/150‑1 (MQ=255)
     aGATGTCGCAAAGAAGTGGCAATATAATCACTACTGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCa                                                                                                                                               <  2:763750/150‑1 (MQ=21)
      gATGGCGCAACGAAGGCGCAATATAATCACAACGGTCATCTTGTTATTACGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAg                                                                                                                                              <  2:144178/150‑1 (MQ=25)
         ggCGAAACGAAGTGGAAATATAATCACAACGGTTAGCTTGTTATTACGTGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGAc                                                                                                                                           <  2:1922069/150‑1 (MQ=255)
          gcgcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACg                                                                                                                                          >  1:1659333/1‑150 (MQ=255)
            gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGAAAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc                                                                                                                                        <  2:130650/150‑1 (MQ=255)
               accaaGTGCCAATATAATCACAACGGTGAGATAGTTATTACGGGTGATAATGCCACAGTAAACAAAAATGGAAAAACCACCGTTTACGGAAAAGATTCCACCGGTACTGAAATCAACTCTAATAACCGTAAAGTGATTCAGCACGGCGaa                                                                                                                                     >  2:639613/4‑149 (MQ=255)
               aCGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt                                                                                                                                     >  2:1083836/1‑150 (MQ=255)
                cGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATc                                                                                                                                    <  1:490848/150‑1 (MQ=255)
                 gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt                                                                                                                                   >  2:1629128/1‑150 (MQ=255)
                  aaGTGGCAATATAATCTCAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAAGAACATTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTg                                                                                                                                  <  1:133499/150‑1 (MQ=255)
                   aGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                                 <  1:1404485/150‑1 (MQ=255)
                      ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGTTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                                                                                                              <  1:1863493/150‑1 (MQ=255)
                      ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTGAGGACGGCGATCTGGATg                                                                                                                              <  2:371200/150‑1 (MQ=255)
                      ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                                                                                                              <  1:1226578/150‑1 (MQ=255)
                      ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAAAGGTAATAACGGGAAAGAGATTCAGGACGGCGATCTGGAtt                                                                                                                              >  2:515912/1‑149 (MQ=255)
                        cAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTc                                                                                                                            >  2:2020777/1‑150 (MQ=255)
                             taATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                       >  1:1339321/1‑150 (MQ=255)
                              aaTCAAAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAACCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggc                                                                                                                      <  1:806536/150‑1 (MQ=255)
                                 cactACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAATATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATGCAGGCCGGCGATCTGGATGTCAgccgcggc                                                                                                                   >  2:362436/1‑150 (MQ=21)
                                   caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                                 <  1:574700/150‑1 (MQ=255)
                                   caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                                                                                                                 >  2:883218/1‑150 (MQ=255)
                                      cGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCTATGAAGCTAAAGATTCAGCAGGAAAGGAAAACACAGATAAAAACTGCAATGTGATGCAGGACAGAGATCTGGATGTGAGCGGCGGCGGTcc                                                                                                              >  1:100587/1‑149 (MQ=255)
                                        gTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                            <  1:1326865/150‑1 (MQ=255)
                                          gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt                                                                                                          <  2:1541292/150‑1 (MQ=255)
                                            gCTTGTTATTACGGGTGATAATGCCACAGTGAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                        >  1:831470/1‑150 (MQ=255)
                                              ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                      >  2:742401/1‑150 (MQ=255)
                                              ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                      <  2:1329713/150‑1 (MQ=255)
                                                gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                                    <  2:1356381/150‑1 (MQ=255)
                                                 ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                   <  1:2020777/150‑1 (MQ=255)
                                                 ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                   <  1:640267/150‑1 (MQ=255)
                                                 ttattaCGGGGGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata                                                                                                   <  1:1111966/150‑1 (MQ=255)
                                                  tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat                                                                                                  >  1:1431295/1‑150 (MQ=255)
                                                  tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat                                                                                                  >  1:1989450/1‑150 (MQ=255)
                                                   attaGGGGTTATAATGCCAGAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAATATTCCACCTGTACGGAAATCAACGGTAATATCTGGAAAGTGATTCAGGACTGCGATCTGGATGTCAGCGGCGGCGGTGACGGTATTGATATc                                                                                                 <  2:115233/150‑1 (MQ=17)
                                                   attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                 <  1:1455217/150‑1 (MQ=255)
                                                   attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                 <  1:1119893/150‑1 (MQ=255)
                                                      aCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                              >  1:190121/1‑150 (MQ=255)
                                                        gggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                            <  1:570139/150‑1 (MQ=255)
                                                          gTGATAATGCCACAGTAAACAACACTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTg                                                                                          <  1:842994/150‑1 (MQ=255)
                                                               aaTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAgc                                                                                     >  1:818930/1‑150 (MQ=255)
                                                                aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa                                                                                                                                      >  1:455977/1‑100 (MQ=255)
                                                                aTGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGAAATGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGa                                                                                                                                      <  2:455977/100‑1 (MQ=11)
                                                                     acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGAGACAGCGCGACg                                                                               >  1:1903547/1‑150 (MQ=255)
                                                                     acaGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                               <  1:859453/150‑1 (MQ=255)
                                                                        gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTg                                                                            <  2:501395/150‑1 (MQ=255)
                                                                         tAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGAAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTgg                                                                           <  2:1438613/150‑1 (MQ=255)
                                                                             caacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                                                                                                                    >  2:1277588/1‑89 (MQ=255)
                                                                              aacaaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAAc                                                                      <  2:1374714/150‑1 (MQ=255)
                                                                                caaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACAGTATTGATATc                                                                                                 <  2:617451/121‑1 (MQ=12)
                                                                                  aTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgca                                                                                  >  1:2067549/1‑133 (MQ=255)
                                                                                    ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGc                                                                <  2:2048566/150‑1 (MQ=255)
                                                                                       aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGTATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCAcc                                                             <  2:1237018/150‑1 (MQ=255)
                                                                                           ccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGa                                                         >  1:1910334/1‑150 (MQ=255)
                                                                                            caccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATc                                                                                                 <  2:2088797/109‑1 (MQ=12)
                                                                                              ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAccc                                                      >  2:922555/1‑149 (MQ=255)
                                                                                              ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                      >  2:1839538/1‑150 (MQ=255)
                                                                                              ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                      <  1:395434/150‑1 (MQ=255)
                                                                                               cGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                            >  2:1386509/1‑111 (MQ=25)
                                                                                                   gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                 <  2:1745932/150‑1 (MQ=255)
                                                                                                    aCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCg                                                <  1:719133/150‑1 (MQ=255)
                                                                                                      ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                              >  1:88578/1‑150 (MQ=255)
                                                                                                          aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                          <  1:1817446/150‑1 (MQ=40)
                                                                                                           aGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                         >  2:1223025/1‑150 (MQ=40)
                                                                                                            gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                        <  2:879677/150‑1 (MQ=39)
                                                                                                            gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                        <  2:2027316/150‑1 (MQ=39)
                                                                                                            gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGGCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                        <  2:1718980/150‑1 (MQ=255)
                                                                                                            gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATGCCGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                        <  2:1856211/150‑1 (MQ=255)
                                                                                                             aTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                       <  1:1801821/150‑1 (MQ=39)
                                                                                                                   ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGc                                                                >  1:1631706/1‑119 (MQ=12)
                                                                                                                   ccGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCAATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                 <  2:1582491/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCACCGGTAATAACGCGAAAGTGATTCAGGACGGCGAGCTGGAGGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  2:599408/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATACCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  2:74594/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGAGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  1:1042838/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  1:666774/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  1:597957/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  1:33815/150‑1 (MQ=255)
                                                                                                                    cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAGTGATATCACCGCTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATggg                                <  1:103696/150‑1 (MQ=255)
                                                                                                                        aCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGAGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATc                                             >  2:1372024/1‑133 (MQ=17)
                                                                                                                           gAAATCAACGGTAATTACGGGAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAg                         <  1:1306251/150‑1 (MQ=255)
                                                                                                                           gAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                                                                                     <  1:1604614/90‑1 (MQ=18)
                                                                                                                               tCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                     <  1:1439866/150‑1 (MQ=255)
                                                                                                                                cAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                      <  2:1958969/68‑1 (MQ=40)
                                                                                                                                cAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                      >  1:1958969/1‑68 (MQ=40)
                                                                                                                                 aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                                 >  1:340020/1‑40 (MQ=255)
                                                                                                                                 aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                   <  2:1659333/150‑1 (MQ=255)
                                                                                                                                  aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCggcg                                                                                          <  1:1038319/78‑3 (MQ=25)
                                                                                                                                  aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCggcg                                                                                          >  2:1038319/1‑76 (MQ=25)
                                                                                                                                  aCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACg                  <  2:2094493/150‑1 (MQ=255)
                                                                                                                                      taataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat                                                                                                    >  1:996381/1‑64 (MQ=37)
                                                                                                                                        ataaCGGGAAAGTGATGCAGGACGGCGATCTGGACGTCAGCGGCGGCGGTCACGGTAGTGATCTCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa            <  1:1256877/150‑1 (MQ=255)
                                                                                                                                         taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGt                                                                             <  2:1717196/84‑1 (MQ=18)
                                                                                                                                         taaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGt                                                                             >  1:1717196/1‑84 (MQ=18)
                                                                                                                                           aCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAgg         >  2:1053429/1‑150 (MQ=255)
                                                                                                                                              ggAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCa      >  1:1497026/1‑150 (MQ=255)
                                                                                                                                              ggAAAGAGATTCAGGACGGCGATCAGGAACTCAGCCGCGGCGGTCACGGTATAGAAATCACCGGTGACAGCGCGACGGTGGATAACAAGGTCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCa      <  1:1303168/150‑1 (MQ=14)
                                                                                                                                               gAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACTGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAt     <  1:1483066/150‑1 (MQ=255)
                                                                                                                                                  aGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGt  >  2:1941259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT  >  NC_000913/1465995‑1466290

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