Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,187,989 A→G 5.5% pseudogene (558/2434 nt) yqiG → putative outer membrane usher protein YqiG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,187,9890AG 5.5% 505.2 / 12.7 163pseudogene (558/2434 nt)yqiGputative outer membrane usher protein YqiG
Reads supporting (aligned to +/- strand):  ref base A (75/79);  new base G (4/5);  total (79/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 165 total aligned reads displayed.

TGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCT  >  NC_000913/3187840‑3188137
                                                                                                                                                     |                                                                                                                                                    
tGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGa                                                                                                                                                      <  2:904093/150‑1 (MQ=255)
tGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGa                                                                                                                                                      <  2:66656/150‑1 (MQ=255)
 gAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGAt                                                                                                                                                     >  2:2050238/1‑150 (MQ=255)
    cTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                                                                  >  1:137211/1‑150 (MQ=255)
     tttGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAcc                                                                                                                                                 <  2:407424/150‑1 (MQ=255)
         gAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTa                                                                                                                                             >  1:1026042/1‑150 (MQ=255)
         gAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTa                                                                                                                                             >  2:1231472/1‑150 (MQ=255)
          aaTCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTAc                                                                                                                                            <  2:1788755/150‑1 (MQ=255)
            tCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGa                                                                                                                                          >  1:2441196/1‑150 (MQ=255)
              aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                        >  2:70008/1‑150 (MQ=255)
              aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                        <  2:819798/150‑1 (MQ=255)
              aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                        >  1:1659582/1‑150 (MQ=255)
              aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                        <  1:2451907/150‑1 (MQ=255)
               tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCa                                                                                                                                       <  2:591321/150‑1 (MQ=255)
               tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCa                                                                                                                                       >  1:262412/1‑150 (MQ=255)
               tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGGACTGATTACCAGTACGATCa                                                                                                                                       <  2:950237/150‑1 (MQ=255)
                ggTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAg                                                                                                                                      >  1:480604/1‑150 (MQ=255)
                 gTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGt                                                                                                                                     <  1:573178/150‑1 (MQ=255)
                     tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                 >  2:1031082/1‑150 (MQ=255)
                     tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGAAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCATGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                 >  2:719315/1‑150 (MQ=255)
                       cTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGAAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTa                                                                                                                               <  1:459162/150‑1 (MQ=255)
                        tGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTaa                                                                                                                              >  2:1296514/1‑150 (MQ=255)
                        tGGGGCATTTCTTGAATACAATTTATATGCCTCTCATTATGCAACACATAAGGCCGATAATTCGAAAAATATAAGTTCCTATGGGCAGGCTGGGGTTTATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATATaa                                                                                                                              >  2:2309442/1‑150 (MQ=255)
                         ggggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAAc                                                                                                                             <  2:1428808/150‑1 (MQ=255)
                          gggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACa                                                                                                                            >  1:1087802/1‑150 (MQ=255)
                           ggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACaa                                                                                                                           >  1:1956024/1‑150 (MQ=255)
                             cATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCGGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATg                                                                                                                         <  2:1974284/150‑1 (MQ=255)
                              aTTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATgg                                                                                                                        >  1:505156/1‑150 (MQ=255)
                               tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                       <  1:1197074/150‑1 (MQ=255)
                                 tCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcagg                                                                                                                                       <  1:926501/132‑5 (MQ=255)
                                        tACAATTTATATTCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTTGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAg                                                                                                              <  2:762871/150‑1 (MQ=255)
                                        tACAATTTATATGCCTCACATTATGCACCACAACAGGGCGATAATACGCAAAATATAAGTTCCTATGGGCACGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGAACGAACAGTCATTTAACAATGGCAAAAGCCAg                                                                                                              >  2:552583/1‑150 (MQ=255)
                                          cAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGc                                                                                                            >  2:754856/1‑150 (MQ=255)
                                          cAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGc                                                                                                            <  2:2327040/150‑1 (MQ=255)
                                               tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCCGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCa                                                                                                       <  2:2447987/150‑1 (MQ=255)
                                                  atGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAAcc                                                                                                    <  2:2166852/150‑1 (MQ=255)
                                                   tGCCTCTCATTATGCCCCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCt                                                                                                   >  1:1500543/1‑150 (MQ=255)
                                                   tGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCt                                                                                                   >  2:134275/1‑150 (MQ=255)
                                                    gCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCAGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTg                                                                                                  <  2:182802/150‑1 (MQ=255)
                                                     ccTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTgg                                                                                                 >  2:539584/1‑150 (MQ=255)
                                                      ctctCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGa                                                                                                >  1:1590062/1‑150 (MQ=255)
                                                           aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                                                                           <  2:1503508/150‑1 (MQ=255)
                                                           aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGGTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                                                                           <  1:1340476/150‑1 (MQ=255)
                                                              aTGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCcgc                                                                                        <  1:1929684/150‑1 (MQ=255)
                                                                   ccacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtatt                                                                                   >  1:812827/1‑150 (MQ=255)
                                                                    cacaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                                                                  >  1:1540140/1‑150 (MQ=255)
                                                                        tCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTatttattt                                                                              >  2:1532216/1‑150 (MQ=255)
                                                                         cAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTg                                                                             <  1:111121/150‑1 (MQ=255)
                                                                          aGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGt                                                                            >  1:888772/1‑150 (MQ=255)
                                                                           gggCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGtt                                                                           <  1:919938/150‑1 (MQ=255)
                                                                            ggCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGttt                                                                          >  2:2295004/1‑150 (MQ=255)
                                                                              cGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCg                                                                        <  1:70008/150‑1 (MQ=255)
                                                                                aTAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGcc                                                                      <  1:1331579/150‑1 (MQ=255)
                                                                                    ttCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAAt                                                                  >  1:1171772/1‑150 (MQ=255)
                                                                                      cGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATcc                                                                >  2:404028/1‑150 (MQ=255)
                                                                                      cGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATcc                                                                >  1:589417/1‑150 (MQ=255)
                                                                                       gaaaaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCATGGCGCCTGCGTACTGATTACCACTACGATCAGTCATTTAACAATGGAAAAAGCCACGCGACCAACATGGATTTTCCGCGTATTTAATTCTTTCGCCAAATccc                                                               >  2:1849081/3‑150 (MQ=255)
                                                                                        cAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCa                                                              >  2:1571394/1‑150 (MQ=255)
                                                                                        cAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATATAACAATGGCAAAAGCCAGGCGACCAACCTGGATATTCCGCGTATTTATTTGTTTCGCCCAATCCCa                                                              >  2:109738/1‑150 (MQ=255)
                                                                                         aaaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAg                                                             <  1:620930/150‑1 (MQ=255)
                                                                                            atatAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaa                                                          <  1:2050238/150‑1 (MQ=255)
                                                                                             tataAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaat                                                         <  2:480604/150‑1 (MQ=255)
                                                                                             tataAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGACGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaat                                                         <  1:2252411/150‑1 (MQ=255)
                                                                                                 aGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGTTTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaat                                                     >  2:350860/1‑150 (MQ=255)
                                                                                                  gTTGCTATGGGCAGGGTGGGGTTACTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAAGGGCAAAAGCCAGTCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                    <  1:2293770/150‑1 (MQ=255)
                                                                                                  gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                    >  2:309436/1‑150 (MQ=255)
                                                                                                  gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                    >  1:2051493/1‑150 (MQ=255)
                                                                                                  gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                    >  1:764678/1‑150 (MQ=255)
                                                                                                   ttCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGc                                                   <  1:1039034/150‑1 (MQ=255)
                                                                                                    tCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCa                                                  <  2:160336/150‑1 (MQ=255)
                                                                                                      cTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCaaa                                                <  1:2295004/150‑1 (MQ=255)
                                                                                                        aTGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaac                                              >  2:925423/1‑150 (MQ=255)
                                                                                                          gggCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                            >  1:2369743/1‑150 (MQ=255)
                                                                                                          gggCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                            <  1:748011/150‑1 (MQ=255)
                                                                                                             cAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaactaact                                         <  2:1956024/150‑1 (MQ=255)
                                                                                                               ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                       >  1:1730271/1‑150 (MQ=255)
                                                                                                               ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                       >  2:2430234/1‑150 (MQ=255)
                                                                                                                gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                      <  1:887144/150‑1 (MQ=255)
                                                                                                                gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                      <  1:968653/150‑1 (MQ=255)
                                                                                                                    gggTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTc                                  <  2:243704/150‑1 (MQ=255)
                                                                                                                         aaTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                             <  1:543099/150‑1 (MQ=255)
                                                                                                                          aTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTcc                                                                                          >  1:1463847/1‑88 (MQ=255)
                                                                                                                          aTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTcc                                                                                          <  2:1463847/88‑1 (MQ=255)
                                                                                                                            cTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                          >  2:1146756/1‑150 (MQ=255)
                                                                                                                            cTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                          <  2:1730271/150‑1 (MQ=255)
                                                                                                                              tGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATAc                        >  2:2408983/1‑150 (MQ=255)
                                                                                                                                gggCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTTGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTg                      <  2:616921/150‑1 (MQ=255)
                                                                                                                                  gCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGaa                    >  1:1249505/1‑150 (MQ=255)
                                                                                                                                  gCCAGGCGCCTGCGGACTGATTACCAGTACGATCAGTCATTAAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCACCGAATGCAAAACTAACTAAAGGTCAATACGATACTGaa                    <  2:2155268/150‑1 (MQ=255)
                                                                                                                                      ggCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCct                >  1:55308/1‑150 (MQ=255)
                                                                                                                                      ggCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCct                <  2:1087802/150‑1 (MQ=255)
                                                                                                                                           gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATcc                 <  2:516489/141‑1 (MQ=255)
                                                                                                                                           cTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtt           >  1:696658/1‑150 (MQ=255)
                                                                                                                                            tGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCTATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAttt          >  1:1217601/1‑150 (MQ=255)
                                                                                                                                             gCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtttt         >  2:689839/1‑150 (MQ=255)
                                                                                                                                             gCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtttt         <  1:533067/150‑1 (MQ=255)
                                                                                                                                             gCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtttt         >  1:1851817/1‑150 (MQ=255)
                                                                                                                                              cGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTc        <  1:1728294/150‑1 (MQ=255)
                                                                                                                                                 aCTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAc     <  2:1659582/150‑1 (MQ=255)
                                                                                                                                                  cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAct    >  1:1688041/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAACTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCT  >  NC_000913/3187840‑3188137

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