Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,856 A→G 5.1% pseudogene (1465/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,8560AG 5.1% 318.9 / 2.5 101pseudogene (1465/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (58/32);  new base G (3/2);  total (63/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

TGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAAT  >  NC_000913/1466711‑1467004
                                                                                                                                                 |                                                                                                                                                    
tGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGa                                                                                                                                                  >  1:460129/1‑150 (MQ=255)
tGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGa                                                                                                                                                  >  1:460109/1‑150 (MQ=255)
   tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                               >  2:1880782/1‑150 (MQ=255)
   tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                               >  2:593786/1‑150 (MQ=255)
   tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                               >  2:593838/1‑150 (MQ=255)
   tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                               >  1:630112/1‑150 (MQ=255)
   tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                                                                                               >  2:1639859/1‑150 (MQ=255)
           aCGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                                                                                                       >  1:41594/1‑150 (MQ=39)
                  ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGt                                                                                                                                <  1:2204998/150‑1 (MQ=40)
                  ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGt                                                                                                                                <  1:2394718/150‑1 (MQ=40)
                   gggATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGACCATCGGTCTCCAGATTGACGGCGACCAGGCCATCGTc                                                                                                                               >  2:2024523/3‑150 (MQ=255)
                   gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                                                                                               >  1:1448528/1‑150 (MQ=40)
                   gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                                                                                               >  2:1465708/1‑150 (MQ=40)
                        tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                                                                                          >  1:1420594/1‑150 (MQ=255)
                        tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                                                                                          >  1:1401855/1‑150 (MQ=255)
                        tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                                                                                          >  1:874520/1‑150 (MQ=255)
                          aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCCTCGTCAAAAACg                                                                                                                        >  1:2297470/1‑150 (MQ=255)
                           ggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGa                                                                                                                       >  1:1546073/1‑150 (MQ=255)
                            gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                                                                                                      >  2:854977/1‑150 (MQ=255)
                             aCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAg                                                                                                                     >  2:2174840/1‑150 (MQ=255)
                                 cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                                                                                                 >  2:979538/1‑150 (MQ=255)
                                 cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                                                                                                 >  1:626302/1‑150 (MQ=255)
                                 cGATCAGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCCGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                                                                                                 >  1:626239/1‑150 (MQ=255)
                                  gATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgag                                                                                                                >  2:185438/1‑150 (MQ=255)
                                   aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgaga                                                                                                               >  1:1192233/1‑150 (MQ=255)
                                   aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgaga                                                                                                               >  1:2299725/1‑150 (MQ=255)
                                    tCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGagag                                                                                                              >  1:496217/1‑150 (MQ=255)
                                    tCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGagag                                                                                                              <  2:1887684/150‑1 (MQ=255)
                                      tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgaga                                                                                                               >  1:889880/1‑147 (MQ=255)
                                      tGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCgaga                                                                                                               <  2:889880/147‑1 (MQ=255)
                                       ggATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                                                                                           >  2:2380262/1‑150 (MQ=255)
                                            tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc                                                                                                                                                 <  2:1908078/107‑1 (MQ=32)
                                            tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc                                                                                                                                                 >  1:1908078/1‑107 (MQ=32)
                                            tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                                                      <  1:1418115/150‑1 (MQ=255)
                                            tCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCa                                                                                                      <  1:185769/150‑1 (MQ=255)
                                             cAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                                                                                                     >  2:613549/1‑150 (MQ=255)
                                                    ggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggc                                                                                              <  2:630112/150‑1 (MQ=255)
                                                        gTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                                          >  2:1874390/1‑150 (MQ=255)
                                                        gTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggca                                                                                          >  2:1874358/1‑150 (MQ=255)
                                                         tCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                                         >  1:203727/1‑150 (MQ=255)
                                                         tCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                                         >  1:526015/1‑150 (MQ=255)
                                                         tCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                                         <  2:1646403/150‑1 (MQ=255)
                                                         tCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                                         >  1:19512/1‑150 (MQ=255)
                                                             ggTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     <  1:2004212/150‑1 (MQ=255)
                                                             ggTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     <  2:1363379/150‑1 (MQ=255)
                                                              gTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggca                                                                                    <  2:2297470/150‑1 (MQ=255)
                                                               tATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                                                   >  2:1241885/1‑150 (MQ=255)
                                                               tATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggcac                                                                                   >  2:1241822/1‑150 (MQ=255)
                                                                       tCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCa                                                                           >  1:345596/1‑150 (MQ=255)
                                                                        cACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACGATCACCAATGGCGGCACCGGCACGCAGATCaa                                                                          >  2:166388/1‑150 (MQ=255)
                                                                         aCCGGCGACAGCGCAACGGTGGATGACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                                                         <  2:19512/150‑1 (MQ=255)
                                                                         aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAAc                                                                         <  2:203727/150‑1 (MQ=255)
                                                                           cGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACgg                                                                       >  1:1118785/1‑150 (MQ=255)
                                                                            ggCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCCTAGTAAACAACGAAGGCGAGAGCACGATCACCAAAGGCGGCACCGGCACTCAGATAAACGGt                                                                      >  2:1863716/1‑150 (MQ=11)
                                                                                       aaCGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCCTCGTCAAGAACGAAGGCGAGAGCCCAATCACCAAAGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAc                                                           >  2:508758/1‑150 (MQ=17)
                                                                                          ggTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACcgc                                                        <  2:1546073/150‑1 (MQ=255)
                                                                                              gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTACCGACGCCACCGCGAAt                                                    >  1:1506408/1‑150 (MQ=255)
                                                                                              gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAAt                                                    >  2:791384/1‑150 (MQ=255)
                                                                                                tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCTGATTGACGGCGACCAGGCCATCGTCAATAACGATGGCGAGAGCACTATCACCAATGGTGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATaa                                                  <  1:284235/150‑1 (MQ=255)
                                                                                                tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCATTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATaa                                                  <  2:2323323/150‑1 (MQ=255)
                                                                                                 aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATgg                                                                                               >  2:1123894/1‑104 (MQ=255)
                                                                                                 aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATgg                                                                                               <  1:1123894/104‑1 (MQ=255)
                                                                                                 aaCAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAAc                                                 >  1:1168034/1‑150 (MQ=255)
                                                                                                    aaGGGCACCAAGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGAGACCAGGCGGAAGTCAATAAAGAAGGCGAGAGCACCATCACCAAAGGAGGCACCGGCACTCAGAAAAACCCTAAAGACGCCACCGCGAAGAACaca                                              >  1:711394/1‑148 (MQ=14)
                                                                                                       ggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCTGTTTGACGGCGACCTGGCGTTTGTAAAAAACGAAGGCGAGAGCTCTATCACCAAGGGGGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                                           <  1:508742/150‑1 (MQ=17)
                                                                                                       ggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCCGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCATTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGa                                           <  1:508758/150‑1 (MQ=255)
                                                                                                        gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     >  2:316736/1‑107 (MQ=255)
                                                                                                        gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     <  1:316736/107‑1 (MQ=255)
                                                                                                        gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     >  2:1967477/1‑107 (MQ=255)
                                                                                                        gCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc                                                                                     <  1:1967477/107‑1 (MQ=255)
                                                                                                           ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGaaaaa                                       >  1:826468/1‑150 (MQ=255)
                                                                                                              tGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCa                                    >  1:2067194/1‑150 (MQ=255)
                                                                                                                aCCGTCACCGATCGGGTGTCCATCGTTATCCTGTTTGACGGCGACCAGGCCTTCGTAAATAACGACGGCGAGAGCACTATCACCATTGGCGGCCCCGGCACTCAGATCACCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACt                                  <  1:1863716/150‑1 (MQ=255)
                                                                                                                    tCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTg                              <  2:874568/150‑1 (MQ=255)
                                                                                                                    tCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTg                              >  2:133402/1‑150 (MQ=255)
                                                                                                                    tCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAGAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTg                              <  2:626302/150‑1 (MQ=255)
                                                                                                                    acACCGATCCGGAGTCCAACGGTATCCAGATTGACGGCGACCAGGCCATCGTAAAAAACGAAGGCGAGACCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTg                              <  2:1015650/149‑1 (MQ=255)
                                                                                                                      aCCGATCCGGAGTCCATCGGTATCCAGATTGCGGGCGACCAGGCCATCGTCAATAACGACGGCGAGAGCACTATCACCAATGGCGGCCCCGGCACTCAGATCACCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGAt                            <  1:1651246/150‑1 (MQ=25)
                                                                                                                      aCCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCa                                                                                                                                         <  2:2330250/41‑1 (MQ=38)
                                                                                                                         gATCCGGAGTCCATCGGTATCCAGTTTGTCGGCTACCAGGCGATCGTATATAACGAAGGCGAGAGCACTATCACCATTGGCGGCACCGGCACTCAGATCTACGGTAACGACGCCACCGCGAATTACAGTGTAAAAACCACTGTTGTTggt                         <  2:711394/150‑2 (MQ=21)
                                                                                                                         gATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCCGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGa                         <  2:711393/150‑1 (MQ=255)
                                                                                                                         gATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGa                         <  2:2166958/150‑1 (MQ=255)
                                                                                                                          aTCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGaa                        >  1:1422868/1‑150 (MQ=255)
                                                                                                                          aTCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCAATCTTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGaa                        <  1:978905/150‑1 (MQ=255)
                                                                                                                             cGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAg                     >  2:1623232/1‑150 (MQ=255)
                                                                                                                               gAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGAt                   >  2:1787698/1‑150 (MQ=255)
                                                                                                                               gAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGAt                   <  2:345596/150‑1 (MQ=255)
                                                                                                                               gAGTCCATCGGTATCCAGAGTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGAt                   <  2:2366145/150‑1 (MQ=255)
                                                                                                                                 gTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATACCAGTGGAAAAACCACTGTTGATGGAAAAGATTc                 <  2:526015/150‑1 (MQ=255)
                                                                                                                                   ccATCGGTATCCAGATTGACGGCGACCAGGCCCTCGTAAATAACGAAGGCGAGAGCCCTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCa               <  2:966492/150‑1 (MQ=255)
                                                                                                                                     aTCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACg             >  1:1696626/1‑150 (MQ=255)
                                                                                                                                     aTCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACg             >  1:1999078/1‑150 (MQ=255)
                                                                                                                                       cGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAATACCACTGTTGATGGAAAAGATTCCACggg           <  2:2061944/150‑1 (MQ=255)
                                                                                                                                       cGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACggg           >  2:83689/1‑150 (MQ=255)
                                                                                                                                         gTATCCAGGTTGACGGCGACCAGGCGTTCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTa         <  2:333214/150‑1 (MQ=255)
                                                                                                                                           aTCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTAcc       >  2:1851663/1‑150 (MQ=255)
                                                                                                                                            tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCa      >  1:2109724/1‑150 (MQ=255)
                                                                                                                                            tCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCa      >  1:1713450/1‑150 (MQ=255)
                                                                                                                                                gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAAt  >  2:1890230/1‑150 (MQ=255)
                                                                                                                                                gATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAAt  >  1:184361/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAAT  >  NC_000913/1466711‑1467004

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: