Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,903,704 C→T intergenic (+330/+1043) ygcE → / ← queE putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,903,7040CT99.6% 1122.9 / ‑4.8 273intergenic (+330/+1043)ygcE/queEputative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (124/148);  total (124/149)
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.83e-01

Only 100 of 273 total aligned reads displayed.

TTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTA  >  NC_000913/2903606‑2903801
                                                                                                  |                                                                                                 
ttATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTc                                                                                                 <  2:3755403/101‑1 (MQ=255)
   tCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCaaa                                                                                              >  2:3755927/1‑101 (MQ=255)
     aaTTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAAc                                                                                            >  1:3755929/1‑101 (MQ=255)
        tttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACaaa                                                                                         <  1:3755530/101‑1 (MQ=255)
         ttCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAg                                                                                        <  2:3755774/101‑1 (MQ=255)
            aCTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACa                                                                                     <  1:3755505/101‑1 (MQ=255)
             cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACaa                                                                                    >  1:3755935/1‑101 (MQ=255)
             cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACaa                                                                                    <  1:3755668/101‑1 (MQ=255)
             cTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACaa                                                                                    >  2:3755934/1‑101 (MQ=255)
               tAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAAc                                                                                  >  1:3755939/1‑101 (MQ=255)
                aaaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACt                                                                                 <  2:3755439/101‑1 (MQ=255)
                 aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACaaa                                                                                   >  2:3755943/1‑98 (MQ=255)
                 aaGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACtt                                                                                >  1:3755942/1‑101 (MQ=255)
                    tATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCa                                                                             >  2:3755949/1‑101 (MQ=255)
                      tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTcc                                                                              >  1:3755956/1‑98 (MQ=255)
                      tGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGa                                                                           >  2:3755953/1‑101 (MQ=255)
                       gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGAt                                                                          <  2:3755540/101‑1 (MQ=255)
                       gAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGAt                                                                          <  2:3755728/101‑1 (MQ=255)
                        aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGAtt                                                                         >  2:3755958/1‑101 (MQ=255)
                        aaCTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGAtt                                                                         <  2:3755586/101‑1 (MQ=255)
                          cTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATtta                                                                       >  1:3755960/1‑101 (MQ=255)
                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATttatt                                                                     <  1:3755687/101‑1 (MQ=255)
                            ggATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATttatt                                                                     <  2:3755341/101‑1 (MQ=255)
                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATttatta                                                                    <  1:3755404/101‑1 (MQ=255)
                             gATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATttatta                                                                    <  1:3755346/101‑1 (MQ=255)
                              atatCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTAc                                                                   >  2:3755962/1‑101 (MQ=255)
                                 tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                >  2:3755965/1‑101 (MQ=255)
                                 tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                <  1:3755706/101‑1 (MQ=255)
                                 tCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTa                                                                <  1:3755277/101‑1 (MQ=255)
                                  cTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTaa                                                               <  2:3755526/101‑1 (MQ=255)
                                   tttATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAAc                                                              <  2:3755248/101‑1 (MQ=255)
                                     tATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCa                                                            <  1:3755463/101‑1 (MQ=255)
                                      aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCaa                                                           >  1:3755969/1‑101 (MQ=255)
                                      aTTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCaa                                                           <  1:3755159/101‑1 (MQ=255)
                                       ttCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAAcc                                                             >  2:3755973/1‑98 (MQ=255)
                                        tCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATg                                                         <  1:3755603/101‑1 (MQ=255)
                                        tCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATg                                                         >  1:3755974/1‑101 (MQ=255)
                                         cATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGa                                                        <  2:3755381/101‑1 (MQ=255)
                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATa                                                     <  2:3755665/101‑1 (MQ=255)
                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATa                                                     >  2:3755979/1‑101 (MQ=255)
                                            aaGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATa                                                     <  2:3755493/101‑1 (MQ=255)
                                             aGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATaa                                                    <  1:3755737/101‑1 (MQ=255)
                                              gAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAAc                                                   >  1:3755980/1‑101 (MQ=255)
                                                aCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                 <  1:3755039/101‑1 (MQ=255)
                                                aCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                 >  1:3755985/1‑101 (MQ=255)
                                                aCTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACaa                                                 >  2:3755984/1‑101 (MQ=255)
                                                 cTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                <  2:3755829/101‑1 (MQ=255)
                                                 cTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAg                                                <  2:3755553/101‑1 (MQ=255)
                                                  tATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGc                                               <  1:3755535/101‑1 (MQ=255)
                                                    ttCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCtt                                             >  2:3755990/1‑101 (MQ=255)
                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaa                                          >  2:3755994/1‑101 (MQ=255)
                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaa                                          <  2:3755095/101‑1 (MQ=255)
                                                       atcaACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaa                                          >  1:3755993/1‑101 (MQ=255)
                                                          aaCAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTaaaaaa                                       >  2:3755996/1‑101 (MQ=255)
                                                            cAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAata                                     <  2:3755637/101‑1 (MQ=255)
                                                            cAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAata                                     <  1:3755784/101‑1 (MQ=255)
                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                    <  1:3755389/101‑1 (MQ=255)
                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                    <  2:3755478/101‑1 (MQ=255)
                                                             aGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAatat                                    >  1:3756000/1‑101 (MQ=255)
                                                              gCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAAtata                                   >  2:3756002/1‑101 (MQ=255)
                                                               cGTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATAc                                  <  1:3755628/101‑1 (MQ=255)
                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                 >  1:3756006/1‑101 (MQ=255)
                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                 <  1:3755572/101‑1 (MQ=255)
                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                 <  2:3755725/101‑1 (MQ=255)
                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                 <  1:3755502/101‑1 (MQ=255)
                                                                gTATAGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACt                                 >  2:3756007/1‑101 (MQ=255)
                                                                 tataGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTg                                >  1:3756011/1‑101 (MQ=255)
                                                                  ataGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGt                               >  2:3756013/1‑101 (MQ=255)
                                                                  ataGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGt                               <  1:3755884/101‑1 (MQ=255)
                                                                   taGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                              >  2:3756016/1‑101 (MQ=255)
                                                                   taGAGGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                              <  1:3755519/101‑1 (MQ=255)
                                                                    agagGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGtt                              >  2:3756019/1‑100 (MQ=255)
                                                                      agGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAtt                           <  2:3755414/101‑1 (MQ=255)
                                                                      agGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAtt                           >  2:3756022/1‑101 (MQ=255)
                                                                      agGCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTAtt                           <  1:3755484/101‑1 (MQ=255)
                                                                        gCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCt                         <  2:3755511/101‑1 (MQ=255)
                                                                        gCGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCt                         <  2:3755368/101‑1 (MQ=255)
                                                                         cGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTa                        <  1:3755452/101‑1 (MQ=255)
                                                                         cGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTa                        <  1:3755449/101‑1 (MQ=255)
                                                                         cGGTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTa                        <  2:3755690/101‑1 (MQ=255)
                                                                          ggTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                       >  1:3756025/1‑101 (MQ=255)
                                                                          ggTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaa                       <  1:3755676/101‑1 (MQ=255)
                                                                           gTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                      <  1:3755640/101‑1 (MQ=255)
                                                                           gTTATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTaaa                      >  2:3756027/1‑101 (MQ=255)
                                                                            ttATGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAg                     >  1:3756030/1‑101 (MQ=255)
                                                                              aTGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaa                   <  1:3755615/101‑1 (MQ=255)
                                                                              aTGTAAAACCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaa                   >  2:3756032/1‑101 (MQ=255)
                                                                                  aaaaCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaaa               <  1:3755562/101‑1 (MQ=255)
                                                                                  aaaaCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGaaaaaa               >  2:3756036/1‑101 (MQ=255)
                                                                                   aaaCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAt              >  1:3756040/1‑101 (MQ=255)
                                                                                   aaaCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAt              >  1:3756037/1‑101 (MQ=255)
                                                                                    aaCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAAtt             <  1:3755337/101‑1 (MQ=255)
                                                                                     aCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTa            <  1:3755653/101‑1 (MQ=255)
                                                                                     aCCACTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTa            <  2:3755794/101‑1 (MQ=255)
                                                                                        aCTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACa         <  2:3755897/101‑1 (MQ=255)
                                                                                        aCTCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACa         <  2:3755613/101‑1 (MQ=255)
                                                                                          tCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAgg       >  1:3756046/1‑101 (MQ=255)
                                                                                          tCATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAgg       >  2:3756045/1‑101 (MQ=255)
                                                                                           cATTAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGc      >  1:3756049/1‑101 (MQ=255)
                                                                                              tAGCTTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCAtt   <  2:3755565/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TTATCAATTTTCACTTAAAGTATGAACTGGATATCTTTATTCATAAGAACTATTCATCAACAGCGTATAGAGGCGGTTATGTAAAACCACTCATTAGCCTCAAAACAAAGACAAACTTCCAGATTTATTACTTAACCAATGAATAACAAGCTTAAAAAATATACTGTTATTCTAAAGAAAAAATTAACAGGCATTA  >  NC_000913/2903606‑2903801

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