Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,542,142 T→G Y36D (TAT→GAT)  fimE → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,542,1420TG99.1% 792.1 / ‑2.1 223Y36D (TAT→GAT) fimEtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (105/116);  total (106/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 223 total aligned reads displayed.

CGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACG  >  NC_000913/4542046‑4542238
                                                                                                |                                                                                                
cgtcgtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATc                                                                                               <  1:3378032/100‑1 (MQ=255)
 gtcgtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGTGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCg                                                                                              <  2:4295742/100‑1 (MQ=255)
   cgtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGc                                                                                            >  1:178938/1‑100 (MQ=255)
    gtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCGGATCGGCa                                                                                           <  2:5273236/100‑1 (MQ=255)
    gtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGcn                                                                                           <  1:1497325/100‑2 (MQ=255)
    gtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCa                                                                                           >  2:2817703/1‑100 (MQ=255)
    gtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCa                                                                                           >  1:3128169/1‑100 (MQ=255)
    gtTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCa                                                                                           >  2:2110618/1‑100 (MQ=255)
      tATCTTACCGGTAAAGAAGTTCAGGCCATTATGCACGCGGTTTATTACAGGCAAACGGGAGCCAGAGATTATTGTCTTCTCCTGTTGGCAGATCGGCATg                                                                                         >  2:2031305/1‑100 (MQ=255)
      tATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATg                                                                                         >  1:593381/1‑100 (MQ=255)
      tATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATg                                                                                         <  2:2555860/100‑1 (MQ=255)
          ttACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCAtcccat                                                                                     >  2:1994002/1‑95 (MQ=255)
           tACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATg                                                                                    <  2:3464593/100‑1 (MQ=255)
            aCCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTAGTCTGTTGGCAGATCGGCATGGGATGc                                                                                   <  2:3623194/100‑1 (MQ=255)
             ccGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCg                                                                                  >  1:4384823/1‑100 (MQ=255)
             ccGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCg                                                                                  >  2:1899446/1‑100 (MQ=255)
             ccGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTATTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCg                                                                                  <  2:1244105/100‑1 (MQ=255)
              cGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGt                                                                                 <  1:4021385/100‑1 (MQ=255)
               ggTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTa                                                                                <  1:222874/100‑1 (MQ=255)
               ggTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTa                                                                                >  2:2458353/1‑100 (MQ=255)
                 tAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTAtt                                                                              >  2:2176289/1‑100 (MQ=255)
                  aaaGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTa                                                                             >  1:2058874/1‑100 (MQ=255)
                   aagaagTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAg                                                                            <  1:1329579/100‑1 (MQ=255)
                   aagaagTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAg                                                                            >  1:3933753/1‑100 (MQ=255)
                    agaagTTCAGGCCATGATGCAGGCGGTATGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGt                                                                           >  1:4643294/1‑100 (MQ=255)
                     gaagTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTg                                                                          >  1:4261557/1‑100 (MQ=255)
                     gaagTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTg                                                                          >  2:2962936/1‑100 (MQ=255)
                       agTACACGCCATGATGCAGGCGGTTTGTTACTGGGCAACGGGAGCCAGAGGTTAGTGTCCTATTCTGTGGGCAGATCGGTATGGGATGCGTATTAGtctc                                                                        >  2:701925/1‑97 (MQ=255)
                        gTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAAc                                                                       >  2:635523/1‑100 (MQ=255)
                        gTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAAc                                                                       <  2:10999/100‑1 (MQ=255)
                          tCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTATGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTg                                                                     >  1:1381245/1‑100 (MQ=255)
                           cAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCGGTTGGCAGATCGGCATGGGCTGCGTATTAGTGAACTGc                                                                    >  2:5307118/1‑100 (MQ=255)
                            nggCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCt                                                                   >  1:4443672/2‑100 (MQ=255)
                            aGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCt                                                                   <  2:5261330/100‑1 (MQ=255)
                             ggCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCtt                                                                  <  2:2863263/100‑1 (MQ=255)
                             ggCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCtt                                                                  >  2:4431133/1‑100 (MQ=255)
                             ggCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCtt                                                                  <  2:5442591/100‑1 (MQ=255)
                              gCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTg                                                                 >  1:4735950/1‑100 (MQ=255)
                               ccATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGa                                                                <  1:3833447/100‑1 (MQ=255)
                                cATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGAt                                                               >  2:1685363/1‑100 (MQ=255)
                                cATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGAt                                                               >  2:506629/1‑100 (MQ=255)
                                 atgatgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATc                                                              <  1:5346412/100‑1 (MQ=255)
                                 atgatgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATc                                                              <  2:1928720/100‑1 (MQ=255)
                                  tgatgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCt                                                             <  1:1970443/100‑1 (MQ=255)
                                   gatgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTg                                                            <  2:4915224/100‑1 (MQ=255)
                                     tgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCa                                                          >  1:1739702/1‑100 (MQ=255)
                                     tgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCa                                                          <  1:423354/100‑1 (MQ=255)
                                     tgCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCa                                                          >  2:3820691/1‑100 (MQ=255)
                                      gCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCAt                                                         <  1:3478493/100‑1 (MQ=255)
                                       cAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCAtt                                                        >  1:5161291/1‑100 (MQ=255)
                                          gCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATc                                                     >  1:5189095/1‑100 (MQ=255)
                                              tttGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGga                                                 >  2:4715377/1‑100 (MQ=255)
                                              tttGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGga                                                 >  1:1185674/1‑100 (MQ=255)
                                              tttGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGga                                                 >  1:3276881/1‑100 (MQ=255)
                                                tGTTACGGTGCAACGGGGGCCAGAGATTATTGTCTAATTCTGTGGGCAGATCCGCATCGGATGCGTATTAGTGCACTGCTTGATCTGCATTATCAGgaac                                               <  1:1713512/100‑3 (MQ=255)
                                                 gTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                              >  2:3182994/1‑100 (MQ=255)
                                                 gTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                              >  1:1532111/1‑100 (MQ=255)
                                                 gTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgacct                                              >  1:1069775/1‑100 (MQ=255)
                                                   tACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttg                                            >  1:2414133/1‑100 (MQ=255)
                                                     cGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgac                                          <  2:3897832/100‑1 (MQ=255)
                                                     cGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGATGGCAGATCGGCATGGGATGCGTAATAGTGAACTACTTGATCTGCATTATCAGggccttgac                                          >  2:4079013/1‑100 (MQ=255)
                                                       gggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgacct                                        <  2:5501083/100‑1 (MQ=255)
                                                       gggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgacct                                        >  1:2565040/1‑100 (MQ=255)
                                                        ggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgacctt                                       >  2:1774946/1‑100 (MQ=255)
                                                        ggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgacctt                                       <  2:329507/100‑1 (MQ=255)
                                                        ggCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGgaccttgacctt                                       <  1:3239254/100‑1 (MQ=255)
                                                         gCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTa                                      <  1:2171458/100‑1 (MQ=255)
                                                          cAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTaa                                     >  2:4623280/1‑100 (MQ=255)
                                                           aaCGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAAt                                    <  2:1390684/100‑1 (MQ=255)
                                                           aaCGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAAt                                    <  1:5109508/100‑1 (MQ=255)
                                                            aCGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATg                                   >  1:997058/1‑100 (MQ=255)
                                                            aCGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATg                                   <  1:2338127/100‑1 (MQ=255)
                                                             cGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGa                                  <  2:2726867/100‑1 (MQ=255)
                                                               ggAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAg                                <  1:1810525/100‑1 (MQ=255)
                                                               ggAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAg                                <  2:3643833/100‑1 (MQ=255)
                                                               ggAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAg                                >  2:5094750/1‑100 (MQ=255)
                                                                gAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGTCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAgg                               <  2:4990334/100‑1 (MQ=255)
                                                                gAGCCAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAgg                               <  1:2932974/100‑1 (MQ=255)
                                                                   ccAGAGATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAg                            >  1:5442233/1‑100 (MQ=255)
                                                                     agagATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaa                          >  2:5161168/1‑100 (MQ=255)
                                                                     agagATTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaa                          >  1:5071243/1‑100 (MQ=255)
                                                                      gagaTTATTGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaat                         <  2:845578/100‑1 (MQ=255)
                                                                         attattGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaa                      <  2:1719609/100‑1 (MQ=255)
                                                                          ttattGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaat                     <  2:4212189/100‑1 (MQ=255)
                                                                             ttGTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATAtt                  >  1:4572710/1‑100 (MQ=255)
                                                                              tGTCTTATTCTGTTGGCAGATCGGCATGGGATTCGTATTAGTGAACTGCTTGATCTGGATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTc                 <  2:4651961/100‑1 (MQ=255)
                                                                               gTCTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCg                >  1:1694566/1‑100 (MQ=255)
                                                                                 cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGcc              >  1:2475581/1‑100 (MQ=255)
                                                                                 cTTATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGcc              >  1:630998/1‑100 (MQ=255)
                                                                                   tATTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGa            >  1:1261240/1‑100 (MQ=255)
                                                                                    aTTCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGAc           <  2:1597592/100‑1 (MQ=255)
                                                                                      tCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTg         <  1:1574238/100‑1 (MQ=255)
                                                                                      tCTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTg         >  1:4165629/1‑100 (MQ=255)
                                                                                       cTGTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTga        <  2:3741173/100‑1 (MQ=255)
                                                                                         gTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaag      >  2:5046644/1‑100 (MQ=255)
                                                                                         gTTGGCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaag      <  1:1546981/100‑1 (MQ=255)
                                                                                           tGGCAGATCGGCATGGGATGCGTATTAGTGACCTGCTTGATATGCATTATCAGGACCTTGACCTTAATGAGGGAAGAATAAATATTCGCCGACTgaagaa    >  2:5229750/1‑100 (MQ=255)
                                                                                            ggCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAAc   >  1:3655339/1‑100 (MQ=255)
                                                                                            ggCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAAc   >  2:230132/1‑100 (MQ=255)
                                                                                            ggCAGATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAAc   >  2:1498625/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACG  >  NC_000913/4542046‑4542238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑