Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,187,993 C→T 5.7% pseudogene (562/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,9930CT 5.7% 498.9 / 12.6 156pseudogene (562/2434 nt)yqiGputative outer membrane usher protein YqiG
Reads supporting (aligned to +/- strand):  ref base C (73/74);  new base T (4/5);  total (77/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 164 total aligned reads displayed.

CTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCA  >  NC_000913/3187844‑3188142
                                                                                                                                                     |                                                                                                                                                     
cTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                                                                       <  1:79720/150‑1 (MQ=255)
cTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                                                                       >  2:111121/1‑150 (MQ=255)
cTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAc                                                                                                                                                       >  1:137211/1‑150 (MQ=255)
 tttGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAcc                                                                                                                                                      <  1:599944/150‑1 (MQ=255)
 tttGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTAcc                                                                                                                                                      <  1:1785245/150‑1 (MQ=255)
     gAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTa                                                                                                                                                  >  1:1026042/1‑150 (MQ=255)
     gAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTa                                                                                                                                                  >  2:1231472/1‑150 (MQ=255)
     gAATCATGGTATTGATGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATATTGGGGCCTGGCGCCTGCGTACTGATTACCAGTa                                                                                                                                                  >  2:2307143/1‑150 (MQ=255)
       aTCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACg                                                                                                                                                >  2:1331579/1‑150 (MQ=255)
        tCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgc                                                                                                                                               <  2:1081554/150‑2 (MQ=255)
         cATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgca                                                                                                                                              <  1:316419/150‑3 (MQ=255)
          aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                             <  1:2315610/150‑1 (MQ=255)
          aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                             >  2:70008/1‑150 (MQ=255)
          aTGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATc                                                                                                                                             <  1:2451907/150‑1 (MQ=255)
           tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCa                                                                                                                                            <  2:591321/150‑1 (MQ=255)
           tGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGGACTGATTACCAGTACGATCa                                                                                                                                            <  2:950237/150‑1 (MQ=255)
            tgTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAg                                                                                                                                           <  1:568890/149‑1 (MQ=255)
            ggTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGGTAATCAGTACgcaggc                                                                                                                                           <  1:65501/150‑6 (MQ=255)
            ggTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAg                                                                                                                                           >  1:480604/1‑150 (MQ=255)
            ggTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAg                                                                                                                                           <  1:1873394/150‑1 (MQ=255)
             gTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGt                                                                                                                                          <  1:573178/150‑1 (MQ=255)
                 tGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGAAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCATGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAtt                                                                                                                                      >  2:719315/1‑150 (MQ=255)
                  gCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCAttt                                                                                                                                     <  2:1923265/150‑1 (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:243704/1‑150 (MQ=255)
                       ggCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACaa                                                                                                                                >  1:1956024/1‑150 (MQ=255)
                         cATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCGGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATg                                                                                                                              <  2:1974284/150‑1 (MQ=255)
                           tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                            <  2:470225/150‑1 (MQ=255)
                           tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                            <  1:1197074/150‑1 (MQ=255)
                           tttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGc                                                                                                                            <  2:1193960/150‑1 (MQ=255)
                            ttCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCa                                                                                                                           <  2:1786125/150‑1 (MQ=255)
                                    tACAATTTATATTCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTTGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAg                                                                                                                   <  2:762871/150‑1 (MQ=255)
                                    tACAATTTATATGCCTCACATTATGCACCACAACAGGGCGATAATACGCAAAATATAAGTTCCTATGGGCACGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGAACGAACAGTCATTTAACAATGGCAAAAGCCAg                                                                                                                   >  2:552583/1‑150 (MQ=255)
                                     aCAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCATGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAgg                                                                                                                  <  2:655884/150‑1 (MQ=255)
                                      cAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGc                                                                                                                 >  2:500808/1‑150 (MQ=255)
                                       aaTTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCg                                                                                                                <  1:922049/150‑1 (MQ=255)
                                           tataTGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCCGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCa                                                                                                            <  2:2447987/150‑1 (MQ=255)
                                               tGCCTCTCATTATGCCCCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCt                                                                                                        >  1:1500543/1‑150 (MQ=255)
                                                gCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCAGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTg                                                                                                       <  2:182802/150‑1 (MQ=255)
                                                  ctctCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGa                                                                                                     >  1:1590062/1‑150 (MQ=255)
                                                  ctctCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGa                                                                                                     >  2:365537/1‑150 (MQ=255)
                                                   tctcATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGAt                                                                                                    <  1:2252607/150‑1 (MQ=255)
                                                       aTTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGGTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTc                                                                                                <  1:1340476/150‑1 (MQ=255)
                                                        ttATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTGATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTcc                                                                                               <  1:1554765/150‑1 (MQ=255)
                                                             caccacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGTCAGGCGACCAACCTGGATTTTCCGCGta                                                                                          >  1:1679506/1‑150 (MQ=255)
                                                              accacATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtat                                                                                         <  2:736040/150‑1 (MQ=255)
                                                                cacaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                                                                       >  1:1540140/1‑150 (MQ=255)
                                                                cacaTCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGtattt                                                                                       >  2:1468784/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)
                                                                            aTAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGcc                                                                           <  1:1331579/150‑1 (MQ=255)
                                                                             tAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGccc                                                                          <  2:680270/150‑1 (MQ=255)
                                                                                ttCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGTTCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAAt                                                                       >  1:906756/1‑150 (MQ=255)
                                                                                ttCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAAt                                                                       >  1:1171772/1‑150 (MQ=255)
                                                                                  cGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATcc                                                                     >  2:404028/1‑150 (MQ=255)
                                                                                   gaaaaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCATGGCGCCTGCGTACTGATTACCACTACGATCAGTCATTTAACAATGGAAAAAGCCACGCGACCAACATGGATTTTCCGCGTATTTAATTCTTTCGCCAAATccc                                                                    >  2:1849081/3‑150 (MQ=255)
                                                                                   gCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAATCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATccc                                                                    >  1:549101/1‑150 (MQ=255)
                                                                                     aaaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAg                                                                  <  1:500808/150‑1 (MQ=255)
                                                                                      aaaTATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGc                                                                 <  2:2389283/150‑1 (MQ=255)
                                                                                         tataAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGCCCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaat                                                              >  2:218807/1‑150 (MQ=255)
                                                                                         tataAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGACGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaat                                                              <  1:2252411/150‑1 (MQ=255)
                                                                                              gTTCCTATGGGGAGGCTGGCCTTAAACTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCTGTCATTTAACAATGGCAACAGCCAGGCGACCACCCTGGATTTTCCGCGTAGTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                         >  2:850660/1‑150 (MQ=255)
                                                                                              gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                         >  1:2051493/1‑150 (MQ=255)
                                                                                              gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                         >  1:764678/1‑150 (MQ=255)
                                                                                              gTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                         >  2:309436/1‑150 (MQ=255)
                                                                                              cttCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCaatgaatg                                                         >  1:985166/2‑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                                                   <  1:2481826/150‑1 (MQ=255)
                                                                                                    aTGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaac                                                   >  2:925423/1‑150 (MQ=255)
                                                                                                      gggCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAaacta                                                 >  1:869698/1‑150 (MQ=255)
                                                                                                           ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                            >  1:1730271/1‑150 (MQ=255)
                                                                                                           ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                            >  2:543099/1‑150 (MQ=255)
                                                                                                           ggCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtat                                            >  2:2430234/1‑150 (MQ=255)
                                                                                                            gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                           >  2:1512935/1‑150 (MQ=255)
                                                                                                            gCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACtata                                           <  1:968653/150‑1 (MQ=255)
                                                                                                                     aaTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATAc                                  <  1:543099/150‑1 (MQ=255)
                                                                                                                      aTCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTcc                                                                                               >  1:1463847/1‑88 (MQ=255)
                                                                                                                        cTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGAt                               <  2:1730271/150‑1 (MQ=255)
                                                                                                                          tGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATAc                             >  2:2408983/1‑150 (MQ=255)
                                                                                                                           ggggCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACt                            >  2:1364120/1‑150 (MQ=255)
                                                                                                                            gggCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTg                           <  2:2441196/150‑1 (MQ=255)
                                                                                                                            gggCCGGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTg                           <  2:2154368/150‑1 (MQ=255)
                                                                                                                              gCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGaa                         >  1:1249505/1‑150 (MQ=255)
                                                                                                                              gCCAGGCGCCTGCGGACTGATTACCAGTACGATCAGTCATTAAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCACCGAATGCAAAACTAACTAAAGGTCAATACGATACTGaa                         <  2:2155268/150‑1 (MQ=255)
                                                                                                                                  ggCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCct                     <  2:1087802/150‑1 (MQ=255)
                                                                                                                                      tgatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAt                 >  2:1579572/5‑150 (MQ=255)
                                                                                                                                       gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATcc                      <  2:516489/141‑1 (MQ=255)
                                                                                                                                       gatcGTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATcc                      >  1:516489/4‑144 (MQ=255)
                                                                                                                                       cTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAtt                >  1:696658/1‑150 (MQ=255)
                                                                                                                                        tGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCTATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAttt               >  1:1217601/1‑150 (MQ=255)
                                                                                                                                           gTACTGGTAATCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTAttttag            >  1:729177/1‑148 (MQ=255)
                                                                                                                                            tACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGCCGACCAACCTGGATTTTCCGCCTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGa           >  1:2019845/1‑150 (MQ=255)
                                                                                                                                             aCTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAc          <  2:1659582/150‑1 (MQ=255)
                                                                                                                                              cTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGAct         >  1:1688041/1‑150 (MQ=255)
                                                                                                                                                  ttACCAGTACGATCAGTCATTTAACAATGTCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAACCTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCttt     <  1:1726813/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTGGAATCATGGTATTGCTGGGGCATTTCTTGACTACAATTTATATGCCTCTCATTATGCACCACATCAGGGCGATAATTCGCAAAATATAAGTTCCTATGGGCAGGCTGGGGTTAATCTTGGGGCCTGGCGCCTGCGTACTGATTACCAGTACGATCAGTCATTTAACAATGGCAAAAGCCAGGCGACCAACCTGGATTTTCCGCGTATTTATTTGTTTCGCCCAATCCCAGCAATGAATGCAAAACTAACTATAGGTCAATACGATACTGAATCCTCTATTTTCGACTCTTTCCA  >  NC_000913/3187844‑3188142

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