Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,971,621 C→T G134G (GGC→GGT rffH → glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,6210CT89.4% 851.2 / 67.5 282G134G (GGC→GGTrffHglucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (11/19);  new base T (116/136);  total (127/155)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-02

Only 100 of 283 total aligned reads displayed.

GGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCG  >  NC_000913/3971531‑3971711
                                                                                          |                                                                                          
ggTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTa                                                                                          >  1:5575488/1‑93 (MQ=255)
ggTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTa                                                                                          <  2:6055760/93‑1 (MQ=255)
  tGCTGGGCGATAACATCTTCTTCGGTCAGGCCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATc                                                                                        <  1:5885757/93‑1 (MQ=255)
   gCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCa                                                                                       <  2:1704785/93‑1 (MQ=255)
   gCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCa                                                                                       <  2:3738707/93‑1 (MQ=255)
   gCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCa                                                                                       <  1:419080/93‑1 (MQ=255)
    cTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAg                                                                                      >  2:6338325/1‑93 (MQ=255)
    cTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAg                                                                                      <  1:6522343/93‑1 (MQ=255)
      gggCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGt                                                                                    >  1:6231143/1‑93 (MQ=255)
      gggCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGt                                                                                    <  1:5890668/93‑1 (MQ=255)
      gggCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGt                                                                                    <  2:5818666/93‑1 (MQ=255)
       ggCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTg                                                                                   >  2:6242837/1‑93 (MQ=255)
        gCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGa                                                                                  <  2:2980056/93‑1 (MQ=255)
        gCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGa                                                                                  <  2:3106233/93‑1 (MQ=255)
        gCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGa                                                                                  <  1:3589247/93‑1 (MQ=255)
         cGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGAt                                                                                 >  1:309422/1‑93 (MQ=255)
         cGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGAt                                                                                 <  2:1103041/93‑1 (MQ=255)
          gATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATg                                                                                >  1:2960082/1‑93 (MQ=255)
          gATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATg                                                                                <  2:885105/93‑1 (MQ=255)
          gATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATg                                                                                >  1:2220226/1‑93 (MQ=255)
           aTAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATgg                                                                               <  1:6563591/93‑1 (MQ=255)
           aTAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATgg                                                                               <  1:233934/93‑1 (MQ=255)
            tAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGa                                                                              >  1:5167056/1‑93 (MQ=255)
             aaCATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAc                                                                             >  1:3414115/1‑93 (MQ=255)
             aaCATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAc                                                                             >  1:6091760/1‑93 (MQ=255)
              aCATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAcc                                                                            <  1:4857293/93‑1 (MQ=255)
              aCATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAcc                                                                            >  1:3212998/1‑93 (MQ=255)
               cATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAccc                                                                           <  1:6480248/93‑1 (MQ=255)
               cATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAccc                                                                           <  2:4111618/93‑1 (MQ=255)
               cATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGAccc                                                                           >  2:1235792/1‑93 (MQ=255)
                aTCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCg                                                                          >  1:45934/1‑93 (MQ=255)
                  cttcttCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGa                                                                        <  1:5228460/93‑1 (MQ=255)
                   ttcttcGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGaa                                                                       >  2:4191285/1‑93 (MQ=255)
                   ttcttcGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGaa                                                                       <  1:5025366/93‑1 (MQ=255)
                   ttcttcGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGaa                                                                       >  2:242770/1‑93 (MQ=255)
                   ttcttcGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGaa                                                                       <  1:3938491/93‑1 (MQ=255)
                      ttcGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGCACGc                                                                    >  1:4764463/1‑93 (MQ=255)
                        cGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCtt                                                                  <  2:329939/93‑1 (MQ=255)
                         ggTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCttt                                                                 <  2:2849235/93‑1 (MQ=255)
                         ggTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCttt                                                                 <  1:334170/93‑1 (MQ=255)
                         ggTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCttt                                                                 >  2:527585/1‑93 (MQ=255)
                         ggTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCttt                                                                 <  2:5660691/93‑1 (MQ=255)
                         ggTCAGGGCTTCAGTCAGAAGCTGCGTGATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCttt                                                                 <  1:149192/93‑1 (MQ=255)
                          gTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTg                                                                >  2:6507131/1‑93 (MQ=255)
                           tCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTgg                                                               <  1:547949/93‑1 (MQ=255)
                            cAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGc                                                              <  2:1687064/93‑1 (MQ=255)
                            cAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGc                                                              >  1:4354600/1‑93 (MQ=255)
                             aGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCg                                                             <  2:1875277/93‑1 (MQ=255)
                              gggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgt                                                            >  1:5703400/1‑93 (MQ=255)
                              gggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgt                                                            >  1:4108935/1‑93 (MQ=255)
                              gggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgt                                                            <  1:4063288/93‑1 (MQ=255)
                              gggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgg                                                            >  2:3209033/1‑92 (MQ=255)
                               ggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgtg                                                           >  2:1927790/1‑93 (MQ=255)
                               ggCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCgtg                                                           <  2:3570040/93‑1 (MQ=255)
                                  ttCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCgtggtg                                                        <  1:1011019/93‑1 (MQ=255)
                                   tCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGtggtgg                                                       <  2:4857489/93‑1 (MQ=255)
                                    cAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGa                                                      >  1:2488809/1‑93 (MQ=255)
                                     aGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAg                                                     <  2:6278732/93‑1 (MQ=255)
                                      gTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGt                                                    <  1:1207839/93‑1 (MQ=255)
                                      gTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGt                                                    >  1:1724352/1‑93 (MQ=255)
                                      gTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGt                                                    <  2:1565366/93‑1 (MQ=255)
                                         cGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTg                                                 >  2:3013516/1‑93 (MQ=255)
                                         cGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGCGTTTg                                                 >  1:1369609/1‑93 (MQ=255)
                                             gCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTgacga                                             >  2:4496283/1‑93 (MQ=255)
                                               tGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACa                                           >  2:4372446/1‑93 (MQ=255)
                                               tGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACa                                           >  1:2615831/1‑93 (MQ=255)
                                               tGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACa                                           >  2:4694310/1‑93 (MQ=255)
                                               tGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACa                                           <  1:1709051/93‑1 (MQ=255)
                                                gCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACaa                                          <  1:5499299/93‑1 (MQ=255)
                                                gCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACaa                                          >  2:279514/1‑93 (MQ=255)
                                                 cGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAAt                                         <  2:5070938/93‑1 (MQ=255)
                                                  gTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAAtt                                        <  1:816192/93‑1 (MQ=255)
                                                   tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAAttt                                       <  1:4697533/93‑1 (MQ=255)
                                                    cATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTc                                      <  2:1598543/93‑1 (MQ=255)
                                                    cATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTc                                      <  2:180514/93‑1 (MQ=255)
                                                     atattGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTcc                                     <  2:4915791/90‑1 (MQ=255)
                                                     aTGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTcc                                     <  2:499269/93‑1 (MQ=255)
                                                          gcggcgCGCACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCt                                >  2:4043759/1‑93 (MQ=255)
                                                             gcgcgcACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTAtc                             >  1:2859701/1‑93 (MQ=255)
                                                               gcgcACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTAtctc                           <  2:4541083/93‑1 (MQ=255)
                                                                cgcACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGGTTGACGACAATTTCCGCGCTATCTCg                          >  2:133886/1‑93 (MQ=255)
                                                                  cACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCt                        <  2:5560/93‑1 (MQ=255)
                                                                  cACCGAAGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCt                        <  1:2442026/93‑1 (MQ=255)
                                                                       aaGGGGCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGgaag                   >  1:4256302/1‑93 (MQ=255)
                                                                           ggCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGaaaa               <  2:258432/93‑1 (MQ=255)
                                                                            gCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAg              >  2:3397282/1‑93 (MQ=255)
                                                                            gCGACGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAg              <  2:677596/93‑1 (MQ=255)
                                                                              gACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGcc            >  1:5784640/1‑93 (MQ=255)
                                                                                cGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCaa          >  2:4037376/1‑93 (MQ=255)
                                                                                cGGTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCaa          >  1:5938598/1‑93 (MQ=255)
                                                                                  gTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCaaaa        >  1:4732601/1‑93 (MQ=255)
                                                                                  gTTTTTGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCaaaa        <  2:5774857/93‑1 (MQ=255)
                                                                                   tttttGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAAc       <  1:1871959/93‑1 (MQ=255)
                                                                                   tttttGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAAc       <  1:1571747/93‑1 (MQ=255)
                                                                                     tttGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAg     <  1:5812097/93‑1 (MQ=255)
                                                                                      ttGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGc    <  2:2104490/93‑1 (MQ=255)
                                                                                      ttGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGc    <  2:5498382/93‑1 (MQ=255)
                                                                                      ttGGCTATCAGGTGATGGACCCGGCACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGc    >  1:1801594/1‑93 (MQ=255)
                                                                                      ttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGc    <  2:2032002/93‑1 (MQ=255)
                                                                                       tGGTTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGcc   >  1:6413768/1‑93 (MQ=255)
                                                                                          |                                                                                          
GGTGCTGGGCGATAACATCTTCTTCGGTCAGGGCTTCAGTCCGAAGCTGCGTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCG  >  NC_000913/3971531‑3971711

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