Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,905,851 T→C 6.8% S58G (AGC→GGC)  PP_4314 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,905,8510TC 6.8% 166.2 / 2.7 59S58G (AGC→GGC) PP_4314hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (2/2);  total (29/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.40e-01

TAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAGGG  >  NC_002947/4905703‑4905994
                                                                                                                                                    |                                                                                                                                                
tagtaCAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgct                                                                                                                                                  <  5:181317/149‑1 (MQ=255)
  gtaCAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgctgc                                                                                                                                                <  3:111898/149‑1 (MQ=255)
     cAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaa                                                                                                                                              <  6:325729/148‑1 (MQ=255)
       ggTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGa                                                                                                                                           >  7:257656/1‑149 (MQ=255)
        gTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGAc                                                                                                                                          <  8:257656/149‑1 (MQ=255)
         tGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaagt                                                                                                                                            <  8:128850/146‑3 (MQ=255)
         tGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCAGCATGTGCCGCaagt                                                                                                                                            >  7:128850/1‑144 (MQ=255)
               ccACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCt                                                                                                                                   >  4:122942/1‑149 (MQ=255)
                   aGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGata                                                                                                                               <  1:112056/149‑1 (MQ=255)
                    gCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGTACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGatat                                                                                                                              <  1:255946/149‑1 (MQ=255)
                      ggCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATgg                                                                                                                            <  8:63593/149‑1 (MQ=255)
                      ggCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATgg                                                                                                                            <  3:179288/149‑1 (MQ=255)
                       gCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATggg                                                                                                                           <  1:196086/149‑1 (MQ=255)
                           caGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACAtgctg                                                                                                                                                 <  2:7956/123‑1 (MQ=255)
                            aGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTg                                                                                                                      >  5:118743/1‑149 (MQ=255)
                                  cACGCTGGGCCCAGTCGTATGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACg                                                                                                                <  7:288395/149‑1 (MQ=255)
                                   aCGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACgg                                                                                                               >  3:267859/1‑149 (MQ=255)
                                          gCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGt                                                                                                        >  5:132855/1‑149 (MQ=255)
                                            ccAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCa                                                                                                      <  1:154150/149‑1 (MQ=255)
                                                  ggACGAATCATAGACGCTGGGTGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgcg                                                                                                >  2:119723/1‑149 (MQ=255)
                                                  ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgcg                                                                                                >  7:66690/1‑149 (MQ=255)
                                                  ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgc                                                                                                 <  3:122942/148‑1 (MQ=255)
                                                  ggACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCcgc                                                                                                 <  6:291088/148‑1 (MQ=255)
                                                   gACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               <  6:118743/149‑1 (MQ=255)
                                                      gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               >  5:174298/1‑146 (MQ=255)
                                                      gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCgt                                                                                               <  6:174298/146‑1 (MQ=255)
                                                      gAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTc                                                                                            >  3:119755/1‑149 (MQ=255)
                                                            tAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  7:268320/1‑149 (MQ=255)
                                                             aGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  1:303773/1‑148 (MQ=255)
                                                              gACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGca                                                                                      >  7:83109/1‑147 (MQ=255)
                                                                 gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcaggc                                                                                 <  4:267859/149‑1 (MQ=255)
                                                                 gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcagg                                                                                  >  4:165254/1‑148 (MQ=255)
                                                                 gCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcagg                                                                                  >  1:153684/1‑148 (MQ=255)
                                                                    gggCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCAcaggcagg                                                                              >  4:70349/1‑149 (MQ=255)
                                                                        gcACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGcc                                                                          >  7:74388/1‑149 (MQ=255)
                                                                              ccTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGttt                                                                    <  3:78345/149‑1 (MQ=255)
                                                                               cTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGtttt                                                                   <  3:165254/149‑1 (MQ=255)
                                                                                   gATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATTCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCa                                                               <  2:249197/149‑1 (MQ=255)
                                                                                     tCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAgg                                                             >  5:327720/1‑149 (MQ=255)
                                                                                      cAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAgg‑g                                                           >  2:164315/1‑149 (MQ=255)
                                                                                       aGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑Gc                                                          <  1:119723/149‑1 (MQ=255)
                                                                                        gCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑Gc                                                          >  6:234304/1‑148 (MQ=255)
                                                                                            aGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATgg                                                     <  8:268320/149‑1 (MQ=255)
                                                                                              gCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGt                                                   <  2:153684/149‑1 (MQ=255)
                                                                                                    tgcAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTcc                                             >  4:135013/1‑149 (MQ=255)
                                                                                                    tgcAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTcc                                             >  2:202508/1‑149 (MQ=255)
                                                                                                       aGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCtt                                           <  6:327720/148‑1 (MQ=255)
                                                                                                       aGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTg                                          <  8:175216/149‑1 (MQ=255)
                                                                                                              caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  4:270251/1‑148 (MQ=255)
                                                                                                              caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  3:21233/1‑148 (MQ=255)
                                                                                                              caccaGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAAAGGa                                    >  2:252471/1‑148 (MQ=255)
                                                                                                                 cagcagCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGa                                    >  1:314115/1‑145 (MQ=255)
                                                                                                                   gcagcGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATg                              <  8:74388/149‑1 (MQ=255)
                                                                                                                     agcGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATggg                            <  4:21233/149‑1 (MQ=255)
                                                                                                                           ccgccgGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCAAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGc                       >  4:5221/1‑148 (MQ=255)
                                                                                                                              ccgGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGAATTGTCGACGGTGACCGTCAGCCGCTTGTCGCCGCACAGGCAGGTGCCTTGTTTTTCCAGC‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAAtt                   <  2:283436/149‑1 (MQ=255)
                                                                                                                                 ttgtcattGCAGCATGTGCCGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATta                  >  5:265507/7‑147 (MQ=255)
                                                                                                                                 gTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATtata                <  3:270251/149‑1 (MQ=255)
                                                                                                                                   ccACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAAc                                       >  8:158935/1‑124 (MQ=255)
                                                                                                                                      cTTGCGGCACATTCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGc           <  6:307329/149‑1 (MQ=255)
                                                                                                                                               caTGCTGGAATGTCAAGCGCTGTTATTGGTTTTGTCGTCGGTGACCGTCAGCCGGGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGGTATCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAgg   <  1:252471/149‑1 (MQ=255)
                                                                                                                                               caTGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCAAggg  <  4:119755/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TAGTACAGGTGCGTACCACAGCGGCCACAGAAAGCACGCTGGGCCCAGTCGGACGAATCATAGACGCTGGGCGCACGCCCTTCGATCAGCGGAGGCTTGCTGCAGTGCACCACCAGCAGCGGGCCGCCGGTCCACTTGCGGCACATGCTGCAATGACAAGCGCTGATATGGGTATTGTCGACGGTGACCGTCAGCCGCGTGTCGCCGCACAGGCAGGTGCCGTGTTTTTCCAGG‑GCCATGGGTGGCTCCTTGAACGGAAAGATGGGCCTGCAATTATAGAAGCCACCCAGGG  >  NC_002947/4905703‑4905994

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