Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,174,134 A→T 6.6% L210Q (CTG→CAG)  PP_1927 ← arsenical resistance protein ArsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,174,1340AT 6.6% 184.3 / 2.7 60L210Q (CTG→CAG) PP_1927arsenical resistance protein ArsH
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base T (2/2);  total (30/30)
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.42e-01

GGGCAAAGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAAACATACG  >  NC_002947/2173994‑2174270
                                                                                                                                            |                                                                                                                                        
gggCAAAGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGt                                                                                                                                   <  3:205828/148‑1 (MQ=255)
    aaaGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCagcagcagc                                                                                                                                     >  4:231346/1‑142 (MQ=255)
    aaaGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCagcagcagc                                                                                                                                     <  3:231346/142‑1 (MQ=255)
    aaaGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCAGCCGACTCTTTGCGTTCGGAGTAGCTGTCTACCAAGTAATCCTGGCGATCACGCTGCTGCTGCGTGatc                                                                                                                               <  6:160646/148‑3 (MQ=255)
      aGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTg                                                                                                                            <  6:112850/149‑1 (MQ=255)
      aGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATACTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTg                                                                                                                            <  7:130013/149‑1 (MQ=255)
                  gTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTtcctc                                                                                                                 >  3:363850/1‑148 (MQ=255)
                                  cTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGcc                                                                                                >  5:45990/1‑149 (MQ=255)
                                   tCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGcc                                                                                                >  3:194986/1‑148 (MQ=255)
                                   tCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGcc                                                                                                >  4:227620/1‑148 (MQ=255)
                                     aTTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGAt                                                                                             >  5:168392/1‑149 (MQ=255)
                                      ttATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATc                                                                                            <  6:45990/149‑1 (MQ=255)
                                           tcAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTg                                                                                                                            >  8:292269/1‑112 (MQ=255)
                                           tcAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTg                                                                                                                            <  7:292269/112‑1 (MQ=255)
                                           tcAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAg                                                                                        >  7:114767/1‑148 (MQ=255)
                                           tcAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGAAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAg                                                                                        >  3:199177/1‑148 (MQ=255)
                                             aaaTAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAg                                                                                        <  4:175485/146‑1 (MQ=255)
                                              aaTAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGa                                                                                       >  7:165558/1‑146 (MQ=255)
                                               aTAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATc                                                                                                             <  2:83539/123‑1 (MQ=255)
                                               aTAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATc                                                                                                             >  1:83539/1‑123 (MQ=255)
                                                    ccGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGacgacg                                                                              >  1:131656/1‑149 (MQ=255)
                                                      gAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGt                                                                            >  4:66850/1‑149 (MQ=255)
                                                         tatTGGCTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTc                                                                         <  3:294586/149‑1 (MQ=255)
                                                          atTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCa                                                                        <  8:193319/149‑1 (MQ=255)
                                                            tGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCAt                                                                       <  8:165558/148‑1 (MQ=255)
                                                                tctctTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTc                                                                  >  5:152738/1‑149 (MQ=255)
                                                                    ttGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGc                                                              >  4:294556/1‑149 (MQ=255)
                                                                      ggATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGATGGTTTCATGCGTCCGGc                                                              <  4:199177/147‑1 (MQ=255)
                                                                      ggATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCgg                                                               <  2:131656/146‑1 (MQ=255)
                                                                         tAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCAt                                                         >  2:75174/1‑149 (MQ=255)
                                                                           gTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATc                                                        <  3:176767/148‑1 (MQ=255)
                                                                           gTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATc                                                        <  1:75174/148‑1 (MQ=255)
                                                                              cttcCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAAc                                                    >  3:145501/1‑149 (MQ=255)
                                                                              cttcCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAAc                                                    >  3:139915/1‑149 (MQ=255)
                                                                               ttcCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGaa                                                     <  6:317010/147‑1 (MQ=255)
                                                                               ttcCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGaa                                                     >  5:317010/1‑147 (MQ=255)
                                                                                    ccgACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAgg                                              <  6:168392/149‑1 (MQ=255)
                                                                                      gACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAgg                                              >  4:160510/1‑147 (MQ=255)
                                                                                            ttGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCttt                                      >  2:49140/1‑149 (MQ=255)
                                                                                              gCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCtttt                                     <  4:363850/148‑1 (MQ=255)
                                                                                               cGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTg                                    <  3:227620/148‑1 (MQ=255)
                                                                                               cGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTg                                    <  4:194986/148‑1 (MQ=255)
                                                                                                     gAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCa                                                                        >  8:283628/1‑106 (MQ=255)
                                                                                                     gAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCa                                                                        <  7:283628/106‑1 (MQ=255)
                                                                                                        tAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACcgacga                          <  6:152738/149‑1 (MQ=255)
                                                                                                                tACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCg                  >  3:326870/1‑149 (MQ=255)
                                                                                                                  ccAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTc                   <  8:212979/146‑1 (MQ=255)
                                                                                                                  ccAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCg                  <  4:145501/147‑1 (MQ=255)
                                                                                                                  ccAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCg                  <  4:139915/147‑1 (MQ=255)
                                                                                                                  ccAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGACGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTc                   >  7:212979/1‑146 (MQ=255)
                                                                                                                     aGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATg             <  1:153110/149‑1 (MQ=255)
                                                                                                                     aGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATg             >  7:111150/1‑149 (MQ=255)
                                                                                                                     aGTAATCCTCGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATg             <  7:70995/149‑1 (MQ=255)
                                                                                                                       tAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATgg            >  5:21932/1‑148 (MQ=255)
                                                                                                                         aTCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTaa         <  8:114767/149‑1 (MQ=255)
                                                                                                                         aTCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTa          >  6:327255/1‑148 (MQ=255)
                                                                                                                          tCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCACTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTaaa        <  2:231578/149‑1 (MQ=255)
                                                                                                                           ccTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTaaa        >  5:213201/1‑148 (MQ=255)
                                                                                                                                agttCACGCTGCTGCTGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAAACATAc   >  7:14286/4‑148 (MQ=255)
                                                                                                                                agttCACGCTGCTGCTGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAAACATAc   <  8:14286/145‑1 (MQ=255)
                                                                                                                                 gttCACGCTGCTGCTGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAATCATACg  >  4:266059/3‑148 (MQ=255)
                                                                                                                                 gATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGCCGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAAACATACg  >  2:124365/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                        
GGGCAAAGATTTCAATGAGTAAGGATGCGACATGCTCATTATCTCAAATAGACCGAATATTGACTCTCTTGGATAGTTCTTCCGCCGACTCTTTGCGTTCGGAGTAGCGGTCTACCAAGTAATCCTGGCGATCACGCAGCAGCAGCGTGAACTTGGTCAGTTCCTCCATCACGTCGACCAGCCGATCGTAGAAAGACGACGGTTTCATGCGTCCGGCTTCATCGAACTCCAGGAAGGCTTTTGGCACCGACGACTGGTTCGGGATGGTAAACATACG  >  NC_002947/2173994‑2174270

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