Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 18.2% 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,0910GC18.2% 349.5 / 69.6 137intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (58/54);  new base C (12/13);  total (70/67)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 149 total aligned reads displayed.

ACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCT  >  NC_000913/4087993‑4088190
                                                                                                  |                                                                                                   
aCTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                    >  1:2573228/1‑100 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  2:2688663/100‑1 (MQ=255)
 ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  1:2786908/100‑1 (MQ=255)
   ctTTCCCGAGAGCATTACGCGAGTACAACTAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTCTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                 <  1:2130279/100‑1 (MQ=255)
    tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAt                                                                                                <  1:3539531/100‑1 (MQ=255)
     ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATc                                                                                               >  2:2870427/1‑100 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                           <  2:216884/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                           <  2:3483546/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:250819/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  1:3461135/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:3363195/100‑1 (MQ=255)
         cGAGAGCATTACGCGAGTACAACGAAATGCTAGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                           <  1:2934532/100‑1 (MQ=255)
          gagaGTAGTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  1:702593/99‑1 (MQ=255)
          gagaGCATTACGCGAGTACAACGAAATGGGTGTTGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  1:676866/99‑1 (MQ=255)
           agagCATTATGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  2:18277/100‑1 (MQ=255)
            gagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGAcaca                                                                                        <  2:2017514/100‑1 (MQ=255)
            gagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacact                                                                                        <  2:1927800/100‑2 (MQ=255)
              gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacactgc                                                                                      <  1:2123367/100‑4 (MQ=255)
                aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                                    <  1:3376040/100‑1 (MQ=255)
                aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTc                                                                                    >  2:1376723/1‑100 (MQ=255)
                   aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                 >  2:2352852/1‑100 (MQ=255)
                   aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                 >  1:1386199/1‑100 (MQ=255)
                    cgcgAGTTCAACGAAAGGCTTGATGACTATAATAGCCGGGGTGAAATGTGTTTGATGATAAACGCAGTGTCGGATGCGGAGCATCCAACAAAGTTCaac                                                                                 >  1:2763151/1‑97 (MQ=255)
                    cgcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAccc                                                                                >  2:876559/1‑100 (MQ=255)
                      cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCt                                                                              >  2:3217789/1‑100 (MQ=255)
                        aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                            <  2:2617633/100‑1 (MQ=255)
                         gTACAACGAAATGGTTGATGACTATAGTAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                           >  1:3214426/1‑100 (MQ=255)
                         gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtgt                                                                           >  2:4081/1‑100 (MQ=255)
                           aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                         >  2:1529842/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  2:1251567/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  2:227649/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  2:2982036/1‑100 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        <  2:3494024/100‑1 (MQ=255)
                            cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                        >  2:902151/1‑100 (MQ=255)
                             aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                       <  1:134716/100‑1 (MQ=255)
                              aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                      >  2:1966132/1‑100 (MQ=255)
                              aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                      >  1:3064438/1‑100 (MQ=255)
                              aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                      >  2:3046032/1‑100 (MQ=255)
                                 aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAAc                                                                   >  2:758713/1‑100 (MQ=255)
                                  aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACt                                                                  <  1:1197610/100‑1 (MQ=255)
                                        ngatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaa                                                            >  2:3520726/2‑100 (MQ=255)
                                           tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                         <  2:26731/100‑1 (MQ=255)
                                             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                       <  2:1096891/100‑1 (MQ=255)
                                             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATCCGGCCCATCCGACACAGTTCACCCCTTGTTAAGTAACTCCATAAATGcc                                                       >  1:273233/1‑100 (MQ=255)
                                              cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATTCGGTGCATCCGACACAGTTTAACCCTTGTGATGTAACTCCGTAAACGcc                                                       >  1:603891/1‑99 (MQ=255)
                                              cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCt                                                      >  1:1531970/1‑100 (MQ=255)
                                              cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAACGCCt                                                      >  2:1491379/1‑100 (MQ=255)
                                                ataataGCCCGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCa                                                    <  2:2892937/100‑1 (MQ=255)
                                                  aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                                                  >  2:1726618/1‑100 (MQ=255)
                                                  aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCata                                                  >  2:1386801/1‑100 (MQ=255)
                                                     aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                               <  2:3353290/100‑1 (MQ=255)
                                                     aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                               >  2:687098/1‑100 (MQ=255)
                                                       ccTGCGTTAATCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                             <  1:2337897/100‑1 (MQ=255)
                                                        cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                            <  2:1386206/100‑1 (MQ=255)
                                                         tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACATAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                           <  1:1966122/100‑1 (MQ=255)
                                                            gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                        >  2:2005131/1‑100 (MQ=255)
                                                            gTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                        <  1:779138/100‑1 (MQ=255)
                                                             ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                       <  2:2775246/100‑1 (MQ=255)
                                                             ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                       <  1:3459720/100‑1 (MQ=255)
                                                                aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                    >  2:2374352/1‑100 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  1:1149221/100‑1 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  1:2219195/100‑1 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  1:1389703/100‑1 (MQ=255)
                                                                 aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                   <  1:396467/100‑1 (MQ=255)
                                                                  cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCACCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                  >  2:60373/1‑100 (MQ=255)
                                                                    tttttGATGATAAATGCAGTGTCGGATGCGGTGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                <  2:2631501/100‑1 (MQ=255)
                                                                       ttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                             >  1:1319953/1‑100 (MQ=255)
                                                                         gatgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                           <  2:3143900/100‑1 (MQ=255)
                                                                           tgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                         <  1:1078303/100‑1 (MQ=255)
                                                                              tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATTCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                      >  1:2071520/1‑100 (MQ=255)
                                                                              tAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                      >  2:3255339/1‑100 (MQ=255)
                                                                                aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                    >  2:554018/1‑100 (MQ=255)
                                                                                aaTGCAGTGTCGGAGGCGGCGCAACCGACACAGTTCCCCCGTTGTGGAGGAACTCCATAAATGCCCCCTATGCTGTTAACCACACAATAATCTccgcc                      >  1:1338574/1‑98 (MQ=255)
                                                                                  tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                  >  2:3505400/1‑100 (MQ=255)
                                                                                  tGCAGTGGCGGATGCGGCGCATCCGACATAGTGCCCCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTGAACCCGGCAATAATCTccgccg                     <  1:3526290/97‑1 (MQ=255)
                                                                                   gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                 >  1:1653245/1‑100 (MQ=255)
                                                                                    ctgtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGGAGTAACTCCATAAATGCCTCTTATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                >  2:2481661/3‑100 (MQ=39)
                                                                                     aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt               <  2:1053202/100‑1 (MQ=255)
                                                                                      gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg              <  2:3250058/100‑1 (MQ=255)
                                                                                      gtgtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg              >  2:1485824/1‑100 (MQ=255)
                                                                                      gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg              >  2:2151423/1‑100 (MQ=255)
                                                                                       tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa             >  1:1232556/1‑100 (MQ=255)
                                                                                       tgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCATGTTAACCCGGCAATAATCTCCGCCGGATGCTGa             >  2:857990/1‑100 (MQ=255)
                                                                                        gtCGGTTGCGGCGCATCCACCACAGTTCCCCCTTTGTGAACTAACTCCATAAATGGCTCAGATCCTGTCAACCCGGCAATAATCTCCGCCAGATGCTGa             <  1:447567/99‑1 (MQ=255)
                                                                                        gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  2:279986/1‑100 (MQ=255)
                                                                                        gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            <  1:3212469/100‑1 (MQ=255)
                                                                                        gtCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  2:1856190/1‑100 (MQ=255)
                                                                                        gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  1:2717086/1‑100 (MQ=255)
                                                                                        gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  1:2929201/1‑100 (MQ=255)
                                                                                        gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  2:1589607/1‑100 (MQ=255)
                                                                                         tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca           <  1:544673/100‑1 (MQ=255)
                                                                                         tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac            >  2:1261009/1‑99 (MQ=255)
                                                                                           ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca         >  2:1153869/1‑100 (MQ=255)
                                                                                             acgcggcgCACCAGACACAGTTTCCCCCTTGGGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat       <  1:3187459/98‑1 (MQ=255)
                                                                                              tGCGGCGCATCCGACACAGTTCACCTCATGTGAAGTAACTCCATAAATGCCTCATATCCTGTTATCCCGGCAATAATCTCCGCCTGATGCTGACACAata      <  1:1680664/100‑1 (MQ=255)
                                                                                               gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat     <  1:253690/100‑1 (MQ=255)
                                                                                                cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg    >  2:2241849/1‑100 (MQ=255)
                                                                                                cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg    >  2:1889865/1‑100 (MQ=255)
                                                                                                 ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc   >  1:1704734/1‑100 (MQ=255)
                                                                                                  gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt  <  1:3373086/100‑1 (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: