Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,032 C→A 16.5% intergenic (+22/‑58) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0320CA16.5% 269.8 / 36.9 102intergenic (+22/‑58)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (45/40);  new base A (6/11);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 104 total aligned reads displayed.

GTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTG  >  NC_000913/3287933‑3288128
                                                                                                   |                                                                                                
gTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGAcc                                                                                                  <  2:2564382/100‑1 (MQ=255)
  ccgatctgTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTg                                                                                                <  1:1238103/93‑1 (MQ=255)
  aTCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTCTCTGACCTg                                                                                                >  2:2503651/1‑100 (MQ=255)
   tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCt                                                                                                 <  1:2544651/98‑1 (MQ=255)
   tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGt                                                                                               <  2:2639544/100‑1 (MQ=255)
      aCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAg                                                                                            <  2:1447498/100‑1 (MQ=255)
      aCCGTAGTTCCATCAGAAATCACCCAAGGGGGAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTc                                                                                              >  1:1011549/1‑98 (MQ=255)
         gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTgagatcgg                                                                                         >  2:1238108/1‑93 (MQ=255)
         gTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATa                                                                                         <  1:1710601/100‑1 (MQ=255)
            gTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACg                                                                                      <  2:2810672/100‑1 (MQ=255)
             ttagagcataaaTCACCCCAGGGGAATCAGGCGCTATCGGTAACTTAACGCTGGAATACCGTTAATACTTTAGGGCGTTATCTGACAGGTCCGATAACgt                                                                                     <  1:2311741/91‑2 (MQ=255)
             ttCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                     <  2:2695895/100‑1 (MQ=255)
              tCCATCAGAATTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATATGTCAGGGCGTTATCTGACAGGTCAGATAACGcc                                                                                    <  1:1911554/100‑1 (MQ=38)
              tCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGcc                                                                                    <  1:2896679/100‑1 (MQ=255)
              tCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGcc                                                                                    <  1:373016/100‑1 (MQ=255)
               ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGccc                                                                                   <  2:154354/100‑1 (MQ=255)
               ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                     <  1:423553/98‑1 (MQ=255)
                cATCAGAAATCACTCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                  <  1:1445251/100‑1 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                  >  1:2209753/1‑100 (MQ=255)
                cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                  >  2:2594650/1‑100 (MQ=255)
                     gAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                             >  2:1890329/1‑100 (MQ=255)
                     gAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                             >  2:2234429/1‑100 (MQ=255)
                      aaaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcct                                                                            >  2:183040/1‑100 (MQ=255)
                      aaaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                             <  1:671564/99‑1 (MQ=255)
                       aaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                           >  1:341434/1‑100 (MQ=255)
                       aaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGACGTTATCTGACCTGTCAGATAACGCCCTTttcctt                                                                           >  1:1480212/1‑100 (MQ=255)
                        aTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGCGCGTTATCTGACCTGTCAGATAACGCCCTTttccttc                                                                          >  2:827112/1‑100 (MQ=255)
                          cACCCCAGGGGAATCAGACGCTATCGTTAACGTAACGCTGGCTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttcc                                                                         <  1:1984292/99‑1 (MQ=255)
                           aCCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                       <  1:2528736/100‑1 (MQ=255)
                            cccaAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctct                                                                      <  1:2289068/100‑1 (MQ=255)
                             cccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCtt                                                                     >  2:2230396/1‑100 (MQ=255)
                                 ggagAAGGAGACTCTGTGGTGAGCGGAAGGATGGATTACCGGGAAGACGTTAGGGCGGGATCTGACCTGTCAGATAACGccc                                                                                   <  1:1684569/79‑1 (MQ=255)
                                  gggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCg                                                                <  1:2811410/100‑1 (MQ=255)
                                    gAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                              <  1:2267993/100‑1 (MQ=255)
                                    gAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                              >  1:1028511/1‑100 (MQ=255)
                                     aaGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTg                                                             <  2:59525/100‑1 (MQ=255)
                                     aaGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTg                                                             <  1:1271138/100‑1 (MQ=255)
                                     aaGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTg                                                             <  2:1377569/100‑1 (MQ=255)
                                          gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCa                                                        >  2:123927/1‑100 (MQ=255)
                                          gACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCa                                                        >  2:1221310/1‑100 (MQ=255)
                                           aCGCTATCGTTAACGTAACGCTGGACTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATACCGTCCTTTTCCTTCCTCTTTCTCGt                                                               <  1:1772430/92‑1 (MQ=255)
                                            cGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAgg                                                      >  1:1061881/1‑100 (MQ=255)
                                             gCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                     <  1:1651411/100‑1 (MQ=255)
                                               tATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTg                                                   <  1:2514521/100‑1 (MQ=255)
                                                 tCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaa                                                 >  2:648909/1‑100 (MQ=255)
                                                  cGTTAACGTACCGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGTCCTGTCAGATCACGCCCTTTTCCTGCCTCTTTCTCGTTGTATCAGGTTGaaa                                                <  1:240239/100‑1 (MQ=255)
                                                  cGTTAACGTAACGCTTGAGTACCTTTAATACGTTAGTGCGTTGTCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTg                                                   <  1:1002533/97‑1 (MQ=255)
                                                  cGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaa                                                <  1:1730173/100‑1 (MQ=255)
                                                      aacgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATa                                            >  1:2145528/1‑100 (MQ=255)
                                                        cgtaacgCTGTATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAAc                                          <  1:1264156/100‑1 (MQ=255)
                                                         gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACt                                         <  2:762993/100‑1 (MQ=255)
                                                         gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACt                                         <  1:1672023/100‑1 (MQ=255)
                                                          taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                        >  2:1025272/1‑100 (MQ=255)
                                                          taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                        <  2:2458204/100‑1 (MQ=255)
                                                          taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                        >  2:2310808/1‑100 (MQ=255)
                                                            acgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                        <  1:2675979/98‑1 (MQ=255)
                                                             cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGt                                     >  2:521121/1‑100 (MQ=255)
                                                             cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGt                                     >  2:517889/1‑100 (MQ=255)
                                                             cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGt                                     >  1:2698224/1‑100 (MQ=255)
                                                             cgCTGGATTACCGTTAATACGTTACGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTAGTTGTATCAGGTTGAAAAATAACTATGt                                     >  1:328774/1‑100 (MQ=255)
                                                              gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTc                                    >  1:984520/1‑100 (MQ=255)
                                                               cTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCa                                   <  2:826214/100‑1 (MQ=255)
                                                               cTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCa                                   <  2:2063487/100‑1 (MQ=255)
                                                                  gATTACCGTTAATACGTTAGGGCGTTATCTGACCTGCCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaaaa                                >  2:2052134/1‑100 (MQ=255)
                                                                   aTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                               <  1:2461241/100‑1 (MQ=255)
                                                                   aTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAAc                               >  1:629976/1‑100 (MQ=255)
                                                                   aTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTATATCAGGTTGAAAAATAACTATGTCAAAAc                               >  1:2716504/1‑100 (MQ=255)
                                                                       ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAAc                           >  2:63996/1‑100 (MQ=255)
                                                                       ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACGATGTCAAAACGAAc                           >  1:1983754/1‑100 (MQ=255)
                                                                            aaTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcg                      <  2:1547455/100‑1 (MQ=255)
                                                                             aTACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgc                     >  2:2337533/1‑100 (MQ=255)
                                                                              tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTggcg                    >  1:605866/1‑100 (MQ=255)
                                                                                cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                  >  1:293511/1‑100 (MQ=255)
                                                                                 gTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTg                 <  1:2736248/100‑1 (MQ=255)
                                                                                 gTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTg                 <  1:2014581/100‑1 (MQ=255)
                                                                                  ttAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                <  1:742101/100‑1 (MQ=255)
                                                                                  ttAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTCGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                >  1:1433132/1‑100 (MQ=255)
                                                                                   tAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGat               <  1:1105813/100‑1 (MQ=255)
                                                                                   tAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGat               <  1:2052123/100‑1 (MQ=255)
                                                                                    aGGGCGTTATCTGACCTGTCAGATAACTCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGata              <  2:2046894/100‑1 (MQ=255)
                                                                                    aGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGata              >  1:1116189/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAGCTATGTCAAAACGAACATTCGCGGTGATAtt            >  2:667842/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt            <  2:2301061/100‑1 (MQ=255)
                                                                                      ggCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt            >  2:2512952/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt            >  2:1127794/1‑100 (MQ=255)
                                                                                      ggCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGCGATAtt            >  1:1947068/1‑100 (MQ=38)
                                                                                       gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa           >  2:1013077/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa           >  2:637907/1‑100 (MQ=255)
                                                                                       gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCCTTGTATCAGGTTGAAAAATACCTATGTCAAAACGAACATTCGCGGTGATATTa           >  1:2707526/1‑100 (MQ=38)
                                                                                        cGTTATCTGACCTGTCAGATGACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTaa          >  2:1320242/1‑100 (MQ=255)
                                                                                        cGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTaa          <  2:342143/100‑1 (MQ=255)
                                                                                        cGTTATCTGACATGTCAGATAACGCCCTTTTCCTTTCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTaa          <  2:1181963/100‑1 (MQ=255)
                                                                                         gTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc         >  1:2877312/1‑100 (MQ=255)
                                                                                         gTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc         >  2:182015/1‑100 (MQ=255)
                                                                                           tATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCt       >  2:1225985/1‑100 (MQ=255)
                                                                                               tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgtt   <  2:2331638/100‑1 (MQ=255)
                                                                                               tGACCTGTCAGATAACCCCCTATTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTCGGTCAAAACGAACATTAGCGGTGATATTAACTttg     >  1:2641145/1‑98 (MQ=255)
                                                                                                gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg  <  1:2502606/100‑1 (MQ=255)
                                                                                                gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg  >  2:1356264/1‑100 (MQ=255)
                                                                                                  ccTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgt    >  1:1737814/1‑96 (MQ=255)
                                                                                                   |                                                                                                
GTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTG  >  NC_000913/3287933‑3288128

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