Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,392 T→A 8.7% intergenic (‑53/‑31) ybjJ ← / → rcdA putative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913887,3920TA 8.7% 204.5 / 7.2 79intergenic (‑53/‑31)ybjJ/rcdAputative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator
Reads supporting (aligned to +/- strand):  ref base T (32/40);  new base A (5/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGG  >  NC_000913/887246‑887540
                                                                                                                                                  |                                                                                                                                                    
cccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGGCCAGGTTCGGCGTTTCAATGCATTACGTGAATAATTTACGGTCATAGATCACGTCAAATTAAGAAGAGGAAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                   <  2:30041/150‑1 (MQ=255)
cccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                   <  2:35481/150‑1 (MQ=255)
cccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                   >  1:1241201/1‑150 (MQ=255)
cccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                   <  2:622476/150‑1 (MQ=255)
cccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                   <  2:39466/150‑1 (MQ=255)
  catgACGCCATTAACAGGCCTTGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCAAGTCAAAATAAGAAGAGGGTAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                 <  2:238789/147‑1 (MQ=255)
    ggACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTaca                                                                                                                                               >  1:130753/1‑150 (MQ=255)
     gACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                                                                     >  1:58757/1‑143 (MQ=255)
     gACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                                                                     <  2:58757/143‑1 (MQ=255)
     gACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                                                                              <  2:1259734/150‑1 (MQ=255)
         ccATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAAtt                                                                                                                                          >  1:865069/1‑150 (MQ=255)
         ccATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                                                                          <  2:1077029/150‑3 (MQ=255)
         ccAATAACACCCCTCGCAAAAAGACGAACAAAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACCTGAACAAATTACGGTCATAGATCACCTCAAAATAAGAAGAGCGAAGACAACATAAGCAAGGTTGTGTACAAATGTACACaacc                                                                                                                                          <  1:1027183/150‑3 (MQ=255)
          cATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATATGTACACAATTg                                                                                                                                         >  1:571950/1‑150 (MQ=255)
          cATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTGCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                                                                         >  1:665956/1‑150 (MQ=255)
              aaCAGGCCTGGCAAAAAGCAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATTAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAGACAATTGCaga                                                                                                                                     >  2:618893/1‑150 (MQ=255)
                 aGGCCGGGCAAAAAGAAGCACCTAAACAGCGCCCCGGCTCGGCGTTTCAATGCTTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGACGACAACCTTCGCAAGGTGGTGTACATTTGTACACAATTGCagagag                                                                                                                                  <  1:793901/150‑1 (MQ=255)
                  ggCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                 >  2:488314/1‑150 (MQ=255)
                   gCCTGGCAAAAAAAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGTGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                                                                                <  1:1182845/150‑1 (MQ=255)
                      tGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAAATTACGGTCAAAGATCAAGTCAAAATAAGAAGAGGTAAGACAACATTAGCAAGGTTGTGAACATTTGTACACAATTGCAGAGACATGaa                                                                                                                             >  2:696948/1‑150 (MQ=255)
                       ggCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                                                                            >  2:790497/1‑150 (MQ=255)
                           aaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACACCATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGTAAATGAc                                                                                                                        >  1:1148256/1‑150 (MQ=255)
                           aaaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                                                                        <  2:1208527/150‑1 (MQ=255)
                             aagaagAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                                                                      >  2:1075431/1‑150 (MQ=255)
                                  gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTAATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTc                                                                                                                 <  2:1119657/150‑1 (MQ=255)
                                          aCAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACAGTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAAATGAGTCGCGCTaaa                                                                                                         >  1:113741/1‑149 (MQ=255)
                                             gcgcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGTAGAGGGAAGACAACATTAGCATGGTTGTGTACATTTTTAAACAATTGCATAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                                                                      >  2:231758/1‑150 (MQ=255)
                                             gcgcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                                                                      >  2:706356/1‑150 (MQ=255)
                                               gcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATcc                                                                                                    <  2:318848/150‑1 (MQ=255)
                                                cccAGGTGCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                                                   <  2:29393/150‑1 (MQ=255)
                                                 ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGc                                                                                                  >  2:1030528/1‑150 (MQ=255)
                                                  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                                                                                 <  2:1198853/150‑1 (MQ=255)
                                                    ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCATAAGAAGAAGAGGGAAGACAACATTATCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCCCTCTAACGATCCTCAAc                                                                                               >  1:85726/1‑150 (MQ=255)
                                                       tCGGCGTTTCAATGCATTACGTTAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                                                                            <  2:206088/150‑1 (MQ=255)
                                                        cgtcGTTTCAATGCATTACGTGAATAATTTACGGTCATAGATCACGTCAAAATACGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGTGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                                           <  1:927302/147‑1 (MQ=255)
                                                        cGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                                           <  1:414880/150‑1 (MQ=255)
                                                          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                                                                         <  1:1076779/150‑1 (MQ=255)
                                                           cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGGCGCGCTAACGATCCGCAACGGCgaga                                                                                        <  1:781985/150‑1 (MQ=255)
                                                               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaa                                                                                    >  1:264066/1‑150 (MQ=255)
                                                               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaa                                                                                    >  1:521677/1‑150 (MQ=255)
                                                                 aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTAAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAt                                                                                  <  2:957664/150‑1 (MQ=255)
                                                                  aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAtt                                                                                 <  1:281645/150‑1 (MQ=255)
                                                                              aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACAATAGAAAGGTTGTGTTCATTTGTAAAAAATTGCATAGAGAGTAAATGACATGCGTCGCGCTAACGATCCGCAACGGAGAGAAAAAATTATCCAGTACACg                                                                     >  1:75914/1‑150 (MQ=255)
                                                                                gaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                                                                   >  1:1201052/1‑150 (MQ=255)
                                                                                  aTTTACGGGCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAATGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg                                                                 <  1:56808/150‑1 (MQ=255)
                                                                                   tttACGTTCATAGCTCACGTCAAAATATGAAGAGGGAAGACTACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGTCATGCGTCGCGCTAACGATCCGCAACGGCGAGAATAAATTATCGAGGCCACGCTGGa                                                                <  2:429150/150‑1 (MQ=255)
                                                                                    ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg                                                               >  1:723701/1‑150 (MQ=255)
                                                                                    ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg                                                               <  1:482078/150‑1 (MQ=255)
                                                                                          gcATAGATCAAGTCAAAATAAGAAGAGGGAAGACAACTTTAGCAAGCTTGTGTACATTTGTACTCAATTGCAGAGAGACGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg                                                         <  2:506412/149‑1 (MQ=255)
                                                                                               tatgaCGTAAAAATAATAAGAGGTAAGACATCATTAGCAAGGTTGTGTACATTTGTACACTATTGCAGAGAGAGGAAATGACATGCGTCGCGCTATCGTTCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACt                                                    <  1:880479/146‑1 (MQ=255)
                                                                                                  cACGTCAAAATAAGAAGAGGGAAGGCAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                                 >  1:960163/1‑150 (MQ=255)
                                                                                                  cACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                                 >  1:1238863/1‑150 (MQ=255)
                                                                                                  cACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                                 >  1:729855/1‑150 (MQ=255)
                                                                                                     gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACgg                                              <  1:670068/150‑1 (MQ=255)
                                                                                                      tCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGa                                             <  2:665956/150‑1 (MQ=255)
                                                                                                          aaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGATTTGAAACTTTACGGAATAc                                         >  1:849752/1‑150 (MQ=255)
                                                                                                             aagaagAGGGAAGACAACATTAGCAAGGTGGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATg                                      <  1:3278/150‑1 (MQ=255)
                                                                                                                        aGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACTCAc                           <  1:1188392/150‑1 (MQ=255)
                                                                                                                          acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCg                         <  2:542480/150‑1 (MQ=255)
                                                                                                                           caacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGc                        >  1:1298000/1‑150 (MQ=255)
                                                                                                                               aTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCaaaa                    <  2:253428/150‑1 (MQ=255)
                                                                                                                                 tAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAtt                  <  1:529458/150‑1 (MQ=255)
                                                                                                                                 tAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAACCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAAtt                  >  1:988393/1‑150 (MQ=255)
                                                                                                                                  aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTg                 <  2:284874/150‑1 (MQ=255)
                                                                                                                                  aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTg                 >  1:681121/1‑150 (MQ=255)
                                                                                                                                   ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCATAATTGc                >  1:452551/7‑150 (MQ=255)
                                                                                                                                     gcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTa              >  1:610004/5‑150 (MQ=255)
                                                                                                                                     gcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTa              >  1:615906/5‑150 (MQ=255)
                                                                                                                                       ggTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGCCATGCGTCGCGCGAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGTATACATGCGGTTACGCACCGCAAAATTGCTAcc            <  1:328422/150‑1 (MQ=255)
                                                                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTAccc           <  2:650673/150‑1 (MQ=255)
                                                                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGACACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTAccc           <  2:1160575/150‑1 (MQ=255)
                                                                                                                                        attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGGAACGGCGAGAAAAAATTATCCAGGCAACGCTGGATCCGGTGAAACTTTACGGAATACATGCTTTTACGCTCCGCAAAATTGCTTccc           >  2:61213/2‑150 (MQ=255)
                                                                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTCTTACGCACCGCAAAATTGCTACCCTTGc       >  2:713339/1‑150 (MQ=255)
                                                                                                                                            tgtACAAATGTACACTATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGAAAAGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGAGAAACTTTACGGAATACTTGCTGTTACTCACCGCAAAAATGCTACCCttcc       >  1:175160/1‑148 (MQ=255)
                                                                                                                                             gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGGCGCGCTAACGATCCGCAACGGCGAGAAAGAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGcc      <  2:248048/150‑1 (MQ=255)
                                                                                                                                                cATTTGTACACAATTGCAGAGAGAGGTAATGACATGCGTCGCGCTAACGATCCGCAACGCCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTAACGGAATACATGCTGTTACGCAACGCAAAATTGCTACCCTTGCCggg   >  1:246090/1‑150 (MQ=255)
                                                                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCggg   <  1:422580/150‑1 (MQ=255)
                                                                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCggg   >  2:731280/1‑150 (MQ=255)
                                                                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCggg   >  2:133723/1‑150 (MQ=255)
                                                                                                                                                 aTTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCgggg  <  2:1238863/150‑1 (MQ=255)
                                                                                                                                                 aTTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCgggg  <  2:1087456/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGG  >  NC_000913/887246‑887540

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