Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,243,818 (T)5→4 coding (710/735 nt) ycgR ← flagellar velocity braking protein, c‑di‑GMP‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,243,8140T.98.2% 258.9 / ‑1.7 56coding (714/735 nt)ycgRflagellar velocity braking protein, c‑di‑GMP‑regulated
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (23/32);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

AAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGG  >  NC_000913/1243665‑1243959
                                                                                                                                                     |                                                                                                                                                 
aaaGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCtt                                                                                                                                                  >  1:668494/1‑151 (MQ=255)
 aaGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCttt                                                                                                                                                 <  2:291540/151‑1 (MQ=255)
  aGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGTTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCtttt                                                                                                                                                <  2:481983/151‑1 (MQ=255)
  aGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCtttt                                                                                                                                                >  1:26099/1‑151 (MQ=255)
   ggCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCttttc                                                                                                                                               <  1:687685/151‑2 (MQ=255)
     cGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTccc                                                                                                                                             >  2:74036/1‑151 (MQ=255)
         cAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCggg                                                                                                                                         <  2:799009/150‑1 (MQ=255)
           gCAAAAACAACGATTTGGAGGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCTCTGGACGCGATGTAAATCAGTCGCGCACTTTTTCCGCTTTTCCCGGGCtt                                                                                                                                      <  2:335780/151‑1 (MQ=255)
             aaaaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTtct                                                                                                                                    <  1:458450/151‑1 (MQ=255)
              aaaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTtctc                                                                                                                                   <  1:659436/151‑1 (MQ=255)
               aaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTctcg                                                                                                                                  <  1:623438/151‑1 (MQ=255)
                 tCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTctcgct                                                                                                                                >  1:108151/1‑151 (MQ=255)
                      gATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTctcgctcg                                                                                                                              >  2:407989/1‑148 (MQ=255)
                       aTTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCgagag                                                                                                                          <  2:188190/151‑1 (MQ=255)
                       aTTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCgag                                                                                                                            >  2:544245/1‑149 (MQ=255)
                             aTGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGaaa                                                                                                                     >  1:180655/1‑150 (MQ=255)
                               gCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAaat                                                                                                                   <  2:329906/150‑1 (MQ=255)
                               gCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAaat                                                                                                                   <  1:326930/150‑1 (MQ=255)
                                cTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAaata                                                                                                                  >  2:382931/1‑150 (MQ=255)
                                          ccTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTa                                                                                                        >  1:661344/1‑150 (MQ=255)
                                          ccTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTa                                                                                                        <  2:710008/150‑1 (MQ=255)
                                            tCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAAtt                                                                                                     >  2:44611/1‑151 (MQ=255)
                                              gAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGc                                                                                                   >  1:447954/1‑151 (MQ=255)
                                                      aGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCAc                                                                                           <  2:74215/151‑1 (MQ=255)
                                                        cgcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCg                                                                                         <  2:163309/151‑1 (MQ=255)
                                                         gcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGt                                                                                        >  2:88971/1‑151 (MQ=255)
                                                         gcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGt                                                                                        <  1:607752/151‑1 (MQ=255)
                                                           gAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGATGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCg                                                                                      >  1:795362/1‑151 (MQ=255)
                                                             aaTCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCggg                                                                                    <  2:493564/151‑1 (MQ=255)
                                                                cACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTg                                                                                 <  2:180631/151‑1 (MQ=255)
                                                                cACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTg                                                                                 <  2:224880/151‑1 (MQ=255)
                                                                 acacTGCGCCGCAGGCTCCTCGATATATCTACAAGCTTGAGCAGGCAGTGGACGCGATGTAAAACAGTCACGCACTTTGGCCGCTTTTCCCGGGCTTCTCGCGCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCt                                                                                  <  2:737645/146‑1 (MQ=255)
                                                                       cgccgcAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTaa                                                                          <  2:897216/151‑1 (MQ=255)
                                                                           gcAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGaaa                                                                      <  1:714992/151‑1 (MQ=255)
                                                                             aGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAAc                                                                    <  1:242648/151‑1 (MQ=255)
                                                                                  ccGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGaa                                                                <  2:400459/150‑1 (MQ=255)
                                                                                     cgcTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTc                                                            >  1:529716/1‑151 (MQ=255)
                                                                                           atCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACgg                                                      >  2:264374/1‑151 (MQ=255)
                                                                                           atCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACg                                                       <  2:457093/150‑1 (MQ=255)
                                                                                                 aaaCTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtg                                                >  2:39750/1‑151 (MQ=255)
                                                                                                  aaCTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtgg                                               <  1:202460/151‑1 (MQ=255)
                                                                                                      tGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtggtg                                             <  1:793604/149‑1 (MQ=255)
                                                                                                       gAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGAt                                           <  1:652208/150‑1 (MQ=255)
                                                                                                         gCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGt                                        <  1:428981/151‑1 (MQ=255)
                                                                                                           aGGCACTGGACGCGATGTAAATCAGCCTCGCACTTTGTCCGCTTTTCCCGGGCTTCTTGCTCGAGAGGGAAAATAATCCGCTGTAATTGCCGCTCCACCTTCGGGCTGACGTTAAGAAAACGCAAGCTCAGATGGGGAGCGGTGATGGt                                        >  2:899031/1‑149 (MQ=255)
                                                                                                               aCTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCAtt                                  <  1:825256/151‑1 (MQ=255)
                                                                                                                      gcgATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGcc                            >  2:575644/1‑150 (MQ=255)
                                                                                                                      gcgATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCa                           <  1:49332/151‑1 (MQ=255)
                                                                                                                        gATGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatc                         >  1:492242/1‑151 (MQ=255)
                                                                                                                         aTGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatca                        <  2:481973/151‑1 (MQ=255)
                                                                                                                         aTGTAAATCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatc                         >  2:826141/1‑150 (MQ=255)
                                                                                                                             aaaTCAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatcaatc                     >  1:858703/1‑150 (MQ=255)
                                                                                                                                 cAGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACtt                 <  1:74058/150‑1 (MQ=255)
                                                                                                                                  aGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACttt                >  2:909194/1‑150 (MQ=255)
                                                                                                                                  aGTCGCGCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTg               <  1:551166/151‑1 (MQ=255)
                                                                                                                                     cgcgCACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTgcg             <  2:275878/150‑1 (MQ=255)
                                                                                                                                      gcgcACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgct           >  1:672557/1‑151 (MQ=255)
                                                                                                                                      gcgcACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgct           >  2:462860/1‑151 (MQ=255)
                                                                                                                                       cgcACTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgctc          >  2:421962/1‑151 (MQ=255)
                                                                                                                                          aCTTTGTCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgctcgct       <  2:795297/151‑1 (MQ=255)
                                                                                                                                                tCCGCTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGCGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATgg  >  2:162151/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AAAGGCGATCAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGCTTTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATGG  >  NC_000913/1243665‑1243959

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