Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,639 C→T 6.5% pseudogene (1248/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,6390CT 6.5% 867.2 / 30.7 284pseudogene (1248/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (155/106);  new base T (10/8);  total (168/116)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

ACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466490‑1466788
                                                                                                                                                     |                                                                                                                                                     
     cGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtga                                                                                                                                                  <  2:3869630/150‑1 (MQ=255)
       aCAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatg                                                                                                                                                >  2:6577704/1‑150 (MQ=255)
        cAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGtgatga                                                                                                                                               >  2:6394072/1‑150 (MQ=255)
         aGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGAc                                                                                                                                              <  2:6454113/150‑1 (MQ=255)
                cGGTGGATAACAAGGGATCCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACTGCGACCAGGCGGTTGTTAAGAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACgg                                                                                                                                       <  2:3095327/150‑1 (MQ=17)
                 ggTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTCGACGGCGACCAGGCGGTCGTGAAGAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGc                                                                                                                                      >  1:2606032/1‑150 (MQ=14)
                    ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCaaa                                                                                                                                   <  2:7574870/150‑1 (MQ=255)
                          cAAGGGCACCATGACCGCCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaa                                                                                                                             >  2:4956866/1‑150 (MQ=255)
                                cACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaa                                                                                                                       >  2:6152721/1‑150 (MQ=255)
                                 aCCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGGTTGACGGCGACCAGGCGTTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaa                                                                                                                      <  2:8157329/150‑1 (MQ=255)
                                  ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                                     <  2:5301306/150‑1 (MQ=255)
                                  ccATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCaaaa                                                                                                                     <  2:5606582/150‑1 (MQ=255)
                                   cATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                    >  1:1744417/1‑150 (MQ=255)
                                       aCCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                                >  2:7357931/1‑150 (MQ=255)
                                         cGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                                                                                                              >  1:8183357/1‑150 (MQ=255)
                                           tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTg                                                                                                            >  1:1038448/1‑150 (MQ=255)
                                            cACCGATCCGGTGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCGCGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGa                                                                                                           <  1:6194108/150‑1 (MQ=255)
                                             aCCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAc                                                                                                          >  1:3876493/1‑150 (MQ=255)
                                              ccGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAAAAACATCGGCCAAACCACCGTTGACg                                                                                                         >  2:4427532/1‑150 (MQ=255)
                                               cGATCCGGAGTCCATGGGTATCCAGGTTGACGGGGACAAGGCGTTTGTTAAAAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAAAGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                        >  1:2209984/1‑150 (MQ=14)
                                                 aTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCa                                                                                                      >  1:6588538/1‑150 (MQ=255)
                                                   ccGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAg                                                                                                    <  1:5205674/150‑1 (MQ=255)
                                                        gTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTc                                                                                               <  2:1636596/150‑1 (MQ=255)
                                                            aTCGGTATCCAGGTTGACGGCGACCAGGCGGTGGTGAACAACGAAGGCGAGAGCGCGATCACCAAAGGAGGCACCGGCACTCAGATAAACGGTGAAGACGCCACCGCGAAGAACAAGGGAAAAACCACCGTTGAAGGAAAGGATTCCAcc                                                                                           >  1:7089551/1‑150 (MQ=255)
                                                              cGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCgg                                                                                         >  1:4374915/1‑150 (MQ=255)
                                                                gTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGtc                                                                                       >  1:6823406/1‑149 (MQ=255)
                                                                   tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                    >  2:6928318/1‑150 (MQ=255)
                                                                   tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACgg                                                                                    >  2:898701/1‑150 (MQ=255)
                                                                     cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  >  1:6665078/1‑150 (MQ=255)
                                                                     cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGAAGACGCCACGGCCAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaa                                                                                  >  2:7128537/1‑150 (MQ=255)
                                                                     cAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAAAGGAGGCACCGGCACTCAGATGAACGGTGACGACGCCACCGCGAAAAACAGAGGAAAAACCACCGTTGAAGGCAAGGATTCCACCGGTACCGaa                                                                                  >  2:1305743/1‑150 (MQ=255)
                                                                       ggTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCTCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                <  2:2606032/150‑1 (MQ=255)
                                                                          tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGc                                                                             <  2:1107778/150‑1 (MQ=255)
                                                                          tGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGc                                                                             >  2:5233867/1‑150 (MQ=255)
                                                                           gACGGTGACAAGGCCATCGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                            <  1:3229824/150‑1 (MQ=11)
                                                                           gACGGCGACCAGGCGGTTGTTAACAACGTAGGCGAGTGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                            <  1:5865317/150‑1 (MQ=255)
                                                                             cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTgg                                                                          >  1:3076758/1‑150 (MQ=255)
                                                                                    cAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAAc                                                                   <  2:4026716/150‑1 (MQ=255)
                                                                                     aGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACATTGCTGGCAATAACg                                                                  >  2:2308193/1‑150 (MQ=255)
                                                                                      ggCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACgg                                                                 <  2:4217590/150‑1 (MQ=255)
                                                                                       gCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACggg                                                                >  1:4834809/1‑150 (MQ=255)
                                                                                         ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                              <  2:6613852/150‑1 (MQ=255)
                                                                                         ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                              >  1:3567296/1‑150 (MQ=255)
                                                                                         ggTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGaa                                                              <  1:8251309/150‑1 (MQ=255)
                                                                                          ttagtaaacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAg                                                             <  2:5066601/144‑1 (MQ=255)
                                                                                            tgttAACAACGAAGGCGAGAGCGCAATCACCAACGGTGGCACCGGCACGCAGATAAACGGTGAAGACGCCACCGCCAAAAACAAAGGCAAAACCACAGTTGAAGGGAACGATACAGCAGGAACCGTAAATCCAGGCACTAACGCGAAGGt                                                           >  1:5621804/1‑150 (MQ=9)
                                                                                             gttAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTg                                                          <  2:3189719/150‑1 (MQ=255)
                                                                                               tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCGAACAACAAAGGCAAAACCACCGTTGAAGGCAATGATTCCACCGGTACCGAAATAGCTGGCAATAACGGGAAGGTGAt                                                        >  1:5003629/1‑150 (MQ=17)
                                                                                               tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                        >  1:4890372/1‑150 (MQ=255)
                                                                                               tAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                        <  1:4458496/150‑1 (MQ=255)
                                                                                               tAACAACGAAGGCGAGAGCGCCAGCACCTACTGTGGGTCCGTCACGCAGTTGAACGGTGTTGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGAt                                                        <  1:4087728/150‑1 (MQ=18)
                                                                                                aacaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGCAGGTGAtt                                                       >  2:5816293/1‑150 (MQ=255)
                                                                                                 acaacGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTc                                                      >  2:2798973/1‑150 (MQ=255)
                                                                                                     cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGa                                                  >  1:6862659/1‑150 (MQ=255)
                                                                                                      gAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACACAGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGAc                                                 >  1:3300849/1‑150 (MQ=255)
                                                                                                       aaGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACg                                                <  1:5803124/150‑1 (MQ=255)
                                                                                                       aaGGCGAGAGCACCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGTAGGTGATTCAGGACg                                                <  1:5476035/150‑1 (MQ=255)
                                                                                                        aGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTCCGGAAATTGCTGGCAATAACGGGACGGTGATTCAGGACgg                                               >  1:85921/1‑150 (MQ=255)
                                                                                                        aGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACgg                                               >  1:7435081/1‑150 (MQ=255)
                                                                                                          gCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCg                                             >  2:7483474/1‑150 (MQ=255)
                                                                                                            gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCACGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                           >  2:5527770/1‑150 (MQ=255)
                                                                                                            gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                           >  2:6710769/1‑150 (MQ=255)
                                                                                                            gagaGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGAt                                           >  1:795619/1‑150 (MQ=255)
                                                                                                              gagCGCCATCACCAACGGTGGCACCGGCACTCAGATAAACGGTAAAGACGCCACGGCAAAAAACAAAGGAAAAACCACAGTTGACGGCAAGGATTCCACCGGTACGAAAATAGCGGGCAATATCGGCAAGGTGATTCAGGACGGCGATCt                                         >  1:1579345/1‑150 (MQ=14)
                                                                                                               agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg                                        <  1:6577704/150‑1 (MQ=255)
                                                                                                               agCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTg                                        <  1:6018946/150‑1 (MQ=255)
                                                                                                               agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGTATTCCACAGGTACGGAAATAGATGGGAATAACGGGAAGGTGATTCAGGACGGCGATCTg                                        >  2:7843473/1‑150 (MQ=11)
                                                                                                                tatccATCACCAACGGTGGCACCGGCACGCAGTTTACTGTTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCTGTACGGAAATTGCTGGATATAACGGGAAGGTGATTCAGGACGGCGATCTgg                                       <  1:7843473/147‑1 (MQ=14)
                                                                                                                  gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  1:120164/150‑1 (MQ=255)
                                                                                                                  gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  2:8183357/150‑1 (MQ=255)
                                                                                                                  gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  2:1241657/150‑1 (MQ=255)
                                                                                                                  gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  1:7357931/150‑1 (MQ=255)
                                                                                                                  gcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGAt                                     <  1:649697/150‑1 (MQ=255)
                                                                                                                  cccAACACCACCGGTGCCACCGGCACTCAGATTAACGGTGATCACCCCACGTCAAACAACAACGGCAAAACCACCGCAGACGCCAACGATACCACGCGTACGGAAATAACTGGCAATAACGGTAATGTTATTCAGGACGGCGATCTGGAt                                     <  2:901613/149‑1 (MQ=2)
                                                                                                                    cATCACCAACGGTGGCACCGGCACGCAGATTAAGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTgg                                                                          >  2:3566831/1‑111 (MQ=255)
                                                                                                                    cATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGt                                   >  2:1918406/1‑150 (MQ=255)
                                                                                                                      tCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCAc                                                                                                                                         <  1:4740026/46‑1 (MQ=255)
                                                                                                                      tCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCa                                 >  1:6796125/1‑150 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaacaacaac                                                                                                                            >  1:899846/1‑55 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                            <  2:3971923/103‑1 (MQ=40)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                             >  1:1958490/1‑150 (MQ=255)
                                                                                                                          cAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcgg                             >  2:2508207/1‑150 (MQ=255)
                                                                                                                               gTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                        <  1:1067844/150‑1 (MQ=255)
                                                                                                                                 ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                          >  2:4285008/1‑146 (MQ=255)
                                                                                                                                 ggcaccggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                          <  1:4285008/146‑1 (MQ=255)
                                                                                                                                 ggCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                      <  2:1999476/150‑1 (MQ=255)
                                                                                                                                    aCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                   >  1:4052958/1‑150 (MQ=255)
                                                                                                                                    aCCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                   >  2:7892016/1‑150 (MQ=255)
                                                                                                                                         cACGCAGATTCACGGTGAGGACGCCTCGGCAAACACCAACGGCAAAACCACCGGTGACGGCAAGGATTCCCCCGGTACGGAAATGGCGGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt              <  2:6240259/150‑1 (MQ=17)
                                                                                                                                          acccAGATTAAGGGTGATGACGCCACGGCAAACAACACCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATACCGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt             <  1:3030144/147‑1 (MQ=255)
                                                                                                                                            gaaGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa           <  1:2916348/148‑1 (MQ=255)
                                                                                                                                            gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTCTTGa           >  2:8071205/1‑150 (MQ=255)
                                                                                                                                            gCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa           <  1:3782764/150‑1 (MQ=255)
                                                                                                                                            ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg            <  2:6367012/148‑1 (MQ=255)
                                                                                                                                            ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg            <  2:466014/148‑1 (MQ=255)
                                                                                                                                            ccAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAGCGGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg            >  1:6367012/2‑149 (MQ=255)
                                                                                                                                             cAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat          >  2:4108831/1‑150 (MQ=255)
                                                                                                                                             cAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCt                                                                            <  1:8119279/84‑1 (MQ=39)
                                                                                                                                                   aaCGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc    >  1:2673628/1‑150 (MQ=255)
                                                                                                                                                     cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTCTTGATATCACCgg  >  1:1243403/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGG  >  NC_000913/1466490‑1466788

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: