Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,180 C→G 44.2% A160P (GCA→CCA)  yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,1800CG44.2% 65.1 / 249.2 165A160P (GCA→CCA) yijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base G (34/39);  total (78/87)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 169 total aligned reads displayed.

TAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGT  >  NC_000913/4137081‑4137279
                                                                                                   |                                                                                                   
tAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGc                                                                                                     >  1:1731372/1‑100 (MQ=255)
 aaaCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGCCGTACCGTCGCGGGCGAAGCGATCTGAAACCACCCCAATATGTGcc                                                                                                    >  1:1035728/1‑100 (MQ=255)
   aCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAg                                                                                                  <  1:3245381/100‑1 (MQ=255)
   aCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                                                                  >  1:2300053/1‑100 (MQ=255)
   aCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTAGCACCATCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                                                                  <  2:788356/100‑1 (MQ=255)
     aGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGcc                                                                                                <  2:177164/100‑1 (MQ=255)
      gCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGccc                                                                                               >  2:1213196/1‑100 (MQ=255)
       cacaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCg                                                                                              >  2:2546090/1‑100 (MQ=255)
       cacaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCg                                                                                              >  2:3310940/1‑100 (MQ=255)
        acaTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGt                                                                                             <  2:506346/100‑1 (MQ=255)
         caTCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGtt                                                                                            >  1:2209691/1‑100 (MQ=255)
         caTCTTCTTCCTCCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGtt                                                                                            >  2:2801697/1‑100 (MQ=255)
           tcttctTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTAt                                                                                          <  2:3798774/100‑1 (MQ=255)
            cttcttCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATc                                                                                         >  1:1165571/1‑100 (MQ=255)
             ttcttcCTGCTCGCCTTCGCCAATGTTATGTATCACCAGCGGCGTATCGTCGCTGGCGAAGCTATCTGAAACCACCCCAATATGTGGCAGCCCGTTATcc                                                                                        <  1:861032/100‑1 (MQ=255)
             ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGCAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATcc                                                                                        <  1:2335194/100‑1 (MQ=255)
             ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATcc                                                                                        <  1:4232887/100‑1 (MQ=255)
             ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATcc                                                                                        <  1:1169414/100‑1 (MQ=255)
             ttcttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATcc                                                                                        <  1:3175146/100‑1 (MQ=255)
                ttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGt                                                                                     >  2:3178043/1‑100 (MQ=255)
                 tcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGACGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTc                                                                                    <  1:456127/100‑1 (MQ=255)
                  ccTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                                                                   >  1:1930671/1‑100 (MQ=255)
                  ccTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                                                                   <  1:4241801/100‑1 (MQ=255)
                   cTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGc                                                                                  >  1:733425/1‑100 (MQ=255)
                       gcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAgg                                                                              <  1:1788923/100‑1 (MQ=255)
                         gcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGgag                                                                            >  2:359792/1‑100 (MQ=255)
                         gcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                                                            >  1:438168/1‑100 (MQ=255)
                           cTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCTATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGgag                                                                            >  1:2866322/1‑98 (MQ=255)
                           cTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGAc                                                                          >  1:2696744/1‑100 (MQ=255)
                            tGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACa                                                                         <  2:3488668/100‑1 (MQ=255)
                            tGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                                                         <  2:1140944/100‑1 (MQ=255)
                              cgcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAat                                                                       >  2:2970610/1‑100 (MQ=255)
                                 cAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATc                                                                    >  1:245466/1‑100 (MQ=255)
                                   aTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGc                                                                  <  1:746183/100‑1 (MQ=255)
                                   aTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGc                                                                  <  2:1004366/100‑1 (MQ=255)
                                       tATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                                                               >  1:2536571/1‑99 (MQ=255)
                                           tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                          >  2:1328536/1‑100 (MQ=255)
                                           tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                          >  2:971581/1‑100 (MQ=255)
                                           tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                          <  2:4055950/100‑1 (MQ=255)
                                           tATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                                                          >  1:2047800/1‑100 (MQ=255)
                                              cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  1:266228/100‑1 (MQ=255)
                                              cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  2:1955286/100‑1 (MQ=255)
                                              cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  2:627225/100‑1 (MQ=255)
                                              cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg                                                       <  2:424380/100‑1 (MQ=255)
                                                ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc                                                     <  2:3695140/100‑1 (MQ=255)
                                                ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc                                                     >  2:3771110/1‑100 (MQ=255)
                                                 cAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa                                                    >  1:782334/1‑100 (MQ=255)
                                                   gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACt                                                  <  1:359790/100‑1 (MQ=255)
                                                    cggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTg                                                 <  2:862014/100‑1 (MQ=255)
                                                     ggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg                                                <  2:1487094/100‑1 (MQ=255)
                                                      gcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg                                               <  1:145068/100‑1 (MQ=255)
                                                        gTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt                                             <  2:2696241/100‑1 (MQ=255)
                                                        gTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt                                             <  2:2290105/100‑1 (MQ=255)
                                                         tACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGtt                                            >  1:3627229/1‑100 (MQ=255)
                                                          aCCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTc                                           >  1:3467063/1‑100 (MQ=255)
                                                           ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt                                          <  2:113156/100‑1 (MQ=255)
                                                           ccGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt                                          >  1:1391802/1‑100 (MQ=255)
                                                            cGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt                                         <  2:4046745/100‑1 (MQ=255)
                                                              tCGCAGGCGAAGCCATCTGAAACCACCCCAATATATGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGAGTTCTTg                                        >  1:3776116/1‑99 (MQ=255)
                                                               cgcgGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCt                                      <  2:2118829/100‑1 (MQ=255)
                                                               cgcgGGCGAAGCCATCTGAAACCACCCCAAGATATGCCAGCCCGTTATCCAGTCGCCAGGAGACTATATCTCCCGCTTGATATTCACTGAGGTTCtt                                         >  1:2801181/1‑97 (MQ=255)
                                                                 cgGGCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtgg                                    <  2:182376/100‑1 (MQ=255)
                                                                  gggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggt                                   >  2:854077/1‑100 (MQ=255)
                                                                   ggCGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtg                                  <  2:942197/100‑1 (MQ=255)
                                                                    gCGAAGCCATCTGAAACCACCCCAATATGTGCCGGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtgg                                 >  1:2782914/1‑100 (MQ=255)
                                                                    gCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtgg                                 <  1:3454539/100‑1 (MQ=255)
                                                                    gCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCtggtgg                                 <  1:124679/100‑1 (MQ=255)
                                                                     cGAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTggg                                <  1:1286658/100‑1 (MQ=255)
                                                                     cGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTggg                                <  2:3255843/100‑1 (MQ=255)
                                                                     cGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTggg                                >  2:1519844/1‑100 (MQ=255)
                                                                      gAAGCCATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTgggg                               <  1:2970598/100‑1 (MQ=255)
                                                                         gCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgt                            >  1:3215434/1‑100 (MQ=255)
                                                                          ccATCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCgtg                           >  1:1130207/1‑100 (MQ=255)
                                                                            aTCTGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGtt                         >  2:3260997/1‑100 (MQ=255)
                                                                               tGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                      >  2:3541366/1‑100 (MQ=255)
                                                                               tGAAACCACCCCAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                      <  1:1878991/100‑1 (MQ=255)
                                                                               tGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                      <  2:4188877/100‑1 (MQ=255)
                                                                               tGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTa                      <  1:917910/100‑1 (MQ=255)
                                                                                 aaaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATc                    <  1:4104273/100‑1 (MQ=255)
                                                                                  aaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATc                    >  2:1313980/1‑99 (MQ=255)
                                                                                   accaccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCAt                  <  2:1663963/100‑1 (MQ=255)
                                                                                     caccCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATgg                >  1:1275542/1‑100 (MQ=255)
                                                                                       ccccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:1494746/1‑100 (MQ=255)
                                                                                       ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              >  1:1049560/1‑100 (MQ=255)
                                                                                       ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg              <  1:959248/100‑1 (MQ=255)
                                                                                       ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggc               >  2:2006730/1‑99 (MQ=255)
                                                                                        cccAATATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcgg             <  1:1955326/100‑1 (MQ=255)
                                                                                         ccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcggc            >  2:2153455/1‑100 (MQ=255)
                                                                                         ccAATATGTGCCAGCCCGTTATCCAGTCGCCAGCAAACAATATCGCCCGCATGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggc               >  1:794625/1‑97 (MQ=255)
                                                                                          cAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCt           >  2:415937/1‑100 (MQ=255)
                                                                                           aaTATGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTa          <  2:3970339/100‑1 (MQ=255)
                                                                                             tatGTGGCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTaaa        >  2:2926866/1‑100 (MQ=255)
                                                                                              atGTGGCAGCCCGATATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc       <  2:2478954/100‑1 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      <  1:1297348/100‑1 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      >  1:615266/1‑100 (MQ=255)
                                                                                               tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAcc      <  2:1930675/100‑1 (MQ=255)
                                                                                                gtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa     <  2:1453033/100‑1 (MQ=255)
                                                                                                 tgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAg    <  2:174441/100‑1 (MQ=255)
                                                                                                  gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAgg   >  2:1703697/1‑100 (MQ=255)
                                                                                                   ccAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGt  <  2:602393/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGT  >  NC_000913/4137081‑4137279

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑