Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,542,123 G→T 31.6% G573C (GGC→TGC)  feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,1230GT31.6% 215.2 / 144.8 166G573C (GGC→TGC) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (52/54);  new base T (22/29);  total (78/88)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

Only 100 of 172 total aligned reads displayed.

GATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAG  >  NC_000913/3542024‑3542222
                                                                                                   |                                                                                                   
gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTg                                                                                                     >  2:3895886/1‑100 (MQ=255)
 aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGTGTGCATGTAGATATCTGGCAGGCAACGGTTgg                                                                                                    <  1:621713/100‑1 (MQ=255)
 aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                                                    <  1:4170435/100‑1 (MQ=255)
  tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                                                   <  2:201092/100‑1 (MQ=255)
     cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTg                                                                                                <  1:2922454/100‑1 (MQ=255)
      aTCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGt                                                                                               >  1:2370024/1‑100 (MQ=255)
       tCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGtt                                                                                              >  1:1036341/1‑100 (MQ=255)
        cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGttt                                                                                             >  2:796456/1‑100 (MQ=255)
        cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                                                             <  1:2406348/100‑1 (MQ=255)
         aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTa                                                                                            <  1:3608194/100‑1 (MQ=255)
         aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTa                                                                                            >  2:3670773/1‑100 (MQ=255)
         aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                                                            <  2:1827678/100‑1 (MQ=255)
         aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                                                            >  2:3178149/1‑100 (MQ=255)
          aCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTAc                                                                                           <  1:4033524/100‑1 (MQ=255)
           cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                                                          <  2:1030878/100‑1 (MQ=255)
            gACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAg                                                                                         >  2:1661087/1‑100 (MQ=255)
             aCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAgg                                                                                        >  1:3644403/1‑100 (MQ=255)
             aCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGCCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                                                                                        <  2:402534/100‑1 (MQ=255)
                 ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCGAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                                                    <  2:454871/100‑1 (MQ=255)
                   cgcTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCAt                                                                                  <  1:557171/100‑1 (MQ=255)
                   cgcTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCAt                                                                                  >  1:3195931/1‑100 (MQ=255)
                    gcTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATg                                                                                 <  2:801944/100‑1 (MQ=255)
                     cTGGCGTCCGTCAGCCGGGTGATCACCCCGGTGTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                                                                <  1:1753767/100‑1 (MQ=255)
                      tGGCGTCCGTCAGCCGGGTGATGACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                                                                >  1:1264930/1‑99 (MQ=255)
                        gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGa                                                                             >  1:3423144/1‑100 (MQ=255)
                         cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGaa                                                                            >  2:3033312/1‑100 (MQ=255)
                          gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGaaa                                                                           <  2:2971591/100‑1 (MQ=255)
                           tccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTAGCCTGTTTACAGGTGCCATGGCGAaag                                                                          <  2:2897701/100‑1 (MQ=255)
                            ccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaaga                                                                         <  1:79654/100‑1 (MQ=255)
                             cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGCaagaa                                                                        >  1:1028773/1‑100 (MQ=255)
                             cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                                                        >  1:3944112/1‑100 (MQ=255)
                              gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaag                                                                       >  2:2216804/1‑100 (MQ=255)
                               tcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGt                                                                      <  1:870878/100‑1 (MQ=255)
                                 aGCCGGGTGATCCCCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                                                    <  2:3578722/100‑1 (MQ=255)
                                 aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                                                    <  2:3629247/100‑1 (MQ=255)
                                 aGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                                                    <  1:1381553/100‑1 (MQ=255)
                                  gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTAGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                                                   <  2:1179064/100‑1 (MQ=255)
                                    cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa                                                                  >  1:439009/1‑99 (MQ=255)
                                        tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggn                                                             <  2:1611924/100‑2 (MQ=255)
                                        tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                                                             <  2:89249/100‑1 (MQ=255)
                                        tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                                                             <  2:3897963/100‑1 (MQ=255)
                                        tGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGATCACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                                                             >  1:747034/1‑100 (MQ=255)
                                         gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                                                            >  2:4214047/1‑100 (MQ=255)
                                         gATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                                                            >  1:571180/1‑100 (MQ=255)
                                          aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTa                                                           >  1:700743/1‑100 (MQ=255)
                                          aTCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGGCTGTCTACAGGTGCCATGGCGAAAGAAGTGGCAGTGGGt                                                            >  1:4142609/1‑99 (MQ=255)
                                           tCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAc                                                          >  2:1313304/1‑100 (MQ=255)
                                           tCACCCCAGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAc                                                          >  2:2249706/1‑100 (MQ=255)
                                              cctcGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACgat                                                       >  1:1146540/4‑98 (MQ=255)
                                              ccccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                                                       >  1:213405/1‑100 (MQ=255)
                                                ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                     <  2:3944537/100‑1 (MQ=255)
                                                ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                     <  1:408110/100‑1 (MQ=255)
                                                ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                     <  1:3152480/100‑1 (MQ=255)
                                                ccGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTAGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                     <  2:3758620/100‑1 (MQ=255)
                                                 cGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                                                     <  1:1047532/99‑1 (MQ=255)
                                                  ggTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                                                   <  1:423335/100‑1 (MQ=255)
                                                   gTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAaca                                                  <  1:1054307/100‑1 (MQ=255)
                                                    tCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAacac                                                 <  1:519708/100‑1 (MQ=255)
                                                      ttCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc                                               >  1:3583861/1‑100 (MQ=255)
                                                        cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                             >  2:2540556/1‑100 (MQ=255)
                                                        cAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc                                             <  1:3339020/100‑1 (MQ=255)
                                                         aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                            <  1:273351/100‑1 (MQ=255)
                                                         aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                            >  1:4235151/1‑100 (MQ=255)
                                                         aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                            <  1:1774559/100‑1 (MQ=255)
                                                         aaGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctct                                            >  2:2855557/1‑100 (MQ=255)
                                                          aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa                                           >  2:2333527/1‑100 (MQ=255)
                                                          aGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa                                           <  2:4121655/100‑1 (MQ=255)
                                                           gCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac                                          <  1:3895646/100‑1 (MQ=255)
                                                             cAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTAGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac                                        >  2:2613469/1‑100 (MQ=255)
                                                              aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc                                       <  2:347325/100‑1 (MQ=255)
                                                              aaTTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACAcc                                       <  1:3895868/100‑1 (MQ=255)
                                                                  ggCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                   <  2:3051851/100‑1 (MQ=255)
                                                                  ggCGTGCATGAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg                                   >  2:3163198/1‑100 (MQ=255)
                                                                     gTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGaaa                                <  2:3051598/100‑1 (MQ=255)
                                                                       gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                              <  1:1423485/100‑1 (MQ=255)
                                                                       gCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAat                              >  1:1672447/1‑100 (MQ=255)
                                                                           gAAGATAACTGGCAGGCAACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTc                          <  1:719867/100‑1 (MQ=255)
                                                                            aaGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                         <  2:664519/100‑1 (MQ=255)
                                                                             aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                        <  1:1980757/100‑1 (MQ=255)
                                                                             aGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTCGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAg                        >  1:183854/1‑100 (MQ=255)
                                                                             aGATAACTGGCAGGCAACGGTTGGACTGTTTACAGGAGCCATGGAGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTAAACCGCAGAAAATATTCAg                        <  1:156319/100‑1 (MQ=255)
                                                                             aGATAACTGGCAGCCAACGGTTTGCCTGTTTACAGGTGCCATGGCGCAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCa                         >  1:2646766/1‑99 (MQ=255)
                                                                                      gcaggcaACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAagag               <  1:2789799/100‑1 (MQ=255)
                                                                                       caggcaACGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt              >  2:3829599/1‑100 (MQ=255)
                                                                                       caggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAGCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt              >  1:3252671/1‑100 (MQ=255)
                                                                                       caggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt              >  2:3243935/1‑100 (MQ=255)
                                                                                       atctcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGt              <  2:3333783/96‑1 (MQ=255)
                                                                                        aggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtt             >  1:3053633/1‑100 (MQ=255)
                                                                                         ggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc            >  1:1011600/1‑100 (MQ=255)
                                                                                         ggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc            <  2:3077987/100‑1 (MQ=255)
                                                                                         ggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc            <  1:2506181/100‑1 (MQ=255)
                                                                                         ggcaACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTc            >  2:2576532/1‑100 (MQ=255)
                                                                                           caACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGtagat          >  1:3333771/1‑96 (MQ=255)
                                                                                           caACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa          <  2:1599614/100‑1 (MQ=255)
                                                                                           caACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCaa          >  2:1212387/1‑100 (MQ=255)
                                                                                              cGGTTTGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc       <  2:2633255/100‑1 (MQ=255)
                                                                                              cGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATcc       <  1:1675208/100‑1 (MQ=255)
                                                                                                 ttGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGc    <  2:1321725/100‑1 (MQ=255)
                                                                                                  tGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCa   >  2:3368039/1‑100 (MQ=255)
                                                                                                   ggCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAg  >  1:3694890/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAGAAAATATTCAGGACGAAGAGTTCAATCCGGCAG  >  NC_000913/3542024‑3542222

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: