Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,693,418 C→T 5.8% P282P (CCG→CCA uidB ← glucuronide:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,693,4180CT 5.8% 499.4 / 11.9 155P282P (CCG→CCAuidBglucuronide:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (59/87);  new base T (5/4);  total (64/91)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 157 total aligned reads displayed.

CCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCACA  >  NC_000913/1693269‑1693566
                                                                                                                                                     |                                                                                                                                                    
ccAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                      >  2:61254/1‑150 (MQ=255)
 cAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCg                                                                                                                                                     <  1:1454514/150‑1 (MQ=255)
  aaTTGAAGCGACGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCgg                                                                                                                                                    <  1:37925/150‑1 (MQ=255)
   aTTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGt                                                                                                                                                   >  1:1508164/1‑150 (MQ=255)
    ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGCCCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTg                                                                                                                                                  >  2:19274/1‑150 (MQ=255)
    ttGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCACCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTg                                                                                                                                                  >  2:9781/1‑150 (MQ=255)
     tGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTCTTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGc                                                                                                                                                 >  1:65484/1‑150 (MQ=255)
     tGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGc                                                                                                                                                 >  2:638678/1‑150 (MQ=255)
         gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggg                                                                                                                                             <  1:2077951/150‑3 (MQ=255)
          cGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATg                                                                                                                                            <  2:368840/150‑1 (MQ=255)
           gATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCgggga                                                                                                                                           <  1:359865/150‑5 (MQ=255)
            aTGGTCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGcc                                                                                                                                          <  2:1769806/150‑1 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggggatgg                                                                                                                                        <  2:1155161/150‑8 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                        <  1:2149455/150‑1 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                        <  2:1420998/150‑1 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                        <  2:2091277/150‑1 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                        >  2:245316/1‑150 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                        >  1:1706401/1‑150 (MQ=255)
              ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcaa                                                                                                                                        <  1:2444065/150‑2 (MQ=255)
                ccAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAg                                                                                                                                      >  1:571096/1‑150 (MQ=255)
                  aaCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTa                                                                                                                                    <  1:182782/150‑1 (MQ=255)
                  aaCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTa                                                                                                                                    <  1:2332348/150‑1 (MQ=255)
                   aCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTac                                                                                                                                   <  2:2009017/150‑1 (MQ=255)
                     gCAACCAACGCCACCGGCAGTGACCCGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTacca                                                                                                                                 <  1:2250311/150‑1 (MQ=255)
                     gCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCTGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTacca                                                                                                                                 >  2:1246400/1‑150 (MQ=255)
                       aaccaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaac                                                                                                                               >  1:2386573/1‑150 (MQ=255)
                        accaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGGAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCGGCACCAGCGGTGCCGCTGCCACAGTcccaacc                                                                                                                              >  2:2472652/1‑150 (MQ=255)
                        accaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaacc                                                                                                                              <  2:1001461/150‑1 (MQ=255)
                        accaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaacc                                                                                                                              >  2:55090/1‑150 (MQ=255)
                         ccaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaacca                                                                                                                             <  2:786105/150‑1 (MQ=255)
                          caccGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATACCCCCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCGCGCGCCCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                            >  1:727755/4‑150 (MQ=255)
                          caacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                            >  1:1823299/1‑150 (MQ=255)
                          caacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                            >  1:1618211/1‑150 (MQ=255)
                               ccACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGtttt                                                                                                                       <  1:989147/150‑1 (MQ=255)
                                cACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTg                                                                                                                      <  2:2316644/150‑1 (MQ=255)
                                cACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTg                                                                                                                      <  2:2234203/150‑1 (MQ=255)
                                 aCCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGc                                                                                                                     <  2:704883/150‑1 (MQ=255)
                                 aCCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGGTTTGc                                                                                                                     <  1:1729522/150‑1 (MQ=255)
                                    ggCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCAcc                                                                                                                  <  1:985974/150‑1 (MQ=255)
                                      cAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAg                                                                                                                <  2:390676/150‑1 (MQ=255)
                                        gTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTa                                                                                                              <  1:2212100/150‑1 (MQ=255)
                                        gTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGCTTTGCACCAGTa                                                                                                              >  2:2247437/1‑150 (MQ=255)
                                            ccAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAg                                                                                                          >  1:2219877/1‑150 (MQ=255)
                                               gACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGcac                                                                                                       >  1:2485314/1‑150 (MQ=255)
                                                 cGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCCCAg                                                                                                     >  1:2400721/1‑150 (MQ=255)
                                                   gAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                   <  1:2384515/150‑1 (MQ=255)
                                                   gAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                   >  2:727834/1‑150 (MQ=255)
                                                   gAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTGTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                   >  1:394763/1‑150 (MQ=255)
                                                    aaaCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTGTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGa                                                                                                  <  2:394763/150‑1 (MQ=255)
                                                      aCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                                                                <  1:1066091/150‑1 (MQ=255)
                                                         cAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGc                                                                                             >  2:105597/1‑150 (MQ=255)
                                                             agaaCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCgg                                                                                         <  1:1512087/150‑1 (MQ=255)
                                                               aaCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTa                                                                                       <  2:964089/150‑1 (MQ=255)
                                                               aaCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTa                                                                                       <  2:603988/150‑1 (MQ=255)
                                                                  agcagTTAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTCCACCAGTACCAGCACAGTGAACAGCCCGGTATCa                                                                                    <  2:1609995/150‑1 (MQ=255)
                                                                    cagATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAtt                                                                                  <  1:245316/150‑1 (MQ=255)
                                                                     agATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAttt                                                                                 >  1:659151/1‑150 (MQ=255)
                                                                          aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAacac                                                                            >  1:2456372/1‑150 (MQ=255)
                                                                           ccGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAAcaca                                                                           >  2:1878771/1‑150 (MQ=255)
                                                                               aGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAgc                                                                       <  2:2016338/150‑1 (MQ=255)
                                                                                ggTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAgcg                                                                      >  2:2039078/1‑150 (MQ=255)
                                                                                     ccAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGt                                                                 <  2:95543/150‑1 (MQ=255)
                                                                                      cAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTa                                                                >  1:1829333/1‑150 (MQ=255)
                                                                                      cAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTa                                                                <  1:1265287/150‑1 (MQ=255)
                                                                                         cAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaa                                                             <  1:967385/150‑1 (MQ=255)
                                                                                           aaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaaca                                                           >  1:91545/1‑150 (MQ=255)
                                                                                              cgcCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacg                                                        <  1:495948/150‑1 (MQ=255)
                                                                                                ccAATCATCAAGGTATTCTTTATACAGATCATCGCTACCAACCCCGGCACCAGCGATCCCGATACCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacgac                                                      >  2:449099/1‑150 (MQ=255)
                                                                                                ccAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacgac                                                      <  1:1032047/150‑1 (MQ=255)
                                                                                                   aTCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAg                                                   <  2:234315/150‑1 (MQ=255)
                                                                                                     cAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGgc                                                 >  2:1418305/1‑150 (MQ=255)
                                                                                                           ggTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGAc                                           >  2:503778/1‑150 (MQ=255)
                                                                                                                 cTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCaaa                                     >  1:795006/1‑150 (MQ=255)
                                                                                                                  tttttACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAg                                    <  1:448423/150‑1 (MQ=255)
                                                                                                                   ttttCCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAgg                                   <  2:647238/150‑1 (MQ=255)
                                                                                                                   ttttACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAgg                                   >  2:1961172/1‑150 (MQ=255)
                                                                                                                      tACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCg                                >  2:923386/1‑150 (MQ=255)
                                                                                                                       aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGa                               <  2:419443/150‑1 (MQ=255)
                                                                                                                          gATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCACCCAGGTTTTGCACCAGTACCAGCACAGTGCACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGCAAt                            >  1:793991/1‑150 (MQ=255)
                                                                                                                          gATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAAt                            >  1:7162/1‑150 (MQ=255)
                                                                                                                          gATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAAt                            >  1:1297777/1‑150 (MQ=255)
                                                                                                                              cTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTTTCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGc                        <  2:459664/150‑1 (MQ=255)
                                                                                                                              cTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGc                        >  1:829600/1‑150 (MQ=255)
                                                                                                                              cTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGc                        <  2:2456372/150‑1 (MQ=255)
                                                                                                                              cTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGc                        <  2:2376923/150‑1 (MQ=255)
                                                                                                                                  cgACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGcacac                    <  1:533800/150‑1 (MQ=255)
                                                                                                                                    aCCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGcacacac                  <  2:1810793/150‑1 (MQ=255)
                                                                                                                                     ccATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGcacacaca                 <  2:2485314/150‑1 (MQ=255)
                                                                                                                                     ccATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGcacacaca                 <  2:727755/150‑1 (MQ=255)
                                                                                                                                        gcatcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgc               >  1:534569/5‑149 (MQ=255)
                                                                                                                                        gcatcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgc               <  2:534569/145‑1 (MQ=255)
                                                                                                                                         catcGGCACCGCTGGTGCCGATGCCCCAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACCCACAgcgc             >  2:1253597/4‑150 (MQ=255)
                                                                                                                                         catcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGCACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgcgc             >  1:1362547/4‑150 (MQ=255)
                                                                                                                                         catcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgcgc             >  2:432481/4‑150 (MQ=255)
                                                                                                                                          cccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCa            >  2:1521665/1‑150 (MQ=255)
                                                                                                                                            cGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAAGAGATAGCTCACTTAGATCTACGACGAGGCGCTGACCTCATAGGTAGAAATCAGCAAACACAgggttcc          <  2:1652924/150‑7 (MQ=255)
                                                                                                                                               cACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGAt       >  1:569372/1‑150 (MQ=255)
                                                                                                                                                 ccAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCAAACACAGCGCTCCGATGc     <  2:1761150/150‑1 (MQ=255)
                                                                                                                                                    gCGGTGCCGATGCCACAGTACCAACCAGTTTTTGCACCAGTACGAGCAGAGTGAACAGCCGGGTATGATTTAACTCATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCTAAGGTCGAAATCAGCACACACTGCGCACCGATGcaca  <  2:1823299/150‑1 (MQ=255)
                                                                                                                                                    gCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGcaca  <  1:727834/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCACA  >  NC_000913/1693269‑1693566

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: