Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 16.9% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC16.9% 382.0 / 73.2 153intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (67/60);  new base C (13/13);  total (80/73)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

Only 100 of 173 total aligned reads displayed.

ACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCT  >  NC_000913/4087993‑4088190
                                                                                                  |                                                                                                   
aCTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                    <  1:2965696/100‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTTATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                    >  2:317466/1‑99 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  1:537633/100‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   >  1:3824177/1‑100 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   >  1:1949493/1‑100 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  2:2068018/100‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                    <  2:1241202/99‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgccg                                                                                                   <  1:2168779/100‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcgccg                                                                                                   <  1:3462360/100‑1 (MQ=255)
   ntttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGgcgc                                                                                                  >  2:614828/2‑99 (MQ=255)
   ctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                 <  1:890377/100‑1 (MQ=255)
      tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                              <  2:177021/100‑1 (MQ=255)
      tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                              <  1:371361/100‑1 (MQ=255)
      tCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATcc                                                                                              >  1:3531863/1‑100 (MQ=255)
       cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                             >  2:974054/1‑100 (MQ=255)
       cccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCg                                                                                             >  1:2013696/1‑100 (MQ=255)
        ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCg                                                                                             >  2:1605580/1‑99 (MQ=255)
        ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGa                                                                                            <  2:2866679/100‑1 (MQ=255)
        ccGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGa                                                                                            >  1:3554692/1‑100 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                           <  1:1068484/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:3872579/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  1:30514/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:2316904/100‑1 (MQ=255)
          gagaGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                          <  1:705269/100‑1 (MQ=255)
          gagaGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:2614992/99‑1 (MQ=255)
           agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  1:1412618/100‑1 (MQ=255)
           agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  2:2013686/100‑1 (MQ=255)
           aaagCATGACGCGAGTACAACGACATGGTTTATGACTATAATAGCCTTCTTTAAACGTTGTTGCGGATAAATGCAGTGTCTGATgccgcg                                                                                                   >  2:1311283/3‑90 (MQ=255)
              gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                                      <  2:3914440/100‑1 (MQ=255)
               cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGtt                                                                                     >  2:2513598/1‑100 (MQ=255)
               cATTACACGAGTACAACGAAATTGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                                      <  2:3531837/99‑1 (MQ=255)
                  tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                  <  1:1436978/100‑1 (MQ=255)
                  tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                  >  1:390310/1‑100 (MQ=255)
                   aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                 <  2:1593246/100‑1 (MQ=255)
                    cgcgAGTACAACGAAATGGTTGATGACTAAAAGAGCCCGCGTTAAACGTGTTTGAGGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                  <  2:3568199/98‑1 (MQ=255)
                     gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                 >  1:791361/1‑98 (MQ=255)
                        aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                            <  2:1377494/100‑1 (MQ=255)
                        aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                            >  1:3408132/1‑100 (MQ=255)
                        aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCAGATGCGGCACATCCGACACAGTTCACCCCtt                                                                             >  2:3523809/1‑99 (MQ=255)
                         gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATACATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                           <  2:300753/100‑1 (MQ=255)
                         gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                           >  1:3694243/1‑100 (MQ=255)
                           aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                         <  2:672512/100‑1 (MQ=255)
                           aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGAAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                         >  1:2614763/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGCTGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  2:2195290/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  1:2633007/1‑100 (MQ=255)
                             aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                       <  2:3694215/100‑1 (MQ=255)
                             aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                       >  2:3428023/1‑100 (MQ=255)
                                  aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACt                                                                  <  1:2333494/100‑1 (MQ=255)
                                         gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                           >  2:2692806/1‑100 (MQ=255)
                                         gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                           >  1:724435/1‑100 (MQ=255)
                                             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                       >  1:649187/1‑100 (MQ=255)
                                             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                       >  1:650075/1‑100 (MQ=255)
                                             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                       >  1:1715894/1‑100 (MQ=255)
                                               tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTc                                                     <  1:3335708/100‑1 (MQ=255)
                                                ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                    >  1:1962171/1‑100 (MQ=255)
                                                ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                    <  2:1538391/100‑1 (MQ=255)
                                                 taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCat                                                   >  2:1915953/1‑100 (MQ=255)
                                                      gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                              >  1:489309/1‑100 (MQ=255)
                                                       ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                             >  1:1243848/1‑100 (MQ=255)
                                                         tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                           <  2:1447129/100‑1 (MQ=255)
                                                           cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                         >  2:1084537/1‑100 (MQ=255)
                                                            gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                        >  1:2307920/1‑100 (MQ=255)
                                                             ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                       >  1:3204372/1‑100 (MQ=255)
                                                              tAAACGTTTTTGATGCTAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                                      >  2:507886/1‑100 (MQ=255)
                                                               aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                     <  2:3077837/100‑1 (MQ=255)
                                                                aaCGTTTTTGATTATAAATGCAGTGTCGGATGCGGCGCATCCGAGACAGTTCACCCCTTGGGAAGTAACTCCATAGATGCCTCATATCCTGTCAACCCgg                                    <  2:1117938/100‑1 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  2:3249162/100‑1 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  2:2260258/100‑1 (MQ=255)
                                                                  cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                  >  2:2821864/1‑100 (MQ=255)
                                                                    tttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGCCACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                       <  1:3120487/93‑1 (MQ=255)
                                                                    tttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                <  1:2692823/100‑1 (MQ=255)
                                                                     ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                               <  1:3685436/100‑1 (MQ=255)
                                                                      tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                              >  1:1820861/1‑100 (MQ=255)
                                                                       ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATACCTCATATCCTGTTAACCCGGCaataa                              >  1:295610/1‑99 (MQ=255)
                                                                       ttGATGATAAATGCAGTGTCGGATGCAGTCCATCCTACACAGTTCACCCCTTTTGAAGTACCTCCATAAATTCCTCATATCCTGTTAACCCGGCaataa                              <  2:3511623/99‑1 (MQ=255)
                                                                         gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                           >  1:2707117/1‑100 (MQ=255)
                                                                            gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                        >  1:2191736/1‑100 (MQ=255)
                                                                            gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                         >  2:525717/1‑99 (MQ=255)
                                                                              tAAATGCAGTGTCTGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                      >  2:2444236/1‑100 (MQ=255)
                                                                                aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                    <  1:883907/100‑1 (MQ=255)
                                                                                 gaactgtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                   >  1:102089/6‑100 (MQ=255)
                                                                                 aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                   >  2:2015242/1‑100 (MQ=255)
                                                                                  tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                  <  1:3893665/100‑1 (MQ=255)
                                                                                   gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                 >  2:2862966/1‑100 (MQ=255)
                                                                                     aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAATTGCCTCATATCCTGCTAACCCGGCAATAATCTCCGCCGGATGCt               <  2:1435657/100‑1 (MQ=255)
                                                                                      gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg              >  2:66831/1‑100 (MQ=255)
                                                                                      gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGACAATAATCTCCGCCGGATGCTg              >  1:2420241/1‑100 (MQ=255)
                                                                                       tgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa             >  1:2507400/1‑100 (MQ=255)
                                                                                       tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa             >  1:1666179/1‑100 (MQ=255)
                                                                                       tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa             >  1:54892/1‑100 (MQ=255)
                                                                                        gtCGGCTGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCAGATGc                >  2:877923/1‑96 (MQ=255)
                                                                                         tCGGATGTGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca           <  2:3675984/100‑1 (MQ=255)
                                                                                         tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGGAACTGCAGAAATGCCTCATGTCCTGTGAACCCGGCAATAATCTCCGCCGGATGCTGac            >  1:179426/1‑99 (MQ=255)
                                                                                          cGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac          >  2:438623/1‑100 (MQ=255)
                                                                                           ggATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca         >  1:3643366/1‑100 (MQ=255)
                                                                                            gATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACaa        <  1:1506243/100‑1 (MQ=255)
                                                                                              tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata      >  1:222308/1‑100 (MQ=255)
                                                                                               gcgccgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTAATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGGCACAatat     <  2:11655/100‑1 (MQ=255)
                                                                                                  gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt  >  1:204926/1‑100 (MQ=255)
                                                                                                  gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt  >  1:1163343/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCT  >  NC_000913/4087993‑4088190

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: