Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,058,717 G→C 6.2% V57V (GTC→GTG)  sodB ← superoxide dismutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,058,7170GC 6.2% 259.1 / 5.2 79V57V (GTC→GTG) sodBsuperoxide dismutase
Reads supporting (aligned to +/- strand):  ref base G (31/43);  new base C (3/2);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATGGGCGTACGG  >  NC_002947/1058570‑1058863
                                                                                                                                                    |                                                                                                                                                   
cTTGTCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGa                                                                                                                                                    <  6:165970/149‑1 (MQ=255)
cTTGTCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGa                                                                                                                                                    <  6:163783/149‑1 (MQ=255)
   gTCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGa                                                                                                                                                 >  7:323750/1‑149 (MQ=255)
    tCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGAt                                                                                                                                                <  4:235626/149‑1 (MQ=255)
      gAAGGAACCGAAAGCGGCGTTGATGGC‑TTCAGCCAGTGAAACGGTAGGCTGGCCGCCGCCGTTTGGAGACAGGCAGTTACAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGACCGAAGAGCTCTTGAAGAtct                                                                                                                                              >  5:324226/1‑149 (MQ=255)
      gAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGATCTCTTGACGAtct                                                                                                                                              >  4:90864/1‑149 (MQ=255)
       aaGGAACCGAAAGCGGCGTTGATGGCGGTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGAtct                                                                                                                                              >  1:2846/1‑149 (MQ=255)
       aaGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGTCGCCCGAAGAGCTCTTGACGAtctc                                                                                                                                             >  6:41805/1‑149 (MQ=255)
       aaGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGAtctc                                                                                                                                             >  4:25280/1‑149 (MQ=255)
       aaGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGAtctc                                                                                                                                             >  3:185168/1‑149 (MQ=255)
            aCCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCa                                                                                                                                        <  6:140397/149‑1 (MQ=255)
              cGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAgg                                                                                                                                      <  2:2846/149‑1 (MQ=255)
              cGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAgg                                                                                                                                      <  3:29302/149‑1 (MQ=255)
              cGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAgg                                                                                                                                      <  3:24205/149‑1 (MQ=255)
              cGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAgg                                                                                                                                      <  1:224647/149‑1 (MQ=255)
              cGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAgg                                                                                                                                      <  8:77906/149‑1 (MQ=255)
                 aaGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTc                                                                                                                                   <  6:251518/149‑1 (MQ=255)
                         ttGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTc                                                                                                                           <  4:89853/149‑1 (MQ=255)
                         ttGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTc                                                                                                                           <  8:323750/149‑1 (MQ=255)
                            aTGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGaa                                                                                                                        <  6:196787/149‑1 (MQ=255)
                            aTGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGaa                                                                                                                        >  1:222076/1‑149 (MQ=255)
                                c‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCg                                                                                                                    <  2:222076/149‑1 (MQ=255)
                                  gTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCgg                                                                                                                   <  1:193446/149‑1 (MQ=255)
                                   tCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGTCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGt                                                                                                                  <  5:41805/149‑1 (MQ=255)
                                   tCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGt                                                                                                                  <  4:185168/149‑1 (MQ=255)
                                      gCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGcc                                                                                                               >  6:191391/1‑149 (MQ=255)
                                       ccAGTCCACCGGTAGGCTGGCCGCCGCAGTTAGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTAGAGCGGCGTTGGTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGAGCCt                                                                                                              >  8:262857/1‑149 (MQ=255)
                                               ccGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGGAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAg                                                                                                      >  6:110352/1‑149 (MQ=255)
                                                      gCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGgttgtt                                                                                                <  3:90864/148‑1 (MQ=255)
                                                      gCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTc                                                                                               <  2:215440/149‑1 (MQ=255)
                                                          gccgccgccGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGGAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAgg                                                                                            <  5:110352/148‑1 (MQ=255)
                                                          gccgccgccGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAgg                                                                                            <  3:134436/148‑1 (MQ=255)
                                                          gccgccgccGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGt                                                                                           >  8:5283/1‑149 (MQ=255)
                                                          gccgccaccGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGt                                                                                           >  3:32501/1‑149 (MQ=255)
                                                              ccgccgTTTGGCGACAGGCAGTTCCAGTATAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGTAGAGCTGTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCAc                                                                                       <  3:31784/149‑1 (MQ=255)
                                                                gccgTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGa                                                                                     <  3:75297/149‑1 (MQ=255)
                                                                     ttGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAg                                                                                >  6:126593/1‑149 (MQ=255)
                                                                        gCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGgtg                                                                             <  5:126593/149‑1 (MQ=255)
                                                                             aGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTtg‑tg                                                                       <  7:5283/149‑1 (MQ=255)
                                                                             aGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTtg‑tg                                                                       >  1:200456/1‑149 (MQ=255)
                                                                                        gTAGAAGGTGTGGTTCCAGAGTTGTGCGGGGTTGTGGAAGTTGCCCCCAGATGCGCTCTTCACGCTCTCTTCCAGGGTCTTGCCTTCTAATTCGGTGCCTGGTACCAGGTTGTTCAGTTTCACGACATAGGTGTTG‑TGGTGCTTGTCgt                                                            <  2:84509/149‑1 (MQ=255)
                                                                                        gTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgt                                                            >  4:2709/1‑149 (MQ=255)
                                                                                        gTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgt                                                            >  3:109957/1‑149 (MQ=255)
                                                                                         tAGAATGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtg                                                           >  4:131486/1‑149 (MQ=255)
                                                                                           gAAGGTGTGGTTCCAGACTTGTGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtgg                                                          <  5:191391/148‑1 (MQ=255)
                                                                                           gAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtggt                                                         <  3:277185/149‑1 (MQ=255)
                                                                                           gAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtgg                                                          <  7:180255/148‑1 (MQ=255)
                                                                                           gAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtgg                                                          <  2:200456/148‑1 (MQ=255)
                                                                                               gtgtGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATa                                                     >  2:80575/1‑149 (MQ=255)
                                                                                                    gTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTc                                                  >  4:314786/1‑147 (MQ=255)
                                                                                                    gTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTc                                                  >  7:119950/1‑147 (MQ=255)
                                                                                                    gTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCa                                                <  3:25280/149‑1 (MQ=255)
                                                                                                           aCTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTc                                         <  6:154478/149‑1 (MQ=255)
                                                                                                             ttGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTg‑tggt                                                                     >  8:273592/1‑119 (MQ=255)
                                                                                                             ttGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTg‑tggt                                                                     <  7:273592/119‑1 (MQ=255)
                                                                                                             ttGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTGCTGGTGCTTGTCGTGGTGATACTCCAg                                               <  1:4681/142‑1 (MQ=255)
                                                                                                             ttGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAg                                               >  2:4681/1‑141 (MQ=255)
                                                                                                                 gcggcgTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGga                                   <  8:119950/149‑1 (MQ=255)
                                                                                                                 gcggcgTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGga                                   <  3:314786/149‑1 (MQ=255)
                                                                                                                   ggcgTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTg                                     >  8:300810/1‑145 (MQ=255)
                                                                                                                            gAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCt                        >  2:168826/1‑149 (MQ=255)
                                                                                                                            gAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCt                        >  5:50918/1‑149 (MQ=255)
                                                                                                                             aaGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtg                                                           <  8:110213/113‑1 (MQ=255)
                                                                                                                             aaGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCgtg                                                           >  7:110213/1‑113 (MQ=255)
                                                                                                                             aaGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTg                       >  6:91317/1‑149 (MQ=255)
                                                                                                                              aGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGc                      <  1:168826/149‑1 (MQ=255)
                                                                                                                               gATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTTGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCa                     <  4:32501/149‑1 (MQ=255)
                                                                                                                                  gccgccCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAggg                  >  6:323076/1‑149 (MQ=255)
                                                                                                                                   ccgccCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGc                 <  5:91317/149‑1 (MQ=255)
                                                                                                                                    cgccCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGca                <  3:218508/149‑1 (MQ=255)
                                                                                                                                    ccctggAAGAGATCGTCAAGAGCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGca                >  7:225867/6‑149 (MQ=255)
                                                                                                                                     cctggAAGAGATCGTCAAGAGCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACAACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGca                >  3:265631/5‑148 (MQ=255)
                                                                                                                                      cccGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTg                                     <  7:206377/126‑1 (MQ=255)
                                                                                                                                      cccGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTg                                     >  8:206377/1‑126 (MQ=255)
                                                                                                                                       tggAAGAGATCGTCAAGAGCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGcatca             >  3:203974/3‑149 (MQ=255)
                                                                                                                                          aaGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATgg          <  3:2709/149‑1 (MQ=255)
                                                                                                                                          aaGAGATCGTCAAGAGCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACAACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATgg          <  4:265631/149‑1 (MQ=255)
                                                                                                                                           agagCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATgg          <  4:109957/148‑1 (MQ=255)
                                                                                                                                            gagCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGc                                                                   >  4:124694/1‑90 (MQ=255)
                                                                                                                                                    gACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATGGGCGTACgg  >  3:337942/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTGTCGAAGGAACCGAAAGCGGCGTTGATGGC‑GTCAGCCAGTGCACCGGTAGGCTGGCCGCCGCCGTTTGGCGACAGGCAGTTCCAGTAGAAGGTGTGGTTCCAGACTTGAGCGGCGTTGTTGAAGATGCCGCCCGAAGAGCTCTTGACGATCTCTTCCAGGGTCTTGCCTTCGAATTCGGTGCCTGGGACCAGGTTGTTCAGGTTCACGACATAGGTGTTG‑TGGTGCTTGTCGTGGTGATACTCCAGGGTTTCCTTGGAGATGTGCGGCTGCAGGGCATCATGGGCGTACGG  >  NC_002947/1058570‑1058863

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