Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,441,060 T→C 10.5% K206E (AAG→GAG)  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,0600TC10.5% 263.9 / 16.6 95K206E (AAG→GAG) rpsD30S ribosomal subunit protein S4
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base C (6/4);  total (52/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440911‑3441209
                                                                                                                                                     |                                                                                                                                                     
cccAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt                                                                                                                                                       >  1:6678/1‑150 (MQ=255)
 ccAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg                                                                                                                                                      >  2:1651572/1‑150 (MQ=255)
 ccAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg                                                                                                                                                      >  1:77875/1‑150 (MQ=255)
  cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATGAACCAGGCGCGGTTTTAGAAACTCTGTCAGAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg                                                                                                                                                     <  1:1105126/150‑1 (MQ=255)
  cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg                                                                                                                                                     <  1:915460/150‑1 (MQ=255)
  cAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg                                                                                                                                                     <  1:136386/150‑1 (MQ=255)
    gagTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                                                                                                                                                   >  1:1155737/1‑150 (MQ=255)
          ggCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc                                                                                                                                             <  2:1470931/150‑1 (MQ=255)
            ccAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTc                                                                                                                                           >  2:117036/1‑150 (MQ=255)
             cAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCTCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCtta                                                                                                                                          <  2:1088397/150‑3 (MQ=255)
              aaaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCttag                                                                                                                                         <  1:1592733/150‑4 (MQ=255)
               aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                        <  1:1070922/150‑1 (MQ=255)
               aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                        <  2:1529598/150‑1 (MQ=255)
               aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                        <  2:1269558/150‑1 (MQ=255)
               aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgac                                                                                                                                        <  1:195090/150‑1 (MQ=255)
               aaGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCAcgac                                                                                                                                        <  1:1087300/150‑1 (MQ=255)
                atccACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCttagta                                                                                                                                       <  1:719971/148‑6 (MQ=255)
                aGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTcgacg                                                                                                                                       <  1:583692/150‑1 (MQ=255)
                aGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTCCAAGTAAAGCttagta                                                                                                                                       <  1:45907/150‑6 (MQ=255)
                    aCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCa                                                                                                                                   <  1:189679/150‑1 (MQ=255)
                         ctAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgtgt                                                                                                                              >  2:1487690/1‑150 (MQ=255)
                           aaaGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGCACACTGCATTGTGTCCTCTGTTAGGTACTAAGCTTTACTTGGAGTAAAGCTCGAAGATCAGGTGTTc                                                                                                                            >  1:1416971/1‑150 (MQ=255)
                             aGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGt                                                                                                                          >  1:597726/1‑150 (MQ=255)
                              ggCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGtt                                                                                                                         >  1:599537/1‑150 (MQ=255)
                                cTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTaa                                                                                                                       >  2:216875/1‑150 (MQ=255)
                                cTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTaa                                                                                                                       >  2:1020907/1‑150 (MQ=255)
                                cTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTaa                                                                                                                       <  1:1016781/150‑1 (MQ=255)
                                  cAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATg                                                                                                                     >  1:1582774/1‑150 (MQ=255)
                                  cAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATg                                                                                                                     <  2:641950/150‑1 (MQ=255)
                                       gTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGc                                                                                                                <  2:626920/150‑1 (MQ=255)
                                       gTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGc                                                                                                                <  1:613823/150‑1 (MQ=255)
                                              tGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGcagacaga                                                                                                         <  2:45062/150‑1 (MQ=255)
                                               ggCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGAt                                                                                                        <  2:1164304/150‑1 (MQ=255)
                                                             tCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACCGATCAGAACGCTCCGGc                                                                                          >  2:1639358/1‑150 (MQ=255)
                                                              cACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCt                                                                                         <  1:203764/150‑1 (MQ=255)
                                                              cACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCt                                                                                         >  2:88475/1‑150 (MQ=255)
                                                              cACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCt                                                                                         >  1:1152195/1‑150 (MQ=255)
                                                              cACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCt                                                                                         >  1:1417883/1‑150 (MQ=255)
                                                                  tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACg                                                                                     <  2:1603950/150‑1 (MQ=255)
                                                                  tGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACg                                                                                     <  2:214191/150‑1 (MQ=255)
                                                                   gCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGc                                                                                    >  2:1288272/1‑150 (MQ=255)
                                                                   gCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGc                                                                                    >  2:918450/1‑150 (MQ=255)
                                                                      cGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTa                                                                                 >  2:972571/1‑150 (MQ=255)
                                                                         tatCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAAc                                                                              <  1:117036/150‑1 (MQ=255)
                                                                          atCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACg                                                                             <  1:977454/150‑1 (MQ=255)
                                                                           tCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGt                                                                            >  1:30186/1‑150 (MQ=255)
                                                                            cAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           >  2:74077/1‑150 (MQ=255)
                                                                            cAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTa                                                                           >  1:925780/1‑150 (MQ=255)
                                                                             aaCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACTCTTAAACGTAc                                                                          >  1:1390347/1‑150 (MQ=255)
                                                                             aaCCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTAc                                                                          >  1:1411338/1‑150 (MQ=255)
                                                                                 aGGCGCGGTTTTAGAAACTCTGTCACAGCACCCTGCATTGTGTCCTGGCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCGGACAGATCAGAACGCTCCGGCGTACGCTTAAGCGGAccgtg                                                                      >  2:1507105/1‑147 (MQ=255)
                                                                                     gcgGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATc                                                                  >  2:1479277/1‑150 (MQ=255)
                                                                                        gTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTg                                                               >  1:168793/1‑150 (MQ=255)
                                                                                        gTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGGTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTg                                                               <  1:1651572/150‑1 (MQ=255)
                                                                                           ttAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCa                                                            >  1:523006/1‑150 (MQ=255)
                                                                                                 aCTCTGTCACAGCCCCCTGCATTGTGGCCTCTCTTTGGTACTACGCTTTCCTTGGAGTACAGCTCGACGCTCAGGTGTTCGTTAATGGCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCa                                                      <  2:704233/150‑1 (MQ=255)
                                                                                                 aCTCTGTCACAGAACCCTGCATTGTGCCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCa                                                      <  1:1426532/150‑1 (MQ=255)
                                                                                                    ctGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACt                                                   >  2:1333347/1‑150 (MQ=255)
                                                                                                      gTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTc                                                 >  1:935396/1‑150 (MQ=255)
                                                                                                        cacaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                               >  2:870541/1‑150 (MQ=255)
                                                                                                        cacaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTcca                                               >  1:427502/1‑150 (MQ=255)
                                                                                                         acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccag                                              <  1:1333347/150‑1 (MQ=255)
                                                                                                         acaGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccag                                              >  1:1114542/1‑150 (MQ=255)
                                                                                                            gAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTccagcca                                           <  2:364182/150‑1 (MQ=255)
                                                                                                               cccTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGt                                        >  1:657382/1‑150 (MQ=255)
                                                                                                                 cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGTTCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTg                                      <  1:936238/150‑1 (MQ=255)
                                                                                                                 cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTg                                      >  2:123050/1‑150 (MQ=255)
                                                                                                                 cTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTg                                      >  2:320845/1‑150 (MQ=255)
                                                                                                                    cATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCt                                   >  1:858623/1‑150 (MQ=255)
                                                                                                                     aTTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCtt                                  <  2:382857/150‑1 (MQ=255)
                                                                                                                        gtgtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTc                               <  2:1155737/150‑1 (MQ=255)
                                                                                                                        gtgtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTc                               >  2:397046/1‑150 (MQ=255)
                                                                                                                        gtgtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTc                               >  2:397041/1‑150 (MQ=255)
                                                                                                                          gtCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCAc                             >  1:1641461/1‑150 (MQ=255)
                                                                                                                             ctctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgct                          <  2:935396/150‑1 (MQ=255)
                                                                                                                              tctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctg                         >  1:527847/1‑150 (MQ=255)
                                                                                                                               ctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgc                        <  1:216875/150‑1 (MQ=255)
                                                                                                                               ctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgc                        <  2:1281989/150‑1 (MQ=255)
                                                                                                                               ctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                           >  1:254570/1‑147 (MQ=255)
                                                                                                                               ctctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgc                           <  2:254570/147‑1 (MQ=255)
                                                                                                                                tctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACgctgct                       >  1:1119554/1‑150 (MQ=255)
                                                                                                                                 ctTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTc                      <  1:1164491/150‑1 (MQ=255)
                                                                                                                                     ggTACTAAGCTTTACTTGGAGTAAAGCTCGCCGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagcc                  >  1:448743/1‑150 (MQ=255)
                                                                                                                                     ggTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagcc                  <  1:1456057/150‑1 (MQ=255)
                                                                                                                                       cgtcgaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccag                >  1:1267082/6‑150 (MQ=255)
                                                                                                                                        gtcgaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc               >  2:1181968/5‑150 (MQ=255)
                                                                                                                                        aCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc               <  1:555747/150‑1 (MQ=255)
                                                                                                                                        aCTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTcagccagc               >  2:883375/1‑150 (MQ=255)
                                                                                                                                         tcgaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCt              >  1:110972/4‑150 (MQ=255)
                                                                                                                                         cTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCt              >  1:1453527/1‑150 (MQ=255)
                                                                                                                                          cgaGCTTTACTCCAAGTAAAGCTTGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTc             >  1:923667/3‑150 (MQ=255)
                                                                                                                                           gcgCTTTACTCAAAGTAAACCTCGACGAACAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACCCTCCGGATAACCCTTAAACGTACCATCCATCTTGCCAGCATCAACTTCCAGCCAGCTTGGATATTCACGCTGCTcagccacatcc            >  1:1330318/3‑145 (MQ=255)
                                                                                                                                           gaGCTTTACTCCAAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTcc            >  2:491648/2‑150 (MQ=255)
                                                                                                                                            aGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCa           <  2:925780/150‑1 (MQ=255)
                                                                                                                                              cTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCaga         >  2:1197446/1‑150 (MQ=255)
                                                                                                                                                   cTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGc    <  1:1197446/150‑1 (MQ=255)
                                                                                                                                                    ttGGAGTAAAGCTCGACGATCAGGTGTTAGTTAATGTCCGCATACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCt   <  2:1475859/150‑1 (MQ=255)
                                                                                                                                                     tGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCtt  <  1:979603/150‑1 (MQ=255)
                                                                                                                                                     tGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCAGCAGACAGATAAGAACGCTCCGGCTTACGCTTAAACGTACCTTCGATCTTGCCAGCATCAACTTACAGCCAGGTTCGCTTTTCACGCTGCTCATCCAGCTCAAGAGCGGCtt  >  2:1045422/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGTGTTCGTTAATGTCCGCAGACAGATCAGAACGCTCCGGCTTACGCTTAAACGTACCTTCCATCTTGCCAGCATCAACTTCCAGCCAGGTTGGCTTTTCACGCTGCTCAGCCAGCTCCAGAGCGGCTT  >  NC_000913/3440911‑3441209

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