Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,630 G→C 10.8% pseudogene (1239/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,6300GC10.8% 388.5 / 35.2 128pseudogene (1239/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (68/41);  new base C (8/6);  total (80/48)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466484‑1466779
                                                                                                                                                  |                                                                                                                                                     
gATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAg                                                                                                                                                    >  1:314298/1‑150 (MQ=255)
 atatCACCGGTGACAGCGCTACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                   <  2:2449151/150‑1 (MQ=255)
 atatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                   <  2:973034/150‑1 (MQ=255)
 atatCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACAAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                   <  2:3196229/150‑1 (MQ=255)
 acaacACCGGCGCCAGCGCAACGGTGGATAACAATGGCACCATGACCGTCACCGATACGGAGTCCATCTGTATCCAGGTTGACGTCGACCATGCGTTTTTTAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGa                                                                                                                                                   <  2:2076609/146‑1 (MQ=17)
       ccGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACg                                                                                                                                             >  1:3422611/1‑150 (MQ=255)
        cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                                                                                            <  2:2152140/150‑1 (MQ=255)
           cGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACTCAGATCAACGgaga                                                                                                                                         >  2:1151295/1‑150 (MQ=255)
                 cgcAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGCTAGTTAACAACGAAGGCGAGAGCGCAATCACCAACGGAGGCACCGGCACTCAGATCAACGGTCAAGACGc                                                                                                                                   >  1:3550283/1‑150 (MQ=14)
                          ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                                          <  2:2600807/150‑1 (MQ=255)
                           gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                                                                                         >  1:520496/1‑150 (MQ=255)
                                  aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACg                                                                                                                  <  2:321988/150‑1 (MQ=255)
                                  aGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACg                                                                                                                  <  2:314298/150‑1 (MQ=255)
                                       aCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                             <  2:1526781/150‑1 (MQ=255)
                                             aCCGTCACCGATCCGGAGTCCATGGGTATCCAGCTCGACGGCGACCAGGCGCTCGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                       <  2:2156259/150‑1 (MQ=11)
                                                 tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                              <  2:318136/107‑1 (MQ=255)
                                                 tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAAc                                                                                                                                              >  1:318136/1‑107 (MQ=255)
                                                 tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTAGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGCGGCACCGGCACGCAGATAAACGGTGAAGACGCCACGGCCAAGAACACCGGCAAAACCACCGATg                                                                                                   >  2:1192318/1‑150 (MQ=18)
                                                  cACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                  <  2:3876111/150‑1 (MQ=255)
                                                   aCCGATCCGTAGTCCAACGGTATCCAGGTTGACGGCGAACAGGCGGTTGTTAACAACGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACTCAGATAAACGGTGATGACGCCACCGCAAAAAACAAAGGAAAAACCACAGTTGaa                                                                                                 >  1:2365692/1‑149 (MQ=18)
                                                    ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACg                                                                                                >  1:1589404/1‑150 (MQ=255)
                                                    ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGGGGCACCGGCACGCGGAGTACCCGGGAAGTCGCAACGGCAACCAACGAAGGAAAATCAACCGTTGACg                                                                                                >  1:3393042/1‑150 (MQ=18)
                                                     cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                               >  1:3015571/1‑150 (MQ=255)
                                                     cGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                               >  1:1912501/1‑150 (MQ=255)
                                                        tCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGcca                                                                                            >  1:1398318/1‑148 (MQ=255)
                                                            gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCCAGGAt                                                                                        >  1:3683380/1‑150 (MQ=255)
                                                            gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                        <  2:3233068/150‑1 (MQ=255)
                                                                 cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                   >  1:2605688/1‑150 (MQ=255)
                                                                 cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGAGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                   >  1:2602176/1‑150 (MQ=255)
                                                                  aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAcc                                                                                  >  2:2383812/1‑150 (MQ=255)
                                                                    cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCGAAAAACAACGGCAAAACCACCGTTGAAGGCAAGGATTCCACCgg                                                                                >  1:1920011/1‑150 (MQ=255)
                                                                    cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                                >  2:2070282/1‑150 (MQ=255)
                                                                          ccAGGTTGACAGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACGGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGa                                                                          >  2:3725077/1‑150 (MQ=255)
                                                                            aGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa                                                                        >  2:2877179/1‑150 (MQ=255)
                                                                              gTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt                                                                      >  2:3851198/1‑150 (MQ=255)
                                                                                  aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAATGATTCCACCGGTACGGAAATAGCTg                                                                  >  1:506810/1‑150 (MQ=255)
                                                                                  aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  >  1:3869970/1‑150 (MQ=255)
                                                                                  aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTg                                                                  >  1:1218750/1‑150 (MQ=255)
                                                                                  aCGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGAAAGGATTCCACCGGTACGGAAATTGCTg                                                                  >  1:3869966/1‑150 (MQ=255)
                                                                                   gggCGACCAGGCTTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTgg                                                                 <  2:3641832/149‑1 (MQ=255)
                                                                                    ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCTGGc                                                                >  1:1193441/1‑150 (MQ=255)
                                                                                    ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                                                >  1:2023420/1‑150 (MQ=255)
                                                                                    ggCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGc                                                                >  1:2121275/1‑150 (MQ=255)
                                                                                         ccAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATaa                                                           >  1:1580479/1‑150 (MQ=255)
                                                                                          cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                          <  1:2617970/150‑1 (MQ=255)
                                                                                           cggCGGTTGTTAACACCGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                         <  2:2658136/149‑1 (MQ=255)
                                                                                           aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                         <  2:671491/150‑1 (MQ=255)
                                                                                           aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACg                                                         >  2:1867711/1‑150 (MQ=255)
                                                                                             gCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACggg                                                       >  1:1393175/1‑150 (MQ=255)
                                                                                               ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                     <  2:2519734/150‑1 (MQ=255)
                                                                                                   gtcaataacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                                                 <  1:3593166/144‑1 (MQ=11)
                                                                                                    atAACAACCAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACGGCAAACAACAACGGCAAAACCANCGATGACGGAAAGGAATCCACCGGTACGGAAAAAGCTGGCAATAACGGGAAGGTGa                                                >  2:3194990/2‑150 (MQ=17)
                                                                                                     taagaccGAAGGCGAGAGCCCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACACCAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                               <  2:2365393/144‑1 (MQ=255)
                                                                                                      aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCACGAACGGGACGGTGagt                                              >  1:3490060/1‑148 (MQ=255)
                                                                                                      aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                              >  1:965495/1‑150 (MQ=255)
                                                                                                      aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                              >  1:3232675/1‑150 (MQ=255)
                                                                                                      aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAtt                                              <  1:2504377/150‑1 (MQ=255)
                                                                                                        caacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCa                                            >  2:2929711/1‑150 (MQ=255)
                                                                                                         aacGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAg                                           >  1:3900977/1‑150 (MQ=255)
                                                                                                          acGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAgg                                          >  2:1733728/1‑150 (MQ=255)
                                                                                                           cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                         <  2:3015571/150‑1 (MQ=255)
                                                                                                            gAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                        <  2:670732/134‑1 (MQ=255)
                                                                                                             aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                         >  1:805173/1‑148 (MQ=255)
                                                                                                             aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                         <  2:805173/148‑1 (MQ=255)
                                                                                                                gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAATCAACAACGGCAAATCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                    >  1:2361269/1‑150 (MQ=255)
                                                                                                                gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                    >  1:2543789/1‑150 (MQ=255)
                                                                                                                    gagCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCTGGCACGAACGGGACGGGGATTCAGGACGGCgcgcg                                >  2:466694/1‑146 (MQ=17)
                                                                                                                      gcaccATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTgg                              <  2:3900977/147‑1 (MQ=255)
                                                                                                                       cgcCATCACCAACGGTGGCACCGGCACTCAGATTAACGGTGATGACGCCACGGCAAACATCAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAATGTGATTCAGGACGGCGATCTGGa                             <  1:2394171/150‑1 (MQ=255)
                                                                                                                       cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  2:3422611/150‑1 (MQ=255)
                                                                                                                       cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  1:1827084/150‑1 (MQ=255)
                                                                                                                       cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGa                             <  1:2394143/150‑1 (MQ=255)
                                                                                                                        accATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                   <  2:2693737/110‑1 (MQ=17)
                                                                                                                        accATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                   >  1:2693737/2‑111 (MQ=17)
                                                                                                                          cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  1:1219838/1‑150 (MQ=255)
                                                                                                                          cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                          >  2:1765196/1‑150 (MQ=255)
                                                                                                                            tCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                        >  1:1726410/1‑150 (MQ=255)
                                                                                                                              aCCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgc                      >  1:706518/1‑150 (MQ=255)
                                                                                                                              aCCAACGGTGGCACCGGCACGCAGATTAAAGGAGATGACGCCACGGCAAACAACAACGGCAAAACCACCGAAGACGGCAAGGATTCCACCGGGACGGAAATTGCTGGCAATAACGGGAAGTAGATTCAGGACGGCGATCTGGATGTCAgc                      >  2:2904565/1‑150 (MQ=17)
                                                                                                                               ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                     >  2:3939820/1‑150 (MQ=255)
                                                                                                                               ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcg                     >  2:3745702/1‑150 (MQ=255)
                                                                                                                               ccAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACAAAGGCAAAACCACAGTTGAAGGCAAGGATTCCACCGGTACGCAAATAGCTGGCAATAACGGGAAGGTCAATCAGGAACGCGATCTCGAAGTCAccg                     >  2:3745230/1‑150 (MQ=17)
                                                                                                                               ccAACGGTGGCACCGGCACGCAGAGTAACGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCGGGCACTCACGGGACGGGGATTCAGGACGGCGATCTGGAGGTCAgcg                     >  2:830699/1‑150 (MQ=17)
                                                                                                                                cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                    >  1:165974/1‑150 (MQ=255)
                                                                                                                                cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                    >  1:3172937/1‑150 (MQ=255)
                                                                                                                                 aaCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggc                   >  1:633436/1‑150 (MQ=255)
                                                                                                                                   cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCATGGATTCCACCGGTACGGAAATTGCTGGAAATAACGGGATGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                 <  1:1973736/150‑1 (MQ=255)
                                                                                                                                   cGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                 <  1:2929711/150‑1 (MQ=255)
                                                                                                                                     gTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg               >  1:1455587/1‑150 (MQ=255)
                                                                                                                                     gTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg               >  1:1029804/1‑150 (MQ=255)
                                                                                                                                        gCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTc            >  1:1729259/1‑150 (MQ=255)
                                                                                                                                         cACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           >  1:2082185/1‑150 (MQ=255)
                                                                                                                                         cACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           <  1:1733728/150‑1 (MQ=255)
                                                                                                                                          aCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           >  2:3471976/1‑149 (MQ=255)
                                                                                                                                          aCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa           <  1:3471976/149‑1 (MQ=255)
                                                                                                                                            cGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACgg        >  1:2559971/1‑150 (MQ=255)
                                                                                                                                             ggCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAATTTGCTGGAAATAACGGTAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGt       <  1:971340/150‑1 (MQ=255)
                                                                                                                                                aCGCAGCTTAACGGTGATGCCGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGAGTCCACCGGTACGGAAATGGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt    <  2:2426862/150‑1 (MQ=17)
                                                                                                                                                  gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  2:2700868/1‑150 (MQ=255)
                                                                                                                                                  gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa  >  2:2741442/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGA  >  NC_000913/1466484‑1466779

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: