Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 720,070 G→C 6.7% H131D (CAT→GAT) ‡ speF ← ornithine decarboxylase, degradative

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913720,0700GC 6.7% 306.7 / 10.0 90H131D (CAT→GAT) ‡speFornithine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base G (36/48);  new base C (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 = 9.25e-01

CTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCAAAG  >  NC_000913/719921‑720218
                                                                                                                                                     |                                                                                                                                                    
cTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCtgatg                                                                                                                                                      <  1:1522676/150‑1 (MQ=255)
  gtgCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGAc                                                                                                                                                    >  2:1054257/1‑150 (MQ=255)
  gtgCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGAc                                                                                                                                                    >  2:1106709/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                                                                                                        >  1:178767/1‑150 (MQ=255)
                                              gCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGAcc                                                                                                        >  2:1988072/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:1850161/1‑150 (MQ=255)
                                                                      gCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGgcgc                                                                                >  2:1857021/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                                                                    <  1:1791111/150‑1 (MQ=255)
                                                                                  gcCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCa                                                                    <  2:974118/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:1850161/150‑1 (MQ=255)
                                                                                                      aaaCTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTc                                                <  1:1857021/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:70893/1‑150 (MQ=255)
                                                                                                                                             gccCTGGTCATCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAAtt         >  1:830527/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)
                                                                                                                                                    tgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCAAAg  >  2:136251/1‑150 (MQ=255)
                                                                                                                                                    tgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCAAAg  <  1:49671/150‑1 (MQ=255)
                                                                                                                                                    tgaCCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCAAAg  <  1:1795144/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGTGCAATGCATGGCGCGCCTTCGTGAATCAGCAGATCGCCCATCGCTACGTCGGCGTTGCACAAGTCGGCACGGAACAGCGCCTCACCAAAGTATTCCACAAACTGATTCCCCGCCGGATGGCGACGGAAAAATTCGCCGCCCTGATGACCAGGGCAATCAAACGCGCTGTTACCTTGATTGACATAATCGACCAGTGCGCGGAAGAAAGGTGGGCGCAGTTGAGTTTCATAATGGCTGGCAGCGGTTTCTAACTGGCGACCATAAAATTCTCGTCGCGATTCGCAATTCTCAAAG  >  NC_000913/719921‑720218

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: