Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,604 G→A 7.9% pseudogene (1213/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,6040GA 7.9% 170.0 / 5.2 66pseudogene (1213/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (38/20);  new base A (3/2);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466457‑1466753
                                                                                                                                                   |                                                                                                                                                     
gATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAgcgtc                                                                                                                                                     >  1:2159342/1‑148 (MQ=255)
    tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCa                                                                                                                                                 <  2:1031647/150‑1 (MQ=255)
      aGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                                                                                                                               >  1:2782581/1‑150 (MQ=255)
          gcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                                                                                                                           <  1:39555/150‑1 (MQ=255)
          gcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACg                                                                                                                                           <  1:2906637/150‑1 (MQ=255)
            ggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCACTATCACCAAAggc                                                                                                                                         >  2:2326499/1‑149 (MQ=255)
                gTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCa                                                                                                                                     <  2:839579/150‑1 (MQ=255)
                    cGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCgg                                                                                                                                 <  2:48430/150‑1 (MQ=255)
                      gTATTGATAGCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCa                                                                                                                               <  2:676438/150‑1 (MQ=255)
                              atCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGAtt                                                                                                                       >  2:2683795/1‑150 (MQ=255)
                                  cgggCGACAGCGGAGCGGTGGATATGTAGGGCACCTTGTCGGTCACGGTTGCGGAGTCCATGGGTATCCAGTTTGACGGCGACCAGGCGGTGGTTAAGAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                   <  2:1696739/148‑1 (MQ=17)
                                   cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                                                                  >  2:2318074/1‑150 (MQ=255)
                                    ggCGACAGCGCAACCGTTGATATCAAGGGCACCATGACCGTCTCCGATCCTGTTCGCATCTGTCTACTGTTTGACGGTGACAAGGCCTTAGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGt                                                                                                                 <  2:982283/150‑1 (MQ=2)
                                        aCAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACTCAGATCAACGGtgagg                                                                                                             >  2:1524379/1‑148 (MQ=14)
                                                 cGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGTTTGACGGAGACAAGGCGTTTGTTAATAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACgg                                                                                                    >  1:1100829/1‑150 (MQ=11)
                                                   gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                  >  2:2883763/1‑150 (MQ=255)
                                                   gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCa                                                                                                  >  2:1579581/1‑150 (MQ=255)
                                                        tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAGGACGCCACGGCAaacaa                                                                                             >  1:2325316/1‑150 (MQ=255)
                                                        gcgatctggCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAAGAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaa                                                                                             <  2:2310068/143‑1 (MQ=14)
                                                           cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCGATCACCAAAGGGGGCACCGGCACGCAGATGAACGGTGAGGACGCCACGGCGaagaacag                                                                                          >  1:1636334/1‑149 (MQ=14)
                                                            aaGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAaacaacaac                                                                                         >  1:1692624/1‑150 (MQ=255)
                                                                   ccATGACCCTCACCGATCCGGAGTCCATCGGTATCCAGCTAGACGGCGACCAGGCGGTAGATAACAACGAAGGCGAGAGCGCCAACACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                  <  2:1100829/150‑1 (MQ=17)
                                                                     aTGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAGGGAGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCGAAAAACAGAGGAAAAacc                                                                                >  2:2343124/1‑150 (MQ=14)
                                                                        aaccacACCGATCATCAGACCATCGGTATCCAGCTAGACGGAGACCAGGCCTTTGTAAAAAACGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                             <  2:642647/145‑1 (MQ=14)
                                                                           gTCACCGATCCGTAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGtt                                                                          <  2:1827275/150‑1 (MQ=255)
                                                                           gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAGGACGCCACCGCCAAAAACACCGGCAAAACCACCGtt                                                                          >  1:602579/1‑150 (MQ=17)
                                                                           gTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTGAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCCAAAAACACAGGCAAAACCACAGtt                                                                          >  1:322413/1‑150 (MQ=17)
                                                                            tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCGATCACCAAAGGCGGCACCGGCACGCAGATAAACGGTGAGGACGCCACGGCGAATAACAGAGGAAAAACCACCGTTg                                                                         >  2:2645233/1‑150 (MQ=14)
                                                                            tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                         >  2:1036046/1‑150 (MQ=255)
                                                                             cACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGTGACAAGGCCTTTGTAAAAAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGTTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                        <  1:1185760/150‑1 (MQ=11)
                                                                                cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                     >  2:2657669/1‑150 (MQ=255)
                                                                                cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                     >  2:2657670/1‑150 (MQ=255)
                                                                                 gATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                                                                    >  2:2364066/1‑150 (MQ=255)
                                                                                   tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCaa                                                                  >  1:995459/1‑150 (MQ=255)
                                                                                     cGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAgg                                                                <  2:216096/150‑1 (MQ=255)
                                                                                       gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTCATGACGCCACCGCAAAAAACAAAGGCAAAACCACCGTTGACGGCAAGGAt                                                              >  2:2868328/1‑150 (MQ=18)
                                                                                          tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTcc                                                           >  1:1700385/1‑150 (MQ=255)
                                                                                          tCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCAAAAAACAGAGGCAAAACAACCGCACACGGCAAGCATAcc                                                           >  1:2521866/1‑150 (MQ=18)
                                                                                             aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        >  1:1528133/1‑150 (MQ=255)
                                                                                             aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        <  2:2903502/150‑1 (MQ=255)
                                                                                             aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                        >  1:2669041/1‑150 (MQ=255)
                                                                                              tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                       <  2:2765882/150‑1 (MQ=255)
                                                                                              tCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCg                                                       <  2:2765889/150‑1 (MQ=255)
                                                                                                         gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAAtt                                            >  1:1514794/1‑150 (MQ=255)
                                                                                                         gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACAtt                                            >  2:2589965/1‑150 (MQ=255)
                                                                                                         gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt                                            >  1:2524243/1‑150 (MQ=255)
                                                                                                            gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCt                                         >  1:1833455/1‑150 (MQ=255)
                                                                                                             tcGGCGACCAGGCGGTTGTTAACCACGCAGGCGCGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTACCGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                        <  2:2613710/149‑1 (MQ=17)
                                                                                                             aCGGAGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                        <  2:2088247/150‑1 (MQ=255)
                                                                                                               ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                      >  1:2459433/1‑150 (MQ=255)
                                                                                                               ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAAAGGAAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACCGAAATAGCTGGc                                      >  2:1213190/1‑150 (MQ=17)
                                                                                                                 cGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCaa                                    >  1:1775783/1‑150 (MQ=255)
                                                                                                                  gACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGGGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCGAAAAACATCGGAAAAACCACCGTTGACGGCAAGGAATCCACCGGTACGGAAATAGCTggacat                                   >  2:675019/1‑146 (MQ=17)
                                                                                                                     cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                >  2:1293272/1‑150 (MQ=255)
                                                                                                                     cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCTAAAAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAAc                                >  2:652748/1‑150 (MQ=17)
                                                                                                                                  acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                   <  2:2572853/150‑1 (MQ=255)
                                                                                                                                  acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                   >  1:172616/1‑150 (MQ=255)
                                                                                                                                   taacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                  <  2:1061607/149‑1 (MQ=255)
                                                                                                                                   caacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                  <  2:1061442/150‑1 (MQ=255)
                                                                                                                                        aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg             >  1:1944512/1‑150 (MQ=255)
                                                                                                                                        aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                         <  2:1800022/122‑1 (MQ=25)
                                                                                                                                        aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                         >  1:1800022/1‑122 (MQ=25)
                                                                                                                                           gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                >  2:250271/1‑112 (MQ=25)
                                                                                                                                           gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                <  1:250271/112‑1 (MQ=39)
                                                                                                                                             gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt        >  1:1442281/1‑150 (MQ=255)
                                                                                                                                             gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt        >  1:1071426/1‑150 (MQ=255)
                                                                                                                                               gagCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCGGGCAATAACGGGAAGGTGATTCAGGACGGCGATCt      >  2:2660364/1‑150 (MQ=255)
                                                                                                                                                  cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa   >  1:677/1‑150 (MQ=255)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCTGGCAATAACGGGAAGGTGATTCAGGACGGAGAGCTGGAt  >  1:2286912/1‑150 (MQ=17)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  <  1:2657670/150‑1 (MQ=255)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  <  1:2657669/150‑1 (MQ=255)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  >  1:2152304/1‑150 (MQ=255)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt  >  1:2152370/1‑150 (MQ=255)
                                                                                                                                                   gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAGCTGGAt  >  1:593270/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAT  >  NC_000913/1466457‑1466753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: