Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,578,340 T→A 6.6% L161Q (CTG→CAG)  rhdA → thiosulfate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,578,3400TA 6.6% 179.0 / 3.0 61L161Q (CTG→CAG) rhdAthiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/2);  total (30/31)
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.44e-01

CTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCAAA  >  NC_002947/5578195‑5578488
                                                                                                                                                 |                                                                                                                                                    
cTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGt                                                                                                                                                   <  1:218291/149‑1 (MQ=255)
cTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGt                                                                                                                                                   >  5:355229/1‑149 (MQ=255)
     tGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCt                                                                                                                                              <  2:130952/149‑1 (MQ=255)
       gCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTgg                                                                                                                                            <  4:284302/149‑1 (MQ=255)
        cGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTggg                                                                                                                                           <  1:57074/149‑1 (MQ=255)
                               aCAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAAtt                                                                                                                     >  7:179976/1‑148 (MQ=255)
                               aCAAACTCTCCAGCGAAGTGCCTGCCATTGGCAACGCCCCTGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTc                                                                                                                    <  6:210198/149‑1 (MQ=255)
                                cAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAAtt                                                                                                                     >  6:299044/1‑147 (MQ=255)
                                         cACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGaa                                                                                                           >  4:5823/1‑148 (MQ=255)
                                          aCCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGaaa                                                                                                          >  7:188532/1‑148 (MQ=255)
                                                ttgcctgccATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCt                                                                                                   <  7:292308/148‑1 (MQ=255)
                                                  gcctgccATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCAAGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGcg                                                                                                   >  8:292308/1‑146 (MQ=255)
                                                        cATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCa                                                                                           <  3:5823/149‑1 (MQ=255)
                                                             gCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAgcg                                                                                                                <  3:360773/123‑1 (MQ=255)
                                                             gCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAgcg                                                                                                                >  4:360773/1‑123 (MQ=255)
                                                             gCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAggc                                                                                      >  1:335377/1‑149 (MQ=255)
                                                                    cccGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCcaca                                                                               >  3:135104/1‑149 (MQ=255)
                                                                    cccGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCcaca                                                                               >  8:108307/1‑149 (MQ=255)
                                                                        gTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATccc                                                                           <  5:299044/149‑1 (MQ=255)
                                                                            aCCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGgc                                                                       <  2:335377/149‑1 (MQ=255)
                                                                               tGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGgcgc                                                                     >  2:293081/1‑148 (MQ=255)
                                                                                gCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGTCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGa                                                                   >  8:93838/1‑149 (MQ=255)
                                                                                gCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGa                                                                   >  4:82723/1‑149 (MQ=255)
                                                                                  acaTAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATc                                                                 <  7:107633/149‑1 (MQ=255)
                                                                                    aTAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCaa                                                               >  7:66452/1‑149 (MQ=255)
                                                                                     tAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAAc                                                              <  1:82446/149‑1 (MQ=255)
                                                                                      aGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                             >  7:345955/1‑149 (MQ=255)
                                                                                      aGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                             >  1:345776/1‑149 (MQ=255)
                                                                                      aGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACt                                                             >  7:194415/1‑149 (MQ=255)
                                                                                         cAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTc                                                           >  8:319798/1‑148 (MQ=255)
                                                                                             gcACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTg                                                      >  7:72888/1‑149 (MQ=255)
                                                                                              cACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTg                                                      >  2:128988/1‑148 (MQ=255)
                                                                                                 ccACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAAc                                                              >  6:388693/1‑137 (MQ=255)
                                                                                                 ccACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAAc                                                              <  5:388693/137‑1 (MQ=255)
                                                                                                       ccACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACAt                                                                              <  1:361572/115‑1 (MQ=255)
                                                                                                       ccACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACAt                                                                              >  2:361572/1‑115 (MQ=255)
                                                                                                       ccACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGc                                            <  1:128988/149‑1 (MQ=255)
                                                                                                         aCCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCAt                                          <  7:319798/149‑1 (MQ=255)
                                                                                                           ccGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCTTg                                         >  4:255986/1‑148 (MQ=255)
                                                                                                           ccGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATgg                                        <  3:82723/149‑1 (MQ=255)
                                                                                                             gcgAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGAt                                      <  8:179976/149‑1 (MQ=255)
                                                                                                              cgAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGAtc                                     <  3:262378/149‑1 (MQ=255)
                                                                                                                    ccTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAAc                                <  3:255986/148‑1 (MQ=255)
                                                                                                                     cTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGa                              <  1:293081/149‑1 (MQ=255)
                                                                                                                     cTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGa                              >  6:19151/1‑149 (MQ=255)
                                                                                                                         aaaGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAATAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGaccac                          <  2:345776/149‑1 (MQ=255)
                                                                                                                         aaaGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGaccac                          <  8:72888/149‑1 (MQ=255)
                                                                                                                           aGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCt                        <  4:135104/149‑1 (MQ=255)
                                                                                                                            gccgccTGGGCGCCGATGACCTGGTCATCTGCGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTg                       <  8:194415/149‑1 (MQ=255)
                                                                                                                                      cccagATGACCAGGTCATCGGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTc                >  7:144655/5‑146 (MQ=255)
                                                                                                                                      cccagATGACCAGGTCATCGGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTc                <  8:144655/142‑1 (MQ=255)
                                                                                                                                       gcCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCa            <  6:355229/149‑1 (MQ=255)
                                                                                                                                       gcCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCa            <  8:345955/149‑1 (MQ=255)
                                                                                                                                        ccGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAg           >  8:24492/1‑149 (MQ=255)
                                                                                                                                        cagATGACCAGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAggcggc                                                                                   <  5:325771/75‑1 (MQ=255)
                                                                                                                                        cagATGACCAGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAggcggc                                                                                   >  6:325771/3‑77 (MQ=255)
                                                                                                                                              aCCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACAt      <  7:93838/148‑1 (MQ=255)
                                                                                                                                               ccTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGCTCTCAACGACCACCTGCGCATTCGCCAGGACATc     >  4:329672/1‑148 (MQ=255)
                                                                                                                                                cTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAgg          >  5:388747/1‑142 (MQ=255)
                                                                                                                                                cTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAgg          <  6:388747/142‑1 (MQ=255)
                                                                                                                                                cTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCa    >  1:180837/1‑148 (MQ=255)
                                                                                                                                                 tGGTCATCTGGGACGCACGCGGCCCGCAAGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCaaa  >  3:94998/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTCAATGGCGGTATCCAGGCCTGGCCGGCAGACAAACTGTCCACCGAAGTGCCTGCCATTGGCAACGCCCCGGTGCACCTGCACATAGACAGCGCACCCACTGCCACCCGCGAGTACCTGCAAAGCCGCCTGGGCGCCGATGACCTGGTCATCTGGGACGCACGCGGCCCGCAGGAATTCAGCGGCGAAAAAGTGCTGGCAGCCAAAGGCGGCCACATCCCCGGCGCGATCAACTTCGAGTGGACCGCCGGCATGGATCTCAACGACCACCTGCGCATTCGCCAGGACATCAAA  >  NC_002947/5578195‑5578488

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: ‑