Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,441,062 G→A 5.7% S205F (TCC→TTC) ‡ rpsD ← 30S ribosomal subunit protein S4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,441,0620GA 5.7% 496.4 / 10.9 158S205F (TCC→TTC) ‡rpsD30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base G (77/70);  new base A (6/3);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 162 total aligned reads displayed.

CAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCT  >  NC_000913/3440913‑3441208
                                                                                                                                                     |                                                                                                                                                  
cAGAGGATGGCCCACGCCACGATCTAAAGGCGCAAGGGGCACCGTGGCGGGCGTCGAACTCACTTGCGCGATATCAACCAGGCGCGGGTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg                                                                                                                                                    <  2:1169921/150‑1 (MQ=255)
 agagTATGGCCAAAGCCACGCTCTAACGGCTCAAGTGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa                                                                                                                                                   <  2:806426/150‑1 (MQ=255)
  gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                                                                                                                                                  <  2:1536049/150‑1 (MQ=255)
  gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                                                                                                                                                  <  2:1651980/150‑1 (MQ=255)
  gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                                                                                                                                                  <  2:416594/150‑1 (MQ=255)
   agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGGCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGt                                                                                                                                                 <  2:982227/150‑1 (MQ=255)
       tGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAg                                                                                                                                             <  2:1105617/150‑1 (MQ=255)
         gCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt                                                                                                                                           >  1:391078/1‑150 (MQ=255)
            aaaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGTCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCttag                                                                                                                                        <  2:1412447/150‑4 (MQ=255)
             aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                       >  2:615734/1‑150 (MQ=255)
             aaGCCAAGCTCTTAAGGCTCCATGGTTCCCTTTGCGGGGGTCGAACCCACTTGCCCGCTATGAACCCGGCCGCGTTTTAGAAACCCTGGCACAGTACCCCGCATTTTGTTCCCTCTTTGGTACTAAGCTTTACTCCAAGTAAAggttagt                                                                                                                                       <  2:1472875/150‑7 (MQ=255)
              aGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgacg                                                                                                                                      <  2:281355/150‑1 (MQ=255)
              aGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCACTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgacg                                                                                                                                      <  1:259529/150‑1 (MQ=255)
                ccACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGAt                                                                                                                                    >  1:31627/1‑150 (MQ=255)
                 cACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGCACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATc                                                                                                                                   >  1:1646790/1‑150 (MQ=255)
                   cGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAg                                                                                                                                 <  1:534988/150‑1 (MQ=255)
                    gCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAgg                                                                                                                                <  1:534745/150‑1 (MQ=255)
                     ctctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt                                                                                                                               >  1:1372172/1‑150 (MQ=255)
                     ctctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt                                                                                                                               >  1:1303884/1‑150 (MQ=255)
                       ctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGCACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgtgt                                                                                                                             >  1:922141/1‑150 (MQ=255)
                         aaaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTc                                                                                                                           >  1:233182/1‑150 (MQ=255)
                           aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCCGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGt                                                                                                                         <  2:1354480/150‑1 (MQ=255)
                           aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGTTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGt                                                                                                                         <  1:1060434/150‑1 (MQ=255)
                            ggCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGtt                                                                                                                        <  1:330343/150‑1 (MQ=255)
                            ggCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGtt                                                                                                                        <  1:809215/150‑1 (MQ=255)
                             gCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTa                                                                                                                       >  1:686409/1‑150 (MQ=255)
                             gCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCGCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCTACGATCAGGTGGTCGTTa                                                                                                                       >  2:1489986/1‑150 (MQ=255)
                              cTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTaa                                                                                                                      >  1:1009122/1‑150 (MQ=255)
                              cTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTaa                                                                                                                      <  1:1259096/150‑1 (MQ=255)
                                 aaGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGATGATCAGGTGTTCGTTAATGt                                                                                                                   <  1:1041542/150‑1 (MQ=255)
                                   gggTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCAGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTcc                                                                                                                 >  1:841482/1‑150 (MQ=255)
                                    ggTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGGTTTAGAAACTCAGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCAATACTTGGAGTAAATCACGACGATCAGGTCTTCGTTAATGTACg                                                                                                                >  2:1285635/1‑150 (MQ=255)
                                      tCACCTTGGCGTGCGTCGAACTCACTTTCTCGATATCAACAAGGCGCGGTTTTAGAAACTCTGTCACAGAAGCCTGCATTGTCTCCTCTCTTTGGGACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGca                                                                                                              >  2:929226/1‑150 (MQ=255)
                                       cACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcag                                                                                                             >  1:1185776/1‑150 (MQ=255)
                                            tGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagacaga                                                                                                        <  2:1855963/150‑1 (MQ=255)
                                               cgtgcgtCGAACTCACTTTCTCGATATCAACCAGTCGCTGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTGTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCa                                                                                                     <  1:267102/150‑1 (MQ=255)
                                               cgtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCa                                                                                                     >  2:899480/1‑150 (MQ=255)
                                                gtgcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAg                                                                                                    >  1:1623568/1‑150 (MQ=255)
                                                  gcgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCATGCATTGTGTCCTCTCTTTGGTACTAAGCTTTAGTTGGAGTAAAGCTCGACGATCAGGGGTTCGTTAATGTCCGCAGCCAGATCAGaa                                                                                                  >  2:1593504/1‑150 (MQ=255)
                                                   cgtCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAAc                                                                                                 >  2:1448571/1‑150 (MQ=255)
                                                       gAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTc                                                                                             <  2:332706/150‑1 (MQ=255)
                                                             aCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCTTTAATGTCCGCAGACAGATCAGAACGCTCCGGCtt                                                                                       >  1:883811/1‑150 (MQ=255)
                                                              cTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTa                                                                                      <  2:1255616/150‑1 (MQ=255)
                                                              cTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGGCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTa                                                                                      <  1:335532/150‑1 (MQ=255)
                                                                tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACg                                                                                    >  2:112539/1‑150 (MQ=255)
                                                                 gCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGc                                                                                   <  2:901845/150‑1 (MQ=255)
                                                                  cTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGCCAGATCAGAACGCTCAGGCTTACGCt                                                                                  >  1:1703184/1‑150 (MQ=255)
                                                                  cTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCt                                                                                  >  2:790919/1‑150 (MQ=255)
                                                                      atatCAACCAGGCGCTGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTaaa                                                                              <  2:99125/150‑1 (MQ=255)
                                                                       tatCAACCAGGCGCGGTTTTAGACACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAAc                                                                             <  1:138035/150‑1 (MQ=255)
                                                                            aCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAcc                                                                        >  1:989588/1‑150 (MQ=255)
                                                                              cAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCtt                                                                      >  2:363174/1‑150 (MQ=255)
                                                                                  cgcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCAt                                                                  >  1:967583/1‑150 (MQ=255)
                                                                                  cgcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCAt                                                                  <  2:1023986/150‑1 (MQ=255)
                                                                                   gcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATc                                                                 <  2:1038576/150‑1 (MQ=255)
                                                                                       ttttAGACACTCTGTCACAGAACCCTGCGTTGTGTACTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGAGTTCGTTAATGTCCGCAGAAAGATCAGAACGCTCCGGCTTACGATTAAACGTACCTTCCATCTTGc                                                             >  1:767527/1‑150 (MQ=255)
                                                                                         ttAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCa                                                           <  1:516332/150‑1 (MQ=255)
                                                                                          tAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAg                                                          >  1:1212300/1‑150 (MQ=255)
                                                                                            gAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCa                                                        <  2:977923/150‑1 (MQ=255)
                                                                                            gAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCa                                                        >  1:639969/1‑150 (MQ=255)
                                                                                             aaaCTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAATGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCAt                                                       <  1:1845208/150‑1 (MQ=255)
                                                                                                  ctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACt                                                  >  2:129380/1‑150 (MQ=255)
                                                                                                  ctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACt                                                  >  1:1655443/1‑150 (MQ=255)
                                                                                                   tGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACtt                                                 >  1:673761/1‑150 (MQ=255)
                                                                                                     tCACAGCACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCCTCTTGCCAGCATCAACTTcc                                               >  2:1752180/1‑150 (MQ=255)
                                                                                                      ctcaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGTTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                              >  1:562114/3‑150 (MQ=255)
                                                                                                      cacaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGAACATTCCATATTTCCAGCATCAACTTcca                                              >  2:174669/1‑150 (MQ=255)
                                                                                                             cccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGt                                       >  1:1102033/1‑150 (MQ=255)
                                                                                                               cTGCATTGTGTCGTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGGTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTg                                     <  1:1534996/150‑1 (MQ=255)
                                                                                                                 tcATTGTTGCCTATCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGATCCGGCGTACGCTTAAACGTACCTTCCATCTTGGCAGCATCAACTTCCAGCCAGGTTGGc                                   <  1:296317/149‑1 (MQ=255)
                                                                                                                 gCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGc                                   <  1:211867/150‑1 (MQ=255)
                                                                                                                 gCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGc                                   <  1:433002/150‑1 (MQ=255)
                                                                                                                 gCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGc                                   <  1:1513866/150‑1 (MQ=255)
                                                                                                                    ttGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAATCTCGACGATCAGGAGTTCGTTAATGTCCGCAGACAGATCTGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCttt                                <  2:240669/150‑1 (MQ=255)
                                                                                                                    ttGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCttt                                <  1:112539/150‑1 (MQ=255)
                                                                                                                      gtgtCCTCTCTTTGGTACTAAGCTTTACTTTGAGTAAAGCTCGACGATCAGGTGTTCGTTACTGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTc                              <  2:1380181/150‑1 (MQ=255)
                                                                                                                        gtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCAc                            >  1:1772161/1‑150 (MQ=255)
                                                                                                                          ccTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                          <  2:609023/150‑1 (MQ=255)
                                                                                                                          ccTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATAAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATAAACATCCAGACAGGTTGGCTTTTCACga                          >  1:1725032/1‑149 (MQ=255)
                                                                                                                           ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgct                         >  2:1416617/1‑150 (MQ=255)
                                                                                                                           ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgct                         >  1:874991/1‑150 (MQ=255)
                                                                                                                            tctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctg                        >  1:1833880/1‑150 (MQ=255)
                                                                                                                               ctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGAAAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGATTTTCACGCTGCTc                     >  1:1815156/1‑150 (MQ=255)
                                                                                                                                 ttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGTTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcag                   <  2:938964/150‑1 (MQ=255)
                                                                                                                                 ttGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcag                   >  1:1131186/1‑150 (MQ=255)
                                                                                                                                   ggTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagcc                 >  2:1678070/1‑150 (MQ=255)
                                                                                                                                   ggTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagcc                 >  2:333958/1‑150 (MQ=255)
                                                                                                                                   ggTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGATCCGGCTTACGCTAAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGCTTGGCTTTTCACGCTGCTcagcc                 >  1:1731894/1‑150 (MQ=255)
                                                                                                                                      gtcgaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc              >  1:177026/5‑150 (MQ=255)
                                                                                                                                      aCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc              >  2:510691/1‑150 (MQ=255)
                                                                                                                                         aaGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc           >  1:855597/1‑150 (MQ=255)
                                                                                                                                         aaGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGGTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc           <  1:1446499/150‑1 (MQ=255)
                                                                                                                                          aGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCa          >  1:1000964/1‑150 (MQ=255)
                                                                                                                                           gCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCag         <  1:1448571/150‑1 (MQ=255)
                                                                                                                                             tttACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCagag       <  2:1473365/150‑1 (MQ=255)
                                                                                                                                             tttACTTGCAGTAAAGCTCGACGATCAGGCGATCGGTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTCGGCTTTTCACGCTGCTCAGCCAGCTCCagag       <  2:922141/150‑1 (MQ=255)
                                                                                                                                               tACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCg     <  2:1719061/150‑1 (MQ=255)
                                                                                                                                                 cTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGc   >  2:63239/1‑150 (MQ=255)
                                                                                                                                                 cTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGc   <  2:868769/150‑1 (MQ=255)
                                                                                                                                                  ttGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCt  <  1:64075/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCT  >  NC_000913/3440913‑3441208

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: ‑