Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,520,438 A→C L89R (CTC→CGC)  hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,4380AC98.8% 854.7 / ‑0.7 244L89R (CTC→CGC) hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (111/130);  total (113/131)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 246 total aligned reads displayed.

CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCG  >  NC_000913/3520342‑3520534
                                                                                                |                                                                                                
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGaa                                                                                               >  1:2752214/1‑100 (MQ=255)
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGaa                                                                                               >  2:2670527/1‑100 (MQ=255)
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGaa                                                                                               >  1:263304/1‑100 (MQ=255)
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGaa                                                                                               <  2:4911374/100‑1 (MQ=255)
 tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAAt                                                                                              >  1:3223100/1‑100 (MQ=255)
  gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATg                                                                                             <  2:4959756/100‑1 (MQ=255)
  gCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGACAGGCAATGGAATGCACTGAGCGAATg                                                                                             <  1:1116145/100‑1 (MQ=255)
    gATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGtt                                                                                           <  2:4452526/100‑1 (MQ=255)
     aTTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCGCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTg                                                                                          <  2:3321537/100‑1 (MQ=255)
        ttCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCt                                                                                       <  1:4610890/100‑1 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTg                                                                                  <  1:236973/100‑1 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTg                                                                                  <  1:3929097/100‑1 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTg                                                                                  >  1:3133378/1‑100 (MQ=255)
              cGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGc                                                                                 >  1:5901/1‑100 (MQ=255)
               ggCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGcc                                                                                >  2:5191735/1‑100 (MQ=255)
                gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCcg                                                                               <  1:711701/100‑1 (MQ=255)
                 cAGATTTGCCTGCGCCCGCACCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGGTGCCTTCCTGCcgc                                                                              >  2:5031788/1‑100 (MQ=255)
                  gatcttGCCGGCGCCCGCGCCTGCTCCGCCTCCGGGCGGGGGATATTGTTATTCTGCCAGGCAATGGATTGCACTGAGCGAATGTTGCCTTCCTGCcgcg                                                                             <  2:4625378/96‑1 (MQ=255)
                     tttGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTa                                                                          <  2:725430/100‑1 (MQ=255)
                      ttGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTAt                                                                         <  1:4802970/100‑1 (MQ=255)
                       tGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCcgcgagtcg                                                                        >  1:4684523/1‑95 (MQ=255)
                       tGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATc                                                                        >  2:4423580/1‑100 (MQ=255)
                       tGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATc                                                                        <  2:3445862/100‑1 (MQ=255)
                        gCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGGCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCa                                                                       >  1:4253957/1‑100 (MQ=255)
                        gCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCa                                                                       <  1:2691887/100‑1 (MQ=255)
                         ccTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATctg                                                                      >  2:559782/1‑98 (MQ=255)
                         ccTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAg                                                                      <  1:5171966/100‑1 (MQ=255)
                         ccTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAg                                                                      <  1:5032209/100‑1 (MQ=255)
                          cTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGt                                                                     >  1:3867680/1‑100 (MQ=255)
                                ccGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGgcg                                                               >  2:2215277/1‑100 (MQ=255)
                                 cgcgCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGgcgc                                                              >  2:1219854/1‑100 (MQ=255)
                                    gcCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGCTATCAGTCCGGCGCttt                                                           >  2:2915081/1‑100 (MQ=255)
                                     ccTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCtttt                                                          >  2:306507/1‑100 (MQ=255)
                                       tGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCa                                                        >  2:167102/1‑100 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCAc                                                       <  1:3707823/100‑1 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCAc                                                       <  1:3020587/100‑1 (MQ=255)
                                        gCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCAc                                                       <  1:5028428/100‑1 (MQ=255)
                                        gCTCCGCCGCCTGGCGGGCGATCTTGTTATTCTGCCAGGCCATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCAc                                                       <  1:1147801/100‑1 (MQ=255)
                                         cTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACt                                                      >  2:4690364/1‑100 (MQ=255)
                                          tCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTa                                                     >  2:4818436/1‑100 (MQ=255)
                                           ccGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTAc                                                    <  2:241725/100‑1 (MQ=255)
                                            cGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACa                                                   >  1:2998500/1‑100 (MQ=255)
                                            cGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACa                                                   >  1:822175/1‑100 (MQ=255)
                                             gCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAg                                                  >  1:2079561/1‑100 (MQ=255)
                                              cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGt                                                 <  1:4678265/100‑1 (MQ=255)
                                              cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGt                                                 >  1:4473980/1‑100 (MQ=255)
                                               ctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGtt                                                <  2:127852/100‑1 (MQ=255)
                                                 cctGGCGGGCGATCTTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTCTCAGTCCGGCGCTTTTCACTACAGTCTg                                              >  2:516217/1‑100 (MQ=255)
                                                  ctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGt                                             <  1:2667290/100‑1 (MQ=255)
                                                  ctGGCGGGCGATATTGTTATTCTGCCAGGCAATCGAATGCACTGACCTAACGTTGCCTTCCTGCCGCGTTATAAGTCCGGCGCTTTTCATTACCCTTTGt                                             <  2:1122655/100‑1 (MQ=255)
                                                    ggcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAg                                           <  1:557892/100‑1 (MQ=255)
                                                     gcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGt                                          <  2:430675/100‑1 (MQ=255)
                                                     gcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGt                                          <  2:2812777/100‑1 (MQ=255)
                                                       gggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                        <  2:3687582/100‑1 (MQ=255)
                                                          cgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTg                                     <  1:3853659/100‑1 (MQ=255)
                                                          cgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTACTGCCGCGTTAACAGTCCGGCGCTTTTCAAGACAGTTGGTAGTGCCTg                                     >  2:2585813/1‑100 (MQ=255)
                                                           gATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGc                                    <  1:1859857/100‑1 (MQ=255)
                                                            atatTGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCt                                   <  2:811506/100‑1 (MQ=255)
                                                               ttGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTcc                                <  1:741226/100‑1 (MQ=255)
                                                               ttGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTcc                                <  1:2473350/100‑1 (MQ=255)
                                                               ttGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTcc                                >  2:1953286/1‑100 (MQ=255)
                                                                tGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                               >  2:2144298/1‑100 (MQ=255)
                                                                tGTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                               >  1:1195188/1‑100 (MQ=255)
                                                                 gTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                              <  1:5531264/100‑1 (MQ=255)
                                                                 gTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                              >  1:986130/1‑100 (MQ=255)
                                                                 gTTATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTCGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGCGCCTGCTTCCAg                              <  2:5284376/100‑1 (MQ=255)
                                                                 cttATTTTGCCACGCAATGGAAAGCACTGAGCGAATGTTCCCTTCCTGCCGCGTTATCACTCCGGCGCTTTTCACTACAGTTTGTAGTGCCCGCTTCCAg                              >  2:3470327/2‑100 (MQ=255)
                                                                  ttATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAgg                             >  1:4250200/1‑100 (MQ=255)
                                                                   tATTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAggg                            >  2:3163773/1‑100 (MQ=255)
                                                                    aTTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGa                           <  1:65458/100‑1 (MQ=255)
                                                                    aTTCTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGa                           <  1:2524409/100‑1 (MQ=255)
                                                                       cTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACa                        <  2:1444739/100‑1 (MQ=255)
                                                                       cTGCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACa                        <  2:1695304/100‑1 (MQ=255)
                                                                         gCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                      >  2:2400036/1‑100 (MQ=255)
                                                                         gCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                      >  1:5426500/1‑100 (MQ=255)
                                                                         gCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATc                      >  2:217942/1‑100 (MQ=255)
                                                                         gCCAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTGCCAGGGAACATc                      >  2:3660397/1‑100 (MQ=255)
                                                                          ccAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCt                     >  1:2828074/1‑100 (MQ=255)
                                                                           cAGGCAATGGACTGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTCCGTAGTGCCTGCTTTCAGGGAACATCTg                    <  1:4903837/100‑1 (MQ=255)
                                                                           cAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                    >  1:1459524/1‑100 (MQ=255)
                                                                           cAGGCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                    >  1:3702713/1‑100 (MQ=255)
                                                                             ggCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                  >  2:3938832/1‑100 (MQ=255)
                                                                             ggCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                  >  2:1352859/1‑100 (MQ=255)
                                                                             ggCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                  >  2:4727254/1‑100 (MQ=255)
                                                                             ggCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                  <  1:2915630/100‑1 (MQ=255)
                                                                              gCAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTa                 <  2:374530/100‑1 (MQ=255)
                                                                               cAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                <  1:4049027/100‑1 (MQ=255)
                                                                               cAATGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAg                <  1:3366917/100‑1 (MQ=255)
                                                                                aaTGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCTTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa               <  2:2967692/100‑1 (MQ=255)
                                                                                aaTGGAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa               >  2:4141430/1‑100 (MQ=255)
                                                                                   ggAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt            <  2:4332850/100‑1 (MQ=255)
                                                                                   ggAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt            <  1:432622/100‑1 (MQ=255)
                                                                                   ggAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt            >  1:2250580/1‑100 (MQ=255)
                                                                                   ggAATGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt            >  2:4021077/1‑100 (MQ=255)
                                                                                    gAATGCACTGAGCGAATGTTGCCATCCTGCCGCGTTATCAGTCCGGAGATTTTCACTACATTTTGTAGTGACTGCTTCGAGGGAACATCTGTGAGATGTa           >  2:927085/1‑100 (MQ=255)
                                                                                      aTGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAAc         >  1:1691377/1‑100 (MQ=255)
                                                                                      aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAAc         <  1:5263406/100‑1 (MQ=255)
                                                                                       tGCACTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACg        <  2:5324652/100‑1 (MQ=255)
                                                                                           cTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac    >  2:420977/1‑100 (MQ=255)
                                                                                           cTGAGCGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac    >  2:2407927/1‑100 (MQ=255)
                                                                                                |                                                                                                
CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCG  >  NC_000913/3520342‑3520534

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: ‑