Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 720,068 A→G 6.7% H131H (CAT→CAC) ‡ speF ← ornithine decarboxylase, degradative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913720,0680AG 6.7% 267.5 / 7.5 90H131H (CAT→CAC) ‡speFornithine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base A (36/48);  new base G (3/3);  total (39/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCA  >  NC_000913/719919‑720215
                                                                                                                                                     |                                                                                                                                                   
tGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCtga                                                                                                                                                     <  2:999969/150‑1 (MQ=255)
tGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCtga                                                                                                                                                     >  1:97807/1‑150 (MQ=255)
 gCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCtgat                                                                                                                                                    <  1:1811912/150‑1 (MQ=255)
  cTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCtgatg                                                                                                                                                   <  1:1522676/150‑1 (MQ=255)
    gtgCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGAc                                                                                                                                                 >  2:1106709/1‑150 (MQ=255)
    gtgCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGAc                                                                                                                                                 >  2:1054257/1‑150 (MQ=255)
      gCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCa                                                                                                                                               >  1:1474870/1‑150 (MQ=255)
       cAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAg                                                                                                                                              >  1:959919/1‑150 (MQ=255)
        aaTGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAgg                                                                                                                                             <  1:1704678/150‑1 (MQ=255)
        aaTGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAgg                                                                                                                                             >  2:1128986/1‑150 (MQ=255)
          tGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGc                                                                                                                                           <  2:1007434/150‑1 (MQ=255)
          tGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGc                                                                                                                                           <  2:991373/150‑1 (MQ=255)
           gCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCa                                                                                                                                          >  2:765902/1‑150 (MQ=255)
            cATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCaa                                                                                                                                         <  2:363026/150‑1 (MQ=255)
            cATGGCGCGCCTTCGTGAATCAGCAGATCGCCAATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCaa                                                                                                                                         >  2:1278353/1‑150 (MQ=255)
             aTGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGGTCATCAGGGcggc                                                                                                                                        <  2:1941345/150‑4 (MQ=255)
               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGGTCATCAGGGcggcga                                                                                                                                      <  1:496111/150‑6 (MQ=255)
               ggCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGGTCATCAGGGcggcga                                                                                                                                      <  1:1588064/150‑6 (MQ=255)
                   cgcCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAAcg                                                                                                                                  >  2:1048080/1‑150 (MQ=255)
                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAAcgcg                                                                                                                                <  2:1384307/150‑1 (MQ=255)
                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAAcgcg                                                                                                                                <  2:256966/150‑1 (MQ=255)
                     ccTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAAcgcg                                                                                                                                >  1:1036417/1‑150 (MQ=255)
                         cGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATc                                                                                                                                       <  1:1795993/139‑1 (MQ=255)
                         cGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATc                                                                                                                                       >  2:1795993/1‑139 (MQ=255)
                           tGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTa                                                                                                                          <  1:1101203/150‑1 (MQ=255)
                           tGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTa                                                                                                                          >  2:188751/1‑150 (MQ=255)
                            gAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTAc                                                                                                                         <  1:1054257/150‑1 (MQ=255)
                              aTCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCt                                                                                                                       >  2:921627/1‑150 (MQ=255)
                               tCAGCAGATCGCCCATCGCTATGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCtt                                                                                                                      <  1:1110441/150‑1 (MQ=255)
                                cagcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttg                                                                                                                     >  1:981706/1‑150 (MQ=255)
                                 agcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttga                                                                                                                    >  2:1818381/1‑150 (MQ=255)
                                 agcagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttga                                                                                                                    >  2:1791111/1‑150 (MQ=255)
                                   cagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttgatt                                                                                                                  <  1:1048080/150‑1 (MQ=255)
                                   cagATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCttgatt                                                                                                                  <  1:765902/150‑1 (MQ=255)
                                      aTCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGAc                                                                                                               <  1:534399/150‑1 (MQ=255)
                                       tCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACa                                                                                                              <  1:1193593/150‑1 (MQ=255)
                                          cccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATaa                                                                                                           >  1:1879363/1‑150 (MQ=255)
                                          cccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATaa                                                                                                           >  2:1498380/1‑150 (MQ=255)
                                           ccATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAAt                                                                                                          <  1:1002771/150‑1 (MQ=255)
                                             aTCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCg                                                                                                        <  1:363489/150‑1 (MQ=255)
                                             aTCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCg                                                                                                        >  1:1114316/1‑150 (MQ=255)
                                                gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAcc                                                                                                     >  2:1988072/1‑150 (MQ=255)
                                                gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAcc                                                                                                     >  1:178767/1‑150 (MQ=255)
                                                gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAcc                                                                                                     <  2:596981/150‑1 (MQ=255)
                                                       cGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTgcgc                                                                                              <  2:1202392/150‑1 (MQ=255)
                                                         gCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCgg                                                                                            <  2:409997/150‑1 (MQ=255)
                                                              gCACAAGTCGGCACGGAACAGCGCCTCCCCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGgaagaa                                                                                       >  1:1665436/1‑150 (MQ=255)
                                                                   aGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTg                                                                                  <  2:1879363/150‑1 (MQ=255)
                                                                      cGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGgc                                                                               >  1:225506/1‑150 (MQ=255)
                                                                       ggCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGgcg                                                                              <  1:1106709/150‑1 (MQ=255)
                                                                        gCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGgcgc                                                                             >  2:1857021/1‑150 (MQ=255)
                                                                        gCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGgcgc                                                                             >  2:1850161/1‑150 (MQ=255)
                                                                         cACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCa                                                                            >  2:1215205/1‑150 (MQ=255)
                                                                          aCGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAg                                                                           >  1:1989148/1‑150 (MQ=255)
                                                                            ggAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGtt                                                                         <  2:225506/150‑1 (MQ=255)
                                                                             gAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTg                                                                        <  1:1497000/150‑1 (MQ=255)
                                                                              aaCAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGa                                                                       <  2:1319186/150‑1 (MQ=255)
                                                                                  gcgcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGttt                                                                   <  1:1044286/150‑1 (MQ=255)
                                                                                  gcgcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGttt                                                                   >  1:995791/1‑150 (MQ=255)
                                                                                    gcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCa                                                                 <  2:974118/150‑1 (MQ=255)
                                                                                    gcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCa                                                                 <  1:1791111/150‑1 (MQ=255)
                                                                                      cTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATa                                                               <  2:97807/150‑1 (MQ=255)
                                                                                       tCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATaa                                                              >  2:228532/1‑150 (MQ=255)
                                                                                        cACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAt                                                             <  1:18038/150‑1 (MQ=255)
                                                                                        cACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAt                                                             >  2:292739/1‑150 (MQ=255)
                                                                                         aCCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGga                                                                                           >  1:813862/1‑119 (MQ=255)
                                                                                         aCCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGga                                                                                           <  2:813862/119‑1 (MQ=255)
                                                                                           cAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggc                                                          <  2:1114316/150‑1 (MQ=255)
                                                                                             aaGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctg                                                        <  1:947786/150‑1 (MQ=255)
                                                                                             aaGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctg                                                        <  1:940687/150‑1 (MQ=255)
                                                                                               gTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctggc                                                      <  2:1112900/150‑1 (MQ=255)
                                                                                               gTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAAtggctggc                                                      <  2:1119962/150‑1 (MQ=255)
                                                                                                   tCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCg                                                  >  2:793323/1‑150 (MQ=255)
                                                                                                        aaaCTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTc                                             <  1:1857021/150‑1 (MQ=255)
                                                                                                        aaaCTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTc                                             <  1:1850161/150‑1 (MQ=255)
                                                                                                             gATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACt                                        <  1:228532/150‑1 (MQ=255)
                                                                                                                 ccccGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCg                                    <  1:1278353/150‑1 (MQ=255)
                                                                                                                  cccGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGa                                   >  2:1103855/1‑150 (MQ=255)
                                                                                                                  cccGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGa                                   <  2:959919/150‑1 (MQ=255)
                                                                                                                          aTGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATaaaa                           <  2:981706/150‑1 (MQ=255)
                                                                                                                          aTGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATaaaa                           >  2:1081629/1‑150 (MQ=255)
                                                                                                                                  ggAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCtcgtc                   <  2:1036417/150‑1 (MQ=255)
                                                                                                                                    aaaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCTGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTcgc                 <  2:1539492/150‑1 (MQ=255)
                                                                                                                                    aaaaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTcgc                 <  2:1474870/150‑1 (MQ=255)
                                                                                                                                      aaaTTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGa               >  1:1387744/1‑150 (MQ=255)
                                                                                                                                         ttCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTc            <  1:1988072/150‑1 (MQ=255)
                                                                                                                                         ttCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTc            >  1:1793377/1‑150 (MQ=255)
                                                                                                                                          tgattgccCTGGTCATCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCg           >  2:125128/6‑150 (MQ=255)
                                                                                                                                          tCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCg           <  1:13395/150‑1 (MQ=255)
                                                                                                                                           cgccgcCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATATATTTCTCGTCGCGATTCGc          <  2:1665436/150‑1 (MQ=255)
                                                                                                                                               gccCTGGTCATCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAAtt      >  1:830527/1‑150 (MQ=255)
                                                                                                                                               gccCTGGTCATCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAAtt      >  1:70893/1‑150 (MQ=255)
                                                                                                                                               gccCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAAtt      >  2:787079/1‑150 (MQ=255)
                                                                                                                                                 ccTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATtct    <  2:755924/150‑1 (MQ=255)
                                                                                                                                                   tgatgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCa  >  1:1126002/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCA  >  NC_000913/719919‑720215

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