Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,548,617 A→C L111L (CTA→CTC)  fimG → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,6170AC98.7% 808.4 / ‑4.3 230L111L (CTA→CTC) fimGminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (107/119);  total (110/120)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

GGGAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGA  >  NC_000913/4548521‑4548713
                                                                                                |                                                                                                
 ggAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAgg                                                                                              <  2:3801273/100‑1 (MQ=255)
 ggAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAgg                                                                                              >  1:4291347/1‑100 (MQ=255)
   aaCGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGAt                                                                                            <  1:1591644/100‑1 (MQ=255)
    aCGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATg                                                                                           <  2:5113607/100‑1 (MQ=255)
    aCGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATg                                                                                           <  2:4252004/100‑1 (MQ=255)
        cGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAg                                                                                       <  1:1137700/100‑1 (MQ=255)
        cGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAg                                                                                       <  2:2322812/100‑1 (MQ=255)
        cGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAg                                                                                       <  1:216718/100‑1 (MQ=255)
         gAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGt                                                                                      >  2:298500/1‑100 (MQ=255)
          aGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTg                                                                                     <  2:2553307/100‑1 (MQ=255)
           gggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTgg                                                                                    <  1:3635955/100‑1 (MQ=255)
           gggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTgg                                                                                    >  1:2678603/1‑100 (MQ=255)
           gggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTgg                                                                                    <  2:1538650/100‑1 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   >  1:4401359/1‑100 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   >  1:4105308/1‑100 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   <  1:557090/100‑1 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   <  1:2769524/100‑1 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   >  2:1873594/1‑100 (MQ=255)
            ggTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGc                                                                                   <  1:2935082/100‑1 (MQ=255)
             gTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCa                                                                                  >  2:5007302/1‑100 (MQ=255)
              tCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCaa                                                                                 >  1:3789564/1‑100 (MQ=255)
                aCTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAaca                                                                               <  2:1786818/100‑1 (MQ=255)
                 cTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAacac                                                                              >  2:468810/1‑100 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            <  2:3239581/100‑1 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            >  1:3516301/1‑100 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            >  1:3645173/1‑100 (MQ=255)
                   gCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACAt                                                                            >  1:3974935/1‑100 (MQ=255)
                     cAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTg                                                                          <  2:2161372/100‑1 (MQ=255)
                     cAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTg                                                                          <  1:155020/100‑1 (MQ=255)
                      aacTTCAGCGCGGCACCCGACACTACCGGATATTATAACAACCATGCGACCGCGCAAAACATCCAGCTAGAGCTCCAGGATGACAGTGGCAACACATTGa                                                                         <  2:4557953/98‑1 (MQ=255)
                      aGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGa                                                                         <  2:1060307/100‑1 (MQ=255)
                       gCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGaa                                                                        <  2:4866142/100‑1 (MQ=255)
                       gCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGaa                                                                        >  2:2887324/1‑100 (MQ=255)
                       gCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGaa                                                                        >  1:4935340/1‑100 (MQ=255)
                        cTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAAt                                                                       <  1:1019412/100‑1 (MQ=255)
                        cTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAAt                                                                       >  1:5048814/1‑100 (MQ=255)
                           cAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACt                                                                    <  2:864472/100‑1 (MQ=255)
                           cAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACt                                                                    <  2:2737181/100‑1 (MQ=255)
                           cAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACt                                                                    <  1:3239061/100‑1 (MQ=255)
                            aGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTg                                                                   >  1:5361239/1‑100 (MQ=255)
                             gCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTgg                                                                  <  1:4190478/100‑1 (MQ=255)
                                 ggCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCa                                                              <  2:1462147/100‑1 (MQ=255)
                                   cAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAAc                                                            <  2:2991338/100‑1 (MQ=255)
                                   cAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAAc                                                            <  1:1267540/100‑1 (MQ=255)
                                    aGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAAcc                                                           <  2:3529308/100‑1 (MQ=255)
                                     gCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCa                                                          >  1:2277574/1‑100 (MQ=255)
                                     gCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCa                                                          >  2:2801635/1‑100 (MQ=255)
                                         aCAGTACCGGATATTATAAAAACCAGGGGACCGTGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAAc                                                      <  1:5515437/100‑1 (MQ=255)
                                         aCAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGTGCTCCAGGATGACAGTGGCACCACATTGAATACTGGCGCAACCAAAAc                                                      <  2:5468548/100‑1 (MQ=255)
                                         aCAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAAc                                                      <  1:2828615/100‑1 (MQ=255)
                                          cAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACa                                                     >  1:5283708/1‑100 (MQ=255)
                                           aGTACCGGATATTATAAAAACGAGGGGACCGCGCAAATCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAg                                                    <  1:5189060/100‑1 (MQ=255)
                                            gTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCTCAACCAAAACAGt                                                   >  2:4371580/1‑100 (MQ=255)
                                             tACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTGAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAACACAGta                                                  >  2:1170077/1‑99 (MQ=255)
                                              aCCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTc                                                 <  2:2821704/100‑1 (MQ=255)
                                              aCCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTc                                                 >  2:3093976/1‑100 (MQ=255)
                                               ccGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCa                                                >  1:1675482/1‑100 (MQ=255)
                                               ccGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCa                                                <  1:2512422/100‑1 (MQ=255)
                                                cGGATGTTATAAAAACCAGGGGAACGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGACAACACATTGAATACTGGCGCAACCAAAACAGTTCAg                                               >  2:331564/1‑100 (MQ=255)
                                                 ggATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTTCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAgg                                              >  2:2678134/1‑100 (MQ=255)
                                                 ggATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAgg                                              <  1:457036/100‑1 (MQ=255)
                                                   atatTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTg                                            <  2:3920703/100‑1 (MQ=255)
                                                     attatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGAGGCAACACATTGAATACTGGCGCAACGAAAACAGTTGAGGTGga                                          >  1:5109884/1‑100 (MQ=255)
                                                      ttatAAAAACCANGAGACCNCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACNAAAACAGTTCAGGTGgat                                         <  1:4215812/100‑1 (MQ=255)
                                                      ttatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgat                                         <  1:4364707/100‑1 (MQ=255)
                                                      ttatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgat                                         >  2:2920012/1‑100 (MQ=255)
                                                      ttatAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgat                                         <  1:3169538/100‑1 (MQ=255)
                                                        ataAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgatga                                       >  1:1630531/1‑100 (MQ=255)
                                                        ataAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgatga                                       >  2:4352167/1‑100 (MQ=255)
                                                         taAAAACCAGGGGACCGCGCAATACATCCAGTTAGAGCTCCAGGATGACAGGGGCAACATATTGAATACTGGCGCAACCAAAACAGTGCAGGTGgatgag                                      >  2:534132/1‑99 (MQ=255)
                                                         taAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgatgat                                      >  2:2286953/1‑100 (MQ=255)
                                                         taAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGgatgat                                      <  2:5182654/100‑1 (MQ=255)
                                                          aaaaaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGAtt                                     >  2:4244766/1‑100 (MQ=255)
                                                           aaaaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTc                                    >  2:1390903/1‑100 (MQ=255)
                                                           aaaaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCCCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTc                                    <  1:914904/100‑1 (MQ=255)
                                                             aaCCAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCt                                  <  1:2906479/100‑1 (MQ=255)
                                                               ccAGGGGACCGCGCAAAACAGCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTca                                <  2:714983/100‑1 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               <  1:3138415/100‑1 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               >  1:5475125/1‑100 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               >  1:2291762/1‑100 (MQ=255)
                                                                cAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcac                               >  2:5323145/1‑100 (MQ=255)
                                                                 aGGGGACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTcaca                              <  1:3014141/100‑1 (MQ=255)
                                                                   gggACCGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAAt                            >  1:645002/1‑100 (MQ=255)
                                                                       ccGCTCAAAACATCCAATTAGAACTCGACGATGACCCTCGCAATACATTGCATACTCGCGCAACCAAAACAGTTCAGCTAGATCATCCATCACAAtaggc                        <  2:450819/100‑5 (MQ=255)
                                                                       ccGCGCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAgc                        <  2:1287825/100‑1 (MQ=255)
                                                                        cgcgCAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAgcg                       <  1:3286257/100‑1 (MQ=255)
                                                                         gcgcAAAACATCCAGTTAGAGCTCTAGTATGACAGCGGCCACCCATATAATCACGGCGCAACCACAACAGTTCAGGTATATGTCTCCTCAAAATGAgcgc                      <  2:4970224/100‑1 (MQ=255)
                                                                           gcAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCAc                    >  1:1398969/1‑100 (MQ=255)
                                                                            cAAAGCATCCAGGTAGAGCGCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACt                   >  2:4416471/1‑100 (MQ=255)
                                                                            cAAAACATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACt                   >  1:466474/1‑100 (MQ=255)
                                                                              aaaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTc                 >  1:2227431/1‑100 (MQ=255)
                                                                               aaCATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTcc                >  1:3439941/1‑100 (MQ=255)
                                                                                 gctcgaGTTAGAGCTACAGGCTGACAGTGGAAACACTTTGAATGCTGGCGCAACCGAAACAGATCAGGTGGATGATTCCTCACAATCAGCGCACATCACg              <  2:2484106/95‑1 (MQ=255)
                                                                                 cATCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCg              <  2:4068136/100‑1 (MQ=255)
                                                                                  aTCCAGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAACCAGCGCACTTCCCGt             <  2:2254250/100‑1 (MQ=255)
                                                                                      aGTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACa         <  1:3249333/100‑1 (MQ=255)
                                                                                       gTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAg        <  1:1875940/100‑1 (MQ=255)
                                                                                       gTTAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAg        <  2:4176477/100‑1 (MQ=255)
                                                                                         tAGAGCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGt      >  2:1837966/1‑100 (MQ=255)
                                                                                             gCTCCAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCaga  >  2:1227137/1‑100 (MQ=255)
                                                                                                |                                                                                                
GGGAACGTCGAGGGTCACTGCCAGCTTCAGCGGGGCAGCCGACAGTACCGGATATTATAAAAACCAGGGGACCGCGCAAAACATCCAGTTAGAGCTACAGGATGACAGTGGCAACACATTGAATACTGGCGCAACCAAAACAGTTCAGGTGGATGATTCCTCACAATCAGCGCACTTCCCGTTACAGGTCAGA  >  NC_000913/4548521‑4548713

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