Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,974,858 T→C 5.4% V131A (GTG→GCG)  wecD → TDP‑fucosamine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,974,8580TC 5.4% 507.3 / 12.7 147V131A (GTG→GCG) wecDTDP‑fucosamine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (68/71);  new base C (4/4);  total (72/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 150 total aligned reads displayed.

GATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  NC_000913/3974710‑3975004
                                                                                                                                                    |                                                                                                                                                  
gATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGt                                                                                                                                                    <  3:45077/149‑1 (MQ=255)
    gCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTgcgc                                                                                                                                                >  4:607191/1‑149 (MQ=255)
      ggTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCgg                                                                                                                                              >  3:1379864/1‑149 (MQ=255)
       gTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGcgcgtc                                                                                                                                             <  3:607191/149‑3 (MQ=255)
       gTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGc                                                                                                                                             <  2:83723/149‑1 (MQ=255)
             tAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCttt                                                                                                                                       >  2:378901/1‑149 (MQ=255)
               gTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGa                                                                                                                                     <  3:1125392/149‑1 (MQ=255)
                tGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTTat                                                                                                                                    <  4:573874/149‑1 (MQ=255)
                tGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGat                                                                                                                                    >  1:980527/1‑149 (MQ=255)
                 ggCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatc                                                                                                                                   <  2:272549/149‑1 (MQ=255)
                  gCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatca                                                                                                                                  <  1:378901/149‑1 (MQ=255)
                    acaAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCt                                                                                                                                         <  3:1125445/140‑1 (MQ=255)
                    acaAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatcatc                                                                                                                                >  2:1273295/1‑149 (MQ=255)
                     caAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGAtcatca                                                                                                                               <  2:404051/149‑1 (MQ=255)
                      aaGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCaa                                                                                                                              <  4:1379864/149‑1 (MQ=255)
                       agagACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAAt                                                                                                                             <  3:218406/149‑1 (MQ=255)
                           aCCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCt                                                                                                                         >  4:346096/1‑149 (MQ=255)
                             cGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  4:548456/1‑149 (MQ=255)
                              gATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatttt                                                                                                                                        <  1:325257/131‑5 (MQ=255)
                                  ttCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  1:1359134/1‑144 (MQ=255)
                                  ttCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTa                                                                                                                  <  3:346096/149‑1 (MQ=255)
                                    cccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatt                                                                                                                                          >  2:888387/1‑121 (MQ=255)
                                    cccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatt                                                                                                                                          <  1:888387/123‑3 (MQ=255)
                                     ccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCttt                                                                                                                                       <  4:1123744/125‑1 (MQ=255)
                                       gCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGAtcatca                                                                                                                               >  1:785516/1‑131 (MQ=255)
                                        cATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatc                                                                                                                                   >  1:698908/1‑126 (MQ=255)
                                        cATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgcg                                                                                                            >  4:525282/1‑149 (MQ=255)
                                           ggcGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       <  1:841942/133‑1 (MQ=255)
                                               tCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatt                                                                                                                                          >  3:775634/1‑110 (MQ=255)
                                                  gTTAGCCAGCGCCGGATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCg                                                                                                  >  3:108455/1‑149 (MQ=255)
                                                   ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTAc                                                                                                                 >  4:1347148/1‑133 (MQ=255)
                                                   ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgc                                                                                                             >  2:958297/1‑137 (MQ=255)
                                                            gccgcATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCgg                                                                                        >  3:1038738/1‑149 (MQ=255)
                                                             ccgcATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGc                                                                                       >  1:215210/1‑149 (MQ=255)
                                                              cgcATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTa                                                                                                                  <  1:229927/121‑1 (MQ=255)
                                                                  tttGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTAc                                                                                                                 <  3:1169943/118‑1 (MQ=255)
                                                                       gcAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCtt                                                                             >  3:323530/1‑149 (MQ=255)
                                                                        cAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGTGATATTCGCGGCTATGTCTCttt                                                                            <  2:269316/149‑1 (MQ=255)
                                                                         aaaGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTa                                                                           <  3:548456/149‑1 (MQ=255)
                                                                                ttttCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGt                                                                                  >  4:889937/1‑135 (MQ=255)
                                                                                           cGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACa                                                         <  3:525282/149‑1 (MQ=255)
                                                                                             tGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTa                                                                           >  2:710411/1‑129 (MQ=255)
                                                                                             tGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGa                                                       <  2:215210/149‑1 (MQ=255)
                                                                                                  tGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGAt                                                      <  4:1159721/145‑1 (MQ=255)
                                                                                                  tGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGAt                                                      >  3:1159721/1‑145 (MQ=255)
                                                                                                     gcCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACt                                                                   >  2:1193236/1‑129 (MQ=255)
                                                                                                     gcCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGaa                                               >  2:804665/1‑149 (MQ=255)
                                                                                                     gcCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGaa                                               >  3:2848/1‑149 (MQ=255)
                                                                                                     gcCTGACGCCAGCACTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGACCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACt                                                                   <  1:1193236/129‑1 (MQ=255)
                                                                                                       cTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAAtt                                             >  1:1012761/1‑149 (MQ=255)
                                                                                                          aCGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAAGGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGc                                          <  4:1230268/149‑1 (MQ=255)
                                                                                                            gCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCct                                        >  1:1350762/1‑149 (MQ=255)
                                                                                                              cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       <  4:343699/68‑1 (MQ=255)
                                                                                                              cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  3:343699/1‑68 (MQ=255)
                                                                                                              cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgc                                      <  1:804665/149‑1 (MQ=255)
                                                                                                              cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgc                                      >  4:1107365/1‑149 (MQ=255)
                                                                                                                gcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgctg                                    >  4:795196/1‑149 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtc                                                                                 >  2:605380/1‑100 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtc                                                                                 >  2:1029977/1‑100 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtc                                                                                 <  1:605380/100‑1 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTAc                                                                          >  2:1041358/1‑107 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTAc                                                                          <  1:1041358/107‑1 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATgc                                                    <  3:489576/129‑1 (MQ=255)
                                                                                                                    tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGAACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtc                                                                                 <  1:1029977/100‑1 (MQ=255)
                                                                                                                     cGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAAc                                                                    >  1:1266938/1‑112 (MQ=255)
                                                                                                                     cGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGctggctgg                               >  2:1238499/1‑149 (MQ=255)
                                                                                                                      gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAAt                                                               >  4:810667/1‑116 (MQ=255)
                                                                                                                      gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAAt                                                               <  3:810667/116‑1 (MQ=255)
                                                                                                                      gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGa                              >  1:451165/1‑149 (MQ=255)
                                                                                                                      gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGa                              <  2:1350762/149‑1 (MQ=255)
                                                                                                                      gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGa                              >  1:643694/1‑149 (MQ=255)
                                                                                                                        ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCAc                                                                                                                                           >  3:737267/1‑38 (MQ=255)
                                                                                                                        ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCAc                                                                                                                                           <  4:737267/38‑1 (MQ=255)
                                                                                                                        ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtctc                                                                               <  1:1087749/98‑1 (MQ=255)
                                                                                                                        ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtctc                                                                               >  2:1087749/1‑98 (MQ=255)
                                                                                                                           tATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCgg                         <  3:1107365/149‑1 (MQ=255)
                                                                                                                            aTGCACAGTGGATGGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGt                        >  4:796404/1‑149 (MQ=255)
                                                                                                                             tGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTg                       <  3:1373362/149‑1 (MQ=255)
                                                                                                                                acaGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAg                    <  3:796404/149‑1 (MQ=255)
                                                                                                                                acaGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAg                    <  4:1461225/149‑1 (MQ=255)
                                                                                                                                  aGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGt                  <  3:210770/149‑1 (MQ=255)
                                                                                                                                   gTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATgcg                                                   <  2:1421768/115‑1 (MQ=255)
                                                                                                                                    tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTgctggct                                 <  4:236991/132‑1 (MQ=255)
                                                                                                                                    tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGc                      >  4:1062346/1‑143 (MQ=255)
                                                                                                                                    tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGc                      <  3:1062346/143‑1 (MQ=255)
                                                                                                                                      gATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctg              <  4:1063064/149‑1 (MQ=255)
                                                                                                                                        ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATg                                                              <  4:637196/99‑1 (MQ=255)
                                                                                                                                        ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTATACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgag            >  3:965102/1‑149 (MQ=255)
                                                                                                                                        ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATAGTCGCGGCTATGTCTCTTTACGGGAACTCAATg                                                              >  3:637196/1‑99 (MQ=255)
                                                                                                                                           aaggtGCCGCGCACGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgc                >  2:524912/5‑142 (MQ=255)
                                                                                                                                           aaaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACg                                                                         <  4:520340/85‑1 (MQ=255)
                                                                                                                                           aaaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgct                                     <  3:955117/121‑1 (MQ=255)
                                                                                                                                            aaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgagctga        >  4:742454/1‑149 (MQ=255)
                                                                                                                                             aaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGAt       <  4:2848/149‑1 (MQ=255)
                                                                                                                                             aaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGAt       >  1:511860/1‑149 (MQ=255)
                                                                                                                                             aaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGAt       <  2:511860/149‑1 (MQ=255)
                                                                                                                                                  cGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCaaa  <  2:643694/149‑1 (MQ=255)
                                                                                                                                                   gTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCct                                        <  4:1149077/110‑1 (MQ=255)
                                                                                                                                                    tGCGCGGCGCCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGc     <  1:70167/144‑1 (MQ=255)
                                                                                                                                                    tGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGc     >  2:70167/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATAGCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  NC_000913/3974710‑3975004

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: