Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,974,861 G→A 5.5% R132H (CGC→CAC)  wecD → TDP‑fucosamine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,974,8610GA 5.5% 488.4 / 13.4 146R132H (CGC→CAC) wecDTDP‑fucosamine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (67/71);  new base A (4/4);  total (71/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 = 3.13e-01

Only 100 of 149 total aligned reads displayed.

GCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  NC_000913/3974714‑3975004
                                                                                                                                                   |                                                                                                                                               
gCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTgcgc                                                                                                                                                >  4:607191/1‑149 (MQ=255)
gCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTgcgc                                                                                                                                                >  1:83723/1‑149 (MQ=255)
   gTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGcgcgtc                                                                                                                                             <  3:607191/149‑3 (MQ=255)
   gTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGc                                                                                                                                             <  2:83723/149‑1 (MQ=255)
           gTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGa                                                                                                                                     <  3:1125392/149‑1 (MQ=255)
            tGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTTat                                                                                                                                    <  4:573874/149‑1 (MQ=255)
             ggCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatc                                                                                                                                   <  2:272549/149‑1 (MQ=255)
              gCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatca                                                                                                                                  <  1:378901/149‑1 (MQ=255)
              gCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatca                                                                                                                                  >  1:404051/1‑149 (MQ=255)
               cacaAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatcat                                                                                                                                 <  4:1444279/149‑1 (MQ=255)
                 caAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGAtcatca                                                                                                                               <  2:404051/149‑1 (MQ=255)
                       aCCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCt                                                                                                                         >  4:346096/1‑149 (MQ=255)
                       aCCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCt                                                                                                                         <  1:654366/149‑1 (MQ=255)
                         cGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  4:548456/1‑149 (MQ=255)
                          gATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatttt                                                                                                                                        >  2:325257/1‑127 (MQ=255)
                              ttCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       <  2:1359134/144‑1 (MQ=255)
                              ttCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  1:1359134/1‑144 (MQ=255)
                                cccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatt                                                                                                                                          <  1:888387/123‑3 (MQ=255)
                                 ccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCttt                                                                                                                                       <  4:1123744/125‑1 (MQ=255)
                                 ccGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCttt                                                                                                                                       >  3:1123744/1‑125 (MQ=255)
                                   gCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGAtcatca                                                                                                                               <  2:785516/131‑1 (MQ=255)
                                    cATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGatc                                                                                                                                   <  2:698908/126‑1 (MQ=255)
                                    cATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgcg                                                                                                            >  4:525282/1‑149 (MQ=255)
                                       ggcGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  2:841942/3‑135 (MQ=255)
                                          gTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTcc                                                                                                      >  1:269316/1‑149 (MQ=255)
                                           tCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGCGCACGGCatt                                                                                                                                          >  3:775634/1‑110 (MQ=255)
                                               ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTAc                                                                                                                 <  3:1347148/133‑1 (MQ=255)
                                               ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTAc                                                                                                                 >  4:1347148/1‑133 (MQ=255)
                                               ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgc                                                                                                             >  2:958297/1‑137 (MQ=255)
                                               ttAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGAGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTgc                                                                                                             <  1:958297/137‑1 (MQ=255)
                                                          cgcATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTa                                                                                                                  >  2:229927/1‑121 (MQ=255)
                                                          cgcATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTa                                                                                                                  <  1:229927/121‑1 (MQ=255)
                                                                     aaaGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTa                                                                           <  3:548456/149‑1 (MQ=255)
                                                                           gTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGaa                                                                     <  1:1273295/149‑1 (MQ=255)
                                                                            ttttCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGt                                                                                  >  4:889937/1‑135 (MQ=255)
                                                                                       cGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACa                                                         <  3:525282/149‑1 (MQ=255)
                                                                                         tGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTa                                                                           <  1:710411/129‑1 (MQ=255)
                                                                                         tGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTa                                                                           >  2:710411/1‑129 (MQ=255)
                                                                                         tGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGa                                                       <  2:215210/149‑1 (MQ=255)
                                                                                              tGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGAt                                                      >  3:1159721/1‑145 (MQ=255)
                                                                                              tGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGAt                                                      <  4:1159721/145‑1 (MQ=255)
                                                                                                 gcCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGaa                                               >  3:2848/1‑149 (MQ=255)
                                                                                                 gcCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGAGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGaa                                               >  3:1063064/1‑149 (MQ=255)
                                                                                                   cTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAAtt                                             >  1:1012761/1‑149 (MQ=255)
                                                                                                    tGACGCCAGCAGTAGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTg                                            <  2:980527/149‑1 (MQ=255)
                                                                                                      aCGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAAGGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGc                                          <  4:1230268/149‑1 (MQ=255)
                                                                                                        gCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCct                                        >  3:335171/1‑149 (MQ=255)
                                                                                                          cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       >  3:343699/1‑68 (MQ=255)
                                                                                                          cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGa                                                                                                                       <  4:343699/68‑1 (MQ=255)
                                                                                                          cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgc                                      >  4:1107365/1‑149 (MQ=255)
                                                                                                          cagcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgc                                      <  1:804665/149‑1 (MQ=255)
                                                                                                            gcagTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgctg                                    >  4:795196/1‑149 (MQ=255)
                                                                                                                tCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtc                                                                                 >  2:1029977/1‑100 (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)
                                                                                                                 cGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCCCCTATGTCTCTTTACGGGAAc                                                                    <  2:1266938/112‑1 (MQ=255)
                                                                                                                  gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAAt                                                               <  3:810667/116‑1 (MQ=255)
                                                                                                                  gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAAt                                                               >  4:810667/1‑116 (MQ=255)
                                                                                                                  gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGa                              >  1:643694/1‑149 (MQ=255)
                                                                                                                  gCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGa                              <  2:1350762/149‑1 (MQ=255)
                                                                                                                    ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCAc                                                                                                                                           <  4:737267/38‑1 (MQ=255)
                                                                                                                    ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCAc                                                                                                                                           >  3:737267/1‑38 (MQ=255)
                                                                                                                    ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtctc                                                                               <  1:1087749/98‑1 (MQ=255)
                                                                                                                    ttttATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGtctc                                                                               >  2:1087749/1‑98 (MQ=255)
                                                                                                                       tATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCgg                         <  2:451165/149‑1 (MQ=255)
                                                                                                                       tATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCgg                         <  3:1107365/149‑1 (MQ=255)
                                                                                                                        aTGCACAGTGGATGGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGt                        >  4:796404/1‑149 (MQ=255)
                                                                                                                         tGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTg                       <  3:1373362/149‑1 (MQ=255)
                                                                                                                          gCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGaa                                                                     >  2:1122412/1‑102 (MQ=255)
                                                                                                                          gCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGaa                                                                     <  1:1122412/102‑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                                                   >  1:1421768/1‑115 (MQ=255)
                                                                                                                                tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTgctggct                                 <  4:236991/132‑1 (MQ=255)
                                                                                                                                tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGc                      >  4:1062346/1‑143 (MQ=255)
                                                                                                                                tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGc                      <  3:1062346/143‑1 (MQ=255)
                                                                                                                                tGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGAGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTgctggct                                 >  3:236991/1‑132 (MQ=255)
                                                                                                                                  gATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctg              <  4:1063064/149‑1 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATg                                                              <  4:637196/99‑1 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGAc                             >  2:1035557/1‑132 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGAc                             <  1:1035557/132‑1 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgag            >  2:568955/1‑149 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTATACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgag            >  3:965102/1‑149 (MQ=255)
                                                                                                                                    ttGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATAGTCGCGGCTATGTCTCTTTACGGGAACTCAATg                                                              >  3:637196/1‑99 (MQ=255)
                                                                                                                                      gAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgagct          >  3:446286/1‑149 (MQ=255)
                                                                                                                                       aaggtGCCGCGCACGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgc                >  2:524912/5‑142 (MQ=255)
                                                                                                                                       aaggtGCCGCGCACGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgc                <  1:524912/138‑1 (MQ=255)
                                                                                                                                       aaaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACg                                                                         <  4:520340/85‑1 (MQ=255)
                                                                                                                                       aaaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCctgct                                     <  3:955117/121‑1 (MQ=255)
                                                                                                                                        aaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgagctga        >  3:1438914/1‑149 (MQ=255)
                                                                                                                                        aaaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTgctgagctga        >  4:742454/1‑149 (MQ=255)
                                                                                                                                         aaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGAt       <  4:2848/149‑1 (MQ=255)
                                                                                                                                         aaTGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGAt       >  1:511860/1‑149 (MQ=255)
                                                                                                                                              cGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCaaa  <  2:643694/149‑1 (MQ=255)
                                                                                                                                               gTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCct                                        >  3:1149077/1‑110 (MQ=255)
                                                                                                                                               gTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCct                                        <  4:1149077/110‑1 (MQ=255)
                                                                                                                                                tGCGCGGCGCCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGc     <  1:70167/144‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCGGTGCTGTAGTGGCACAAGAGACCGATATTCCCGCATTACGTCAGTTAGCCAGCGCCGCATTTGCGCAAAGCCGTTTTCGTGCGCCGTGGTATGCGCCTGACGCCAGCAGTCGCTTTTATGCACAGTGGATTGAAAATGCCGTGCGCGGCACCTTTGATCATCAATGTCTGATTTTACGTGCGGCGTCCGGCGATATTCGCGGCTATGTCTCTTTACGGGAACTCAATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAA  >  NC_000913/3974714‑3975004

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