Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,079,004 G→C 8.2% intergenic (+73/+28) yoeF → / ← yeeX pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,079,0040GC 8.2% 210.0 / 8.8 61intergenic (+73/+28)yoeF/yeeXpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base C (2/3);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATA  >  NC_000913/2078856‑2079148
                                                                                                                                                    |                                                                                                                                                
aGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCg                                                                                                                                                  >  2:523524/1‑149 (MQ=255)
 gTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                 <  1:405597/149‑1 (MQ=255)
   aaCCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGc                                                                                                                                               <  3:514526/149‑1 (MQ=255)
    aCCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                                                                                                                                              >  1:589656/1‑149 (MQ=255)
           cAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCa                                                                                                                                       <  4:405393/149‑1 (MQ=255)
           cAGTTAATAACAATGAAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCa                                                                                                                                       <  3:488480/149‑1 (MQ=255)
            aGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAt                                                                                                                                      <  3:183995/149‑1 (MQ=255)
            aGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCAt                                                                                                                                      <  3:346712/149‑1 (MQ=255)
             gTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAtt                                                                                                                                     >  3:547109/1‑149 (MQ=255)
              ttAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                 <  2:286700/136‑1 (MQ=255)
              ttAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCgg                                                                                                                                                 >  1:286700/1‑136 (MQ=255)
                aataaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGaa                                                                                                                                  <  2:556737/149‑1 (MQ=255)
                   aaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAt                                                                                                                                      <  3:1448/142‑1 (MQ=255)
                   aaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCAt                                                                                                                                      >  4:1448/1‑142 (MQ=255)
                   aaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATaa                                                                                                                               >  3:266530/1‑149 (MQ=255)
                       aTGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACg                                                                                                                           >  1:184604/1‑149 (MQ=255)
                       aTGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACg                                                                                                                           <  1:19649/149‑1 (MQ=255)
                        tGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACgg                                                                                                                          <  4:266530/149‑1 (MQ=255)
                        tGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACgg                                                                                                                          <  2:327635/149‑1 (MQ=255)
                         ggAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACggg                                                                                                                         <  2:232183/149‑1 (MQ=255)
                              tACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGACTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTa                                                                                                                    >  2:394278/1‑149 (MQ=255)
                                cAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAct                                                                                                                  <  2:547842/149‑1 (MQ=255)
                                  gAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtc                                                                                                                <  2:210613/149‑1 (MQ=255)
                                  gAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtc                                                                                                                >  3:513036/1‑149 (MQ=255)
                                             cTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAc                                                                                                                   <  1:253737/135‑1 (MQ=255)
                                             cTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAc                                                                                                                   >  2:253737/1‑135 (MQ=255)
                                                  tgcATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCa                                                                                                >  1:403465/1‑149 (MQ=255)
                                                      tACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTc                                                                                            <  4:208709/149‑1 (MQ=255)
                                                         agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGTGTTGCCGGCCAACCCTCTTTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGcc                                                                                         >  4:349513/1‑149 (MQ=255)
                                                         agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttt                                                                                                    >  3:421058/1‑138 (MQ=255)
                                                         agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGtat                                                                                                    <  4:421058/138‑3 (MQ=255)
                                                             gaagGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATa                                                                                     <  2:184604/149‑1 (MQ=255)
                                                             gaagGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATa                                                                                     <  4:513036/149‑1 (MQ=255)
                                                                   tgttgtTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                         >  2:202410/1‑123 (MQ=255)
                                                                   tgttgtTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                         <  1:202410/123‑1 (MQ=255)
                                                                    gttgttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCa                                                                              <  4:547109/149‑1 (MQ=255)
                                                                    gttgttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCa                                                                              <  1:523524/149‑1 (MQ=255)
                                                                     ttgttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAg                                                                             >  4:114757/1‑149 (MQ=255)
                                                                     ttgttAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAg                                                                             >  3:71287/1‑149 (MQ=255)
                                                                           aTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCttttt                                                                       >  1:251122/1‑149 (MQ=255)
                                                                           aTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAAcc                                                                                                                                           >  4:33880/1‑81 (MQ=255)
                                                                           aTTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAAcc                                                                                                                                           <  3:33880/81‑1 (MQ=255)
                                                                            tttCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTg                                                                      >  1:435992/1‑149 (MQ=255)
                                                                              tCGATTGGTTCTCGGATGAGAATGATGTCGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGga                                                                    <  1:166734/149‑1 (MQ=255)
                                                                                  ttGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGata                                                                >  2:440765/1‑149 (MQ=255)
                                                                                           ggATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGCTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTCGAGATATCGCGGCGt                                                       >  2:173440/1‑149 (MQ=255)
                                                                                             aTGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTc                                                     >  4:373292/1‑149 (MQ=255)
                                                                                                 gaATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTg                                                 <  3:373292/149‑1 (MQ=255)
                                                                                                 gaATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTg                                                 >  4:185351/1‑149 (MQ=255)
                                                                                                                tACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCGTAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAttt                                  >  4:181266/1‑149 (MQ=255)
                                                                                                                    caGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtga                              <  4:268701/149‑1 (MQ=255)
                                                                                                                    caGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtga                              <  1:173440/149‑1 (MQ=255)
                                                                                                                    caGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtga                              <  1:440765/149‑1 (MQ=255)
                                                                                                                     aGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtgat                             <  2:435992/149‑1 (MQ=255)
                                                                                                                        tAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTgatgat                          >  4:500971/1‑149 (MQ=255)
                                                                                                                         aaGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGatgatg                         <  1:411411/149‑1 (MQ=255)
                                                                                                                          aGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGga                                               >  3:411763/1‑126 (MQ=255)
                                                                                                                          aGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGga                                               <  4:411763/126‑1 (MQ=255)
                                                                                                                           gAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCGTAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTa                       <  3:181266/149‑1 (MQ=255)
                                                                                                                           gAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTa                       >  1:226525/1‑149 (MQ=255)
                                                                                                                                        atgagGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGa          >  3:357935/1‑149 (MQ=255)
                                                                                                                                        atgagGGTTGCCCGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGa          >  4:152682/1‑149 (MQ=255)
                                                                                                                                          gagGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATc        <  4:71287/149‑1 (MQ=255)
                                                                                                                                             ggTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTc     >  3:294049/1‑149 (MQ=255)
                                                                                                                                                tGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATa  <  3:584514/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATA  >  NC_000913/2078856‑2079148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: