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)
         cAGCAAAATCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCggg                                                                                                                                         <  2:799009/150‑1 (MQ=255)
           gCAAAAACAACGATTTGGAGGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCTCTGGACGCGATGTAAATCAGTCGCGCACTTTTTCCGC‑TTTTCCCGGGCtt                                                                                                                                      <  2:335780/151‑1 (MQ=255)
             aaaaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTtct                                                                                                                                    <  1:458450/151‑1 (MQ=255)
              aaaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTtctc                                                                                                                                   <  1:659436/151‑1 (MQ=255)
               aaTCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTctcg                                                                                                                                  <  1:623438/151‑1 (MQ=255)
                 tCAACGATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTctcgct                                                                                                                                >  1:108151/1‑151 (MQ=255)
                      gATTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTctcgctcg                                                                                                                              >  2:407989/1‑148 (MQ=255)
                       aTTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCgagag                                                                                                                          <  2:188190/151‑1 (MQ=255)
                       aTTTGGATGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCgag                                                                                                                            >  2:544245/1‑149 (MQ=255)
                             aTGCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGaaa                                                                                                                     >  1:180655/1‑150 (MQ=255)
                               gCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAaat                                                                                                                   <  2:329906/150‑1 (MQ=255)
                               gCTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAaat                                                                                                                   <  1:326930/150‑1 (MQ=255)
                                cTGACGAGTTCCTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAaata                                                                                                                  >  2:382931/1‑150 (MQ=255)
                                          ccTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTa                                                                                                        >  1:661344/1‑150 (MQ=255)
                                          ccTCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTa                                                                                                        <  2:710008/150‑1 (MQ=255)
                                            tCGAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAAtt                                                                                                     >  2:44611/1‑151 (MQ=255)
                                              gAACACGTAGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGc                                                                                                   >  1:447954/1‑151 (MQ=255)
                                                      aGCGCGAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCAc                                                                                           <  2:74215/151‑1 (MQ=255)
                                                        cgcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCg                                                                                         <  2:163309/151‑1 (MQ=255)
                                                         gcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGt                                                                                        >  2:88971/1‑151 (MQ=255)
                                                         gcgAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGt                                                                                        <  1:607752/151‑1 (MQ=255)
                                                           gAAATCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGATGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCg                                                                                      >  1:795362/1‑151 (MQ=255)
                                                             aaTCACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCggg                                                                                    <  2:493564/151‑1 (MQ=255)
                                                                cACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTg                                                                                 <  2:180631/151‑1 (MQ=255)
                                                                cACCCTGCGCCGCAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTg                                                                                 <  2:224880/151‑1 (MQ=255)
                                                                 acacTGCGCCGCAGGCTCCTCGATATATCTACAAGCTTGAGCAGGCAGTGGACGCGATGTAAAACAGTCACGCACTTTGGCCGC‑TTTTCCCGGGCTTCTCGCGCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCt                                                                                  <  2:737645/146‑1 (MQ=255)
                                                                       cgccgcAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTaa                                                                          <  2:897216/151‑1 (MQ=255)
                                                                           gcAGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGaaa                                                                      <  1:714992/151‑1 (MQ=255)
                                                                             aGGCTCCGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAAc                                                                    <  1:242648/151‑1 (MQ=255)
                                                                                  ccGCGCTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGaa                                                                <  2:400459/150‑1 (MQ=255)
                                                                                     cgcTATATCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTc                                                            >  1:529716/1‑151 (MQ=255)
                                                                                           atCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACgg                                                      >  2:264374/1‑151 (MQ=255)
                                                                                           atCTACAAACTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACg                                                       <  2:457093/150‑1 (MQ=255)
                                                                                                 aaaCTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtg                                                >  2:39750/1‑151 (MQ=255)
                                                                                                  aaCTTGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtgg                                               <  1:202460/151‑1 (MQ=255)
                                                                                                      tGAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAgtggtg                                             <  1:793604/149‑1 (MQ=255)
                                                                                                       gAGCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGAt                                           <  1:652208/150‑1 (MQ=255)
                                                                                                         gCAGGCACTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGt                                        <  1:428981/151‑1 (MQ=255)
                                                                                                           aGGCACTGGACGCGATGTAAATCAGCCTCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTTGCTCGAGAGGGAAAATAATCCGCTGTAATTGCCGCTCCACCTTCGGGCTGACGTTAAGAAAACGCAAGCTCAGATGGGGAGCGGTGATGGt                                        >  2:899031/1‑149 (MQ=255)
                                                                                                               aCTGGACGCGATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCAtt                                  <  1:825256/151‑1 (MQ=255)
                                                                                                                      gcgATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGcc                            >  2:575644/1‑150 (MQ=255)
                                                                                                                      gcgATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCa                           <  1:49332/151‑1 (MQ=255)
                                                                                                                        gATGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatc                         >  1:492242/1‑151 (MQ=255)
                                                                                                                         aTGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatca                        <  2:481973/151‑1 (MQ=255)
                                                                                                                         aTGTAAATCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatc                         >  2:826141/1‑150 (MQ=255)
                                                                                                                             aaaTCAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCatcaatc                     >  1:858703/1‑150 (MQ=255)
                                                                                                                                 cAGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACtt                 <  1:74058/150‑1 (MQ=255)
                                                                                                                                  aGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACttt                >  2:909194/1‑150 (MQ=255)
                                                                                                                                  aGTCGCGCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTg               <  1:551166/151‑1 (MQ=255)
                                                                                                                                     cgcgCACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTgcg             <  2:275878/150‑1 (MQ=255)
                                                                                                                                      gcgcACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgct           >  1:672557/1‑151 (MQ=255)
                                                                                                                                      gcgcACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgct           >  2:462860/1‑151 (MQ=255)
                                                                                                                                       cgcACTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgctc          >  2:421962/1‑151 (MQ=255)
                                                                                                                                          aCTTTGTCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGTGGTGATGGTTTCATTCTTGCCATCAATCACTTTGcgctcgct       <  2:795297/151‑1 (MQ=255)
                                                                                                                                                tCCGC‑TTTTCCCGGGCTTCTCGCTCGAGAGAGAAAATAATCCGCTGTAATTGCCGCTCCACCGTCGGGCTGACGTTAAGAAAACGGAAGCTCAGACGGGGAGCGGTGATGGTTTCATTCTTGCCATCAATCACTTTGCGCTCGCTGATgg  >  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: ‑