Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,572,201 C→T 7.8% intergenic (‑71/+21) yicG ← / ← orn inner membrane protein/oligoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,572,2010CT 7.8% 224.2 / 6.1 77intergenic (‑71/+21)yicG/orninner membrane protein/oligoribonuclease
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (4/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

CATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGG  >  NC_002947/5572053‑5572346
                                                                                                                                                    |                                                                                                                                                 
cctGCCGCGCCGCCCGGCAGACAGCCCGCCAGTCATGGCTTCGGCGGTAATGTCGATAACGTCGCGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATAc                                                                                                                                                   <  8:92604/147‑1 (MQ=255)
 aTGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGCGGGCGCAGATAc                                                                                                                                                   >  4:320900/1‑148 (MQ=255)
  tGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGTGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                 <  3:105346/149‑1 (MQ=255)
   gCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                 >  1:238839/1‑148 (MQ=255)
   gCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCACTCAAAACAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACaa                                                                                                                                                 >  1:315358/1‑148 (MQ=255)
      gcgACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACCACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATc                                                                                                                                             >  7:299988/1‑149 (MQ=255)
      gcgACGCCCGGCAGACAGCGCGCCACTCAGGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGCAGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAataac                                                                                                                                             >  2:217967/1‑147 (MQ=255)
        gACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTg                                                                                                                                           >  6:76530/1‑149 (MQ=255)
         aCGCCGGGCAGAATGCGCTTCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCTTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCATGGCACCAAAAGGGGGCGCAGATACAATATCTgc                                                                                                                                          <  2:315358/149‑1 (MQ=255)
         aCGCCCGGCAGATAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgc                                                                                                                                          <  1:131012/149‑1 (MQ=255)
         aCGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgc                                                                                                                                          <  3:340486/149‑1 (MQ=255)
          cGCCCGTCAGACCTCCCCCCCCTCATGGCTTCGGCGGTAATGGCGATACGGTAGAGCATCAACAACATGGGGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAACAGCGGGCGCAGATATTGTATCTgcg                                                                                                                                         <  8:323483/149‑1 (MQ=255)
           gCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATATTGTATCTgcgc                                                                                                                                        <  4:292465/149‑1 (MQ=255)
           gCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgcgc                                                                                                                                        >  2:299575/1‑149 (MQ=255)
           gCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGACCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTgcgc                                                                                                                                        <  7:152138/149‑1 (MQ=255)
              cGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGt                                                                                                                                     <  1:135701/149‑1 (MQ=255)
                gCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGttt                                                                                                                                   <  5:232707/149‑1 (MQ=255)
                 cAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTc                                                                                                                                  >  1:4551/1‑149 (MQ=255)
                   gACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGa                                                                                                                                <  2:4551/149‑1 (MQ=255)
                   gACAGCGCGCCAGTCATGGCTTCGGCGGTAATCGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGa                                                                                                                                >  7:184615/1‑149 (MQ=255)
                    aCAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGAc                                                                                                                               <  5:371367/149‑1 (MQ=255)
                    aCAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGAc                                                                                                                               >  1:71036/1‑149 (MQ=255)
                                   tGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGaa                                                                                                                <  2:72787/149‑1 (MQ=255)
                                       ttCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGc                                                                                                            >  5:302379/1‑149 (MQ=255)
                                         cggcggTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTtga                                                                                                                    <  5:356709/139‑1 (MQ=255)
                                         cggcggTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTtga                                                                                                                    >  6:356709/1‑139 (MQ=255)
                                             ggTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGa                                                                                                      <  6:302379/149‑1 (MQ=255)
                                                aaTGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAg                                                                                                   >  6:198228/1‑149 (MQ=255)
                                                    gCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGa                                                                                                <  7:224992/148‑1 (MQ=255)
                                                     cGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACg                                                                                              <  4:71038/149‑1 (MQ=255)
                                                             tAGAGCATCAACAACATGGTGCGGGTCCTTGCAGGAAAGGCTCGCAGTCTAACCAGTTGCCCAATTCACCAAAAGGGGTCGCAGATACAATATCTGCGCCGTTTCGACTCACATCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgt                                                                                      <  3:238428/149‑2 (MQ=255)
                                                             tAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCgg                                                                                      <  6:129899/149‑1 (MQ=255)
                                                                    tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATa                                                                                                    <  8:170251/128‑1 (MQ=255)
                                                                    tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATa                                                                                                    >  7:170251/1‑128 (MQ=255)
                                                                    tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcga                                                                               >  4:167057/1‑149 (MQ=255)
                                                                    tCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGCGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGcgatcga                                                                               <  8:269435/149‑1 (MQ=255)
                                                                      aacaacATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAtt                                                                             <  2:71036/149‑1 (MQ=255)
                                                                           caTGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                        <  5:76530/149‑1 (MQ=255)
                                                                                   gggTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTc                                                                <  8:299988/149‑1 (MQ=255)
                                                                                      tCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAg                                                             >  8:29967/1‑149 (MQ=255)
                                                                                           gCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAgcgc                                                          >  8:333907/1‑147 (MQ=255)
                                                                                            cAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAg                                                       <  7:29967/149‑1 (MQ=255)
                                                                                                aaaGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATAGTCCAGCGCCAGGTgg                                                   <  1:299575/149‑1 (MQ=255)
                                                                                                  aGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCTACAGGGGGCGCAGATTCAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGTCGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGgtgtgt                                                 <  1:217967/149‑1 (MQ=255)
                                                                                                  aGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGt                                                 <  7:333907/149‑1 (MQ=255)
                                                                                                      gcgcAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTg                                             >  5:90000/1‑149 (MQ=255)
                                                                                                            tCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                       >  8:294831/1‑117 (MQ=255)
                                                                                                            tCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGa                                                                       <  7:294831/117‑1 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTAtt                                              >  8:248293/1‑140 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTAtt                                              <  7:248293/140‑1 (MQ=255)
                                                                                                              tAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttctt                                     >  6:106674/1‑149 (MQ=255)
                                                                                                               aaCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTcg                                                                                                         <  1:281762/80‑1 (MQ=255)
                                                                                                               aaCCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTcg                                                                                                         >  2:281762/1‑80 (MQ=255)
                                                                                                                 ccAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAtt                                                                             <  7:128547/106‑1 (MQ=255)
                                                                                                                 ccAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAtt                                                                             >  8:128547/1‑106 (MQ=255)
                                                                                                                 ccAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGaa                                  <  2:238839/149‑1 (MQ=255)
                                                                                                                   aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                        >  3:131414/1‑141 (MQ=255)
                                                                                                                   aGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCctt                                        <  4:131414/141‑1 (MQ=255)
                                                                                                                        cccAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTca                           <  8:184615/149‑1 (MQ=255)
                                                                                                                           aaGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCGTGCTTGAAGATGTcacg                         >  7:87442/1‑148 (MQ=255)
                                                                                                                             ggCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                        <  4:300677/99‑1 (MQ=255)
                                                                                                                             ggCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCgg                                                                        >  3:300677/1‑99 (MQ=255)
                                                                                                                              gCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAg                                 <  3:235698/137‑1 (MQ=255)
                                                                                                                              gCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAg                                 >  4:235698/1‑137 (MQ=255)
                                                                                                                                aCCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTc                   >  7:203256/1‑149 (MQ=255)
                                                                                                                                   aaaaGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                               >  2:44871/1‑54 (MQ=255)
                                                                                                                                   aaaaGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAg                                                                                                               <  1:44871/54‑1 (MQ=255)
                                                                                                                                   aaaaGGGGGCGCAGATACAATATCTGCGCCGATTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATGTCGTCCAGCGGCAGGTGCGTATTGCCATTCTTGAAGCTGTCACGCACGTCCgg                >  8:145785/1‑149 (MQ=255)
                                                                                                                                        ggggCGCAGATACAATATCTGCGCGGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGcc            >  6:293177/1‑148 (MQ=255)
                                                                                                                                         gggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCa          <  3:167057/149‑1 (MQ=255)
                                                                                                                                         gggCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCa          <  6:90000/149‑1 (MQ=255)
                                                                                                                                           gcgcAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgc        >  2:306137/1‑149 (MQ=255)
                                                                                                                                           gcgcAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgc        >  6:79570/1‑149 (MQ=255)
                                                                                                                                           gcgcAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGAGTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCGTGCTTGAAGCTGTCACGCACGTCCGGCGCCCAgc        >  8:239498/1‑149 (MQ=255)
                                                                                                                                           gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttct                                      <  3:188555/119‑1 (MQ=255)
                                                                                                                                           gcgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCcttct                                      >  4:188555/1‑119 (MQ=255)
                                                                                                                                            cgcAGATATTGTATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcg       >  2:3403/1‑149 (MQ=255)
                                                                                                                                            cgcAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAgcg       >  1:90675/1‑149 (MQ=255)
                                                                                                                                                 aTACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCgg  >  4:319264/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATGCCGCGACGCCCGGCAGACAGCGCGCCAGTCATGGCTTCGGCGGTAATGGCGATAAGGTAGAGCATCAACAACATGGTGCGGGTCCGTGCAGGAAAGGCGCGCAGTCTAACCAGTTGCCCAAGGCACCAAAAGGGGGCGCAGATACAATATCTGCGCCGTTTCGACTCACAGCTTGATGAAGTGCTCGCGATAGTGACGCAGTTCGGCGATCGATTCGCGGATATCGTCCAGCGCCAGGTGGGTATTGCCCTTCTTGAAGCTGTCACGCACGTCCGGCGCCCAGCGCGCGG  >  NC_002947/5572053‑5572346

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: