Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,564,987 A→T 6.8% intergenic (‑62/+19) potF‑II ← / ← PP_3148 putrescine‑binding protein/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,564,9870AT 6.8% 175.6 / 3.1 59intergenic (‑62/+19)potF‑II/PP_3148putrescine‑binding protein/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base T (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

TACACACTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTACTGCTGCACGGTTTCG  >  NC_002947/3564840‑3565135
                                                                                                                                                    |                                                                                                                                                     
tACACACTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCa                                                                                                                                                       >  6:134749/1‑148 (MQ=255)
  cacacTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAg                                                                                                                                                    <  6:180412/149‑1 (MQ=255)
   acacTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGc                                                                                                                                                   <  4:38426/149‑1 (MQ=255)
     acTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGCGCAGGCTTGCCTGTGCATGGGGTGCAAAGCAc                                                                                                                                                 >  4:327154/1‑149 (MQ=255)
     acTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCAc                                                                                                                                                 <  7:60855/149‑1 (MQ=255)
     acTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCAc                                                                                                                                                 >  7:133748/1‑149 (MQ=255)
       tGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCAGGGGGTGCTTTGCAccc                                                                                                                                               >  3:362705/1‑149 (MQ=255)
         aCCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCAccccc                                                                                                                                             >  7:339027/1‑149 (MQ=255)
          ccGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGGGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCTTTGCAccccat                                                                                                                                            <  5:19454/149‑3 (MQ=255)
             tGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCgg                                                                                                                                        >  7:345354/1‑149 (MQ=255)
              gggCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGCGCTTGCAGATGACGCATGGGAACACCTCGTTGT‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGTGGTGCAAAGCACCCCC‑GCGGt                                                                                                                                       >  6:375763/1‑149 (MQ=255)
               ggCGCGGGCTGGGTGGTGGTGGCGTACAGCAGGGTGAATGCAGCGGGGATGAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGTAGCAGGCTTGCCTGTGCATTGGGTGCAAAGCACCCCC‑GCGGTg                                                                                                                                      <  4:356496/149‑1 (MQ=255)
               ggCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTg                                                                                                                                      <  4:152678/149‑1 (MQ=255)
                   cgGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGGATGGGGTGCAAAGCACCCCC‑GCGGTGGATc                                                                                                                                  <  1:31413/149‑1 (MQ=255)
                    ggCCTGGGTGGTGGTGGCGAACAGCAGGGGGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCAGGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCa                                                                                                                                 <  6:220546/149‑1 (MQ=255)
                                 gtggCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACTCTCAGCATc                                                                                                                    <  8:345354/149‑1 (MQ=255)
                                  tggCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCa                                                                                                                   <  5:134749/149‑1 (MQ=255)
                                  tggCGAACACCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGCAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCAGGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCa                                                                                                                   <  8:236511/149‑1 (MQ=255)
                                        aacttccggTGAATGGACGGGGGATCATGGCTTGTACTTGACGCATGGGAACAACCCGTTGGTTCTTG‑TGTTCGTTTTCGGCCCATGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAAtt                                                                                                             <  6:249893/142‑1 (MQ=255)
                                            cagGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACg                                                                                                         >  6:66029/1‑149 (MQ=255)
                                                 tGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCTTTGCAccccat                                                                                                                                            >  3:366678/1‑108 (MQ=255)
                                                 tGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCTTTGCAccccat                                                                                                                                            <  4:366678/110‑3 (MQ=255)
                                                 tGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTccc                                                                                                    <  7:282273/149‑1 (MQ=255)
                                                  gAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTcc                                                                                                     >  8:167581/1‑147 (MQ=255)
                                                      gCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAg                                                                                               >  7:118860/1‑149 (MQ=255)
                                                          cGGTGATCAGGGCTTTCAGATTACGCATGGGAACACCTCGTTGG‑TCTTTTTCTTCGAGTTCGGAGCATGCTTGCCTGTTCATGGGGTGCAAAGCACCCCC‑GCGGTTGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGc                                                                                              <  5:375763/146‑1 (MQ=255)
                                                          cGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCa                                                                                                                                                       <  4:55162/90‑1 (MQ=255)
                                                          cGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCa                                                                                                                                                       >  3:55162/1‑90 (MQ=255)
                                                                tCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGAcgc                                                                                     >  5:155023/1‑149 (MQ=255)
                                                                   gggCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCt                                                                                  >  7:307154/1‑149 (MQ=255)
                                                                       ttGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTa                                                                              <  8:133748/149‑1 (MQ=255)
                                                                        tGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAg                                                                             >  4:125883/1‑149 (MQ=255)
                                                                         gCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCCTTCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTa                                                                              <  1:293470/148‑1 (MQ=255)
                                                                           aGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAg                                                                                                                  >  5:42135/1‑109 (MQ=255)
                                                                           aGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAg                                                                                                                  <  6:42135/109‑1 (MQ=255)
                                                                            gATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTtct                                                                         <  8:118860/149‑1 (MQ=255)
                                                                                     tGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCgg                                                                >  4:225581/1‑149 (MQ=255)
                                                                                       ggAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGgcg                                                              >  4:267230/1‑149 (MQ=255)
                                                                                        gAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGgcgc                                                             <  3:267230/149‑1 (MQ=255)
                                                                                           cacCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCt                                                          <  3:125883/149‑1 (MQ=255)
                                                                                              cTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGa                                                       >  3:71719/1‑149 (MQ=255)
                                                                                                 gTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCg                                                    >  1:310973/1‑149 (MQ=255)
                                                                                                        cTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACg                                               >  8:314359/1‑148 (MQ=255)
                                                                                                          tGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAg                                            >  7:221928/1‑149 (MQ=255)
                                                                                                               ttcGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTg                                       >  6:179264/1‑149 (MQ=255)
                                                                                                                 cGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTgg                                      <  6:155023/148‑1 (MQ=255)
                                                                                                                  gAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGAc                                    <  3:225581/149‑1 (MQ=255)
                                                                                                                     ttCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGaa                                 <  5:66029/149‑1 (MQ=255)
                                                                                                                     gagggAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGaa                                 <  1:320520/146‑1 (MQ=255)
                                                                                                                       cGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAAc                                <  3:327154/148‑1 (MQ=255)
                                                                                                                        ggAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTa                                                                              >  2:358380/1‑101 (MQ=255)
                                                                                                                        ggAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTa                                                                              <  1:358380/101‑1 (MQ=255)
                                                                                                                           gCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCtt                           <  8:307154/149‑1 (MQ=255)
                                                                                                                           gCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCtt                           >  5:17814/1‑149 (MQ=255)
                                                                                                                            cAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTg                          >  2:95012/1‑149 (MQ=255)
                                                                                                                                   gCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAg                     <  7:314359/147‑1 (MQ=255)
                                                                                                                                   gCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAg                     <  5:179264/147‑1 (MQ=255)
                                                                                                                                     cTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTAc                 >  7:21941/1‑149 (MQ=255)
                                                                                                                                        tgCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTActgc              >  6:19432/1‑149 (MQ=255)
                                                                                                                                           aTGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTACtgctgc           <  6:17814/149‑1 (MQ=255)
                                                                                                                                              gggTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTACTGCTGCACg        >  2:202015/1‑149 (MQ=255)
                                                                                                                                                   cAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTACTGCTGCACGGTTTc   <  2:310973/149‑1 (MQ=255)
                                                                                                                                                    aaaGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGAGAGCTTGACCGCCACGTAGCCCTGGACGAACTGCGTTCACAGGTACTGCTGCACgggtacg  >  8:307463/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TACACACTGACCGTGGGCGCGGCCTGGGTGGTGGTGGCGAACAGCAGGGTGAATGCAGCGGGGATCAGGGCTTGCAGATGACGCATGGGAACACCTCGTTGG‑TCTTGTTCTTCGAGTTCGGAGCAGGCTTGCCTGTGCATGGGGTGCAAAGCACCCCC‑GCGGTGGATCAGACACTCAGCATCAGGAATTCACGCTCCCATGAGCTGATGACGCGCTTGTAGTTCTCGTGCTCGGCGCGCTTGACCGCCACGTAGCCCTGGACGAACTGCTTGCCCAGGTACTGCTGCACGGTTTCG  >  NC_002947/3564840‑3565135

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