Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,549,932 T→A 15.1% intergenic (+222/+21) fimH → / ← gntP minor component of type 1 fimbriae/fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,9320TA15.1% 499.3 / 61.1 184intergenic (+222/+21)fimH/gntPminor component of type 1 fimbriae/fructuronate transporter
Reads supporting (aligned to +/- strand):  ref base T (74/81);  new base A (10/18);  total (85/99)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

CTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTC  >  NC_000913/4549833‑4550031
                                                                                                   |                                                                                                   
cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                     <  2:3112550/100‑1 (MQ=255)
cTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAt                                                                                                     <  1:553201/100‑1 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                                <  2:443767/100‑1 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                                >  1:2859405/1‑100 (MQ=255)
     cTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCa                                                                                                <  1:2827696/100‑1 (MQ=255)
      tGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGTACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAg                                                                                               <  2:1787414/100‑1 (MQ=255)
       ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAgg                                                                                              >  1:2103993/1‑100 (MQ=255)
       ggTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAgg                                                                                              <  2:3239064/100‑1 (MQ=255)
       ggTCTGTTCCCCTCACCCTAACCCTCGCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAgg                                                                                              <  2:3338261/100‑1 (MQ=255)
          cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGctc                                                                                           <  1:2406563/100‑1 (MQ=255)
          cTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctc                                                                                           <  1:1431509/100‑1 (MQ=255)
           tgtgCCCCTTACCCTAACCCTCTCCCCAGATGAGCGAGAGGATTGTCCGTGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctca                                                                                          <  2:2961450/100‑2 (MQ=25)
           tgtgCCCCTCACCCTAAGCATCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACGTTGTCATCAGTCTGCGCCTGGCAATGGCAGGctca                                                                                          <  2:2340331/100‑2 (MQ=21)
           tGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctca                                                                                          <  1:1025648/100‑2 (MQ=255)
            gtgcacttCACGCTAGCCCTCTGCCAAGAGGGCCGATGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGGCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  2:2210144/93‑3 (MQ=255)
            gTTCCCCTCACCCTAGCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCAGAGCCTGCCAATGGCAGGctcag                                                                                         <  1:3768982/100‑3 (MQ=37)
            gTTCCCCTCACCCTAACCCTCTCCTCGGATGTGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  2:1894652/100‑3 (MQ=37)
            gTTCCCCTCACCCTAACCCTCTCCCCGGAGGGTCGAGGGGAGTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  2:3197431/100‑3 (MQ=37)
            gTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  1:1779680/100‑3 (MQ=255)
            gTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  1:3377421/100‑3 (MQ=255)
            gTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  1:1545473/100‑3 (MQ=255)
            gTTACCCTCACCCTAACCCTCTTCCCGGAGGGGCGAGGGGACTGGCCGGGCACATTTGTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcag                                                                                         <  2:2264207/100‑3 (MQ=25)
             ttCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTgg                                                                                        <  1:878686/100‑1 (MQ=255)
              tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                       <  1:3095865/100‑1 (MQ=255)
              tCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgt                                                                                       <  1:61350/100‑1 (MQ=255)
                cccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGgtgt                                                                                     >  2:3554923/1‑100 (MQ=255)
                 ccTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTc                                                                                    <  1:244566/100‑1 (MQ=255)
                  cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTcc                                                                                   <  1:3577338/100‑1 (MQ=255)
                  cTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCAATGGCAGGctcagac                                                                                       <  2:3338807/96‑5 (MQ=255)
                     aCCCTAACCCTCTCCCCGGAGGGGGGAGGGGATTGTCCGGGTACATTTTTAGACTTTGTCATCAGTCTGATCCTGCCATTGGCAGGCTCTGGTGTCCtt                                                                                 <  2:158547/99‑1 (MQ=38)
                      cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtttt                                                                               <  1:2208957/100‑1 (MQ=255)
                      cccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCtttt                                                                               >  1:2799466/1‑100 (MQ=255)
                       ccTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTa                                                                              <  1:2753322/100‑1 (MQ=255)
                       ccTAACCCTCTCCCCGGAGGGGCGAGGCGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTa                                                                              <  2:2127302/100‑1 (MQ=255)
                        cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTAc                                                                             <  2:34803/100‑1 (MQ=255)
                        cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTAc                                                                             <  1:2026419/100‑1 (MQ=255)
                        cTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTAc                                                                             >  1:2948555/1‑100 (MQ=255)
                         tgaCCCTCTCCCCGGAGGGGCGATGGCACTGTACGGGCACATTTTTAGACTTTGTCCTCAGTCTGAGCCTGCCATTGGCAGGCTCTGTTGTCCTTTTAc                                                                             <  2:3088223/97‑1 (MQ=255)
                          aaCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGCTGTCCTTTTACg                                                                            <  1:2504183/99‑1 (MQ=255)
                           aCCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCt                                                                          <  2:125593/100‑1 (MQ=255)
                             ccTCTCCCCGGATGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAc                                                                        <  1:2379943/100‑1 (MQ=255)
                              ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAcc                                                                       <  1:1584447/100‑1 (MQ=255)
                              ctctCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTAcc                                                                       <  1:3379065/100‑1 (MQ=255)
                               tctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                      <  1:139239/100‑1 (MQ=255)
                               tctcCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCa                                                                      <  1:1346198/100‑1 (MQ=255)
                                  ccccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                   <  2:2844965/100‑1 (MQ=255)
                                  ccccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGc                                                                   <  1:1895294/100‑1 (MQ=255)
                                   cccGGAGGGGCGAGGGGGCTGTCCGGGCACATTTTTAGACTGTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGGCCTTTTACGCTACCATGCt                                                                  >  2:2276639/1‑100 (MQ=255)
                                   cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                  >  1:1837992/1‑100 (MQ=255)
                                   cccGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCt                                                                  <  1:684009/100‑1 (MQ=255)
                                      ggAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCtaat                                                               >  1:380341/1‑100 (MQ=255)
                                           ggCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATATTCa                                                          <  1:2651440/100‑1 (MQ=255)
                                            gCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                         <  1:1230677/100‑1 (MQ=255)
                                            gCGAGGGGACTGTCCGGGCACATTTTTAGACTTGGTCATCAGTCTGAGCCTGCCATTTGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                         <  2:500714/100‑1 (MQ=255)
                                             cGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAg                                                         <  2:2428419/99‑1 (MQ=255)
                                              gAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGca                                                       >  2:1661365/1‑100 (MQ=255)
                                                ggggACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTCGACGCGACCATGGTAATAATCAGcaca                                                     <  2:1656730/100‑1 (MQ=255)
                                                    aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataa                                                 <  1:2962083/100‑1 (MQ=255)
                                                    aCTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataa                                                 >  1:3223274/1‑100 (MQ=255)
                                                     cTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACaataat                                                >  2:2238489/1‑100 (MQ=255)
                                                       gTCCGGGCACATATTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCCTGCTAATAACCAGAACAATAATCa                                              >  2:3025557/1‑100 (MQ=255)
                                                           gggCACATTTTTAGACTTTGTCATCAGTCTGAGGCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGccc                                          <  1:792028/100‑1 (MQ=255)
                                                              cacaTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAAc                                       >  1:31266/1‑100 (MQ=255)
                                                                 aTTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCAc                                    <  1:2968490/100‑1 (MQ=255)
                                                                   ttttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACgg                                  <  1:1853104/100‑1 (MQ=255)
                                                                     ttAGGCTTTGTCATCAGTCTGAACCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATTCTAATAATCAGCACTATAATCAGCCCAAcc                                      >  2:3425998/1‑94 (MQ=255)
                                                                     ttAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAg                                <  1:500752/100‑1 (MQ=255)
                                                                      tAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGt                               <  1:424497/100‑1 (MQ=255)
                                                                        gACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                             <  1:3515270/100‑1 (MQ=255)
                                                                        gACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTg                             <  2:637110/100‑1 (MQ=255)
                                                                         aCTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGa                            <  2:3592068/100‑1 (MQ=255)
                                                                          cTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGAc                           <  1:1181946/100‑1 (MQ=255)
                                                                            gatgacaaaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCa                         >  2:2729266/9‑100 (MQ=255)
                                                                             tGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAg                        <  1:1743237/100‑1 (MQ=255)
                                                                              gTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGc                       >  2:1916344/1‑100 (MQ=255)
                                                                              gTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGc                       <  1:3624521/100‑1 (MQ=255)
                                                                               tcatcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCGCAACCACGGAGTTGACCAGCt                      <  2:1361651/100‑1 (MQ=255)
                                                                                 atcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAGCCACGGAGTTGACCAGCTcc                    >  1:2080702/1‑100 (MQ=255)
                                                                                  tcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCca                   >  1:1938938/1‑100 (MQ=255)
                                                                                  tcaGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCca                   <  2:1131828/100‑1 (MQ=255)
                                                                                   caGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcag                  >  1:2230279/1‑100 (MQ=255)
                                                                                    caccagAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagc                 >  1:1073370/6‑100 (MQ=255)
                                                                                    caccagAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagc                 >  1:2149923/6‑100 (MQ=255)
                                                                                    aGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGATAATAATAAGCACAATAATCAGAACAACCACGGAGATGACCAGCTCaagc                 >  2:936122/1‑100 (MQ=255)
                                                                                     gTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAACAATCAGCCCAACCACGGAGTTGACCAGCTCcagca                >  1:255615/1‑100 (MQ=255)
                                                                                      tCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCcagcag               >  1:1043292/1‑100 (MQ=255)
                                                                                       cTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGa              >  2:2529233/1‑100 (MQ=255)
                                                                                        agAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAc             >  2:3663363/2‑100 (MQ=255)
                                                                                         gAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcc            >  1:3254884/1‑100 (MQ=255)
                                                                                         gAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcc            <  2:1746262/100‑1 (MQ=255)
                                                                                         gAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCTAACCACGGAGTTGACCAGCTCCAGCAGAcc            >  1:3312746/1‑100 (MQ=255)
                                                                                         gAGCCTGCCAATGGCAGGCTCTGGTGTCCTTTTACGCTACCAAGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACAAGCTACAGCAGAc             >  2:785919/1‑99 (MQ=255)
                                                                                          aGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAccc           >  2:332550/1‑100 (MQ=255)
                                                                                           gCCTGCCATTGGTAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc          >  1:2494558/1‑100 (MQ=255)
                                                                                           gCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc          >  2:3315532/1‑100 (MQ=255)
                                                                                           gCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGAcccc          <  1:1010670/100‑1 (MQ=255)
                                                                                            ccTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCa         <  1:3286816/100‑1 (MQ=255)
                                                                                             cTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg        >  2:2397870/1‑100 (MQ=255)
                                                                                             cTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAg        >  1:3034522/1‑100 (MQ=255)
                                                                                                ccATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGtt     >  1:2179183/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTTTGCTGGTCTGTTCCCCTCACCCTAACCCTCTCCCCGGAGGGGCGAGGGGACTGTCCGGGCACATTTTTAGACTTTGTCATCAGTCTGAGCCTGCCATTGGCAGGCTCTGGTGTCCTTTTACGCTACCATGCTAATAATCAGCACAATAATCAGCCCAACCACGGAGTTGACCAGCTCCAGCAGACCCCAGGTTTTC  >  NC_000913/4549833‑4550031

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: