Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,555,144 A→G 6.1% N406D (AAC→GAC)  sbcB → exodeoxyribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,555,1440AG 6.1% 192.0 / 2.2 65N406D (AAC→GAC) sbcBexodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAG  >  NC_002947/1554997‑1555290
                                                                                                                                                   |                                                                                                                                                  
cagcagTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTaa                                                                                                                                                   <  3:289686/149‑1 (MQ=255)
 agcagTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAAc                                                                                                                                                  <  6:324302/149‑1 (MQ=255)
     gTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                                                                                                              >  3:83608/1‑149 (MQ=255)
      ttGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTccc                                                                                                                                             >  2:72486/1‑149 (MQ=255)
            gAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                                                                                                       >  8:34563/1‑149 (MQ=255)
            gAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                                                                                                       >  6:233074/1‑149 (MQ=255)
            gAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAAc                                                                                                                                       >  6:227730/1‑149 (MQ=255)
                     ataGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCgg                                                                                                                              >  8:205237/1‑149 (MQ=255)
                            aCCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAagag                                                                                                                       <  5:338471/149‑1 (MQ=255)
                              cgcgACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCtt                                                                                                                                    <  1:93213/134‑1 (MQ=255)
                              cgcgACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCtt                                                                                                                                    >  2:93213/1‑134 (MQ=255)
                                cgACCGGCGCTTAAGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTAGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTAGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCgg                                                                                                                    >  5:303364/1‑148 (MQ=255)
                                     ggCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGt                                                                                                              >  8:214981/1‑149 (MQ=255)
                                      gcgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGTCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCCACGtt                                                                                                             >  4:159501/1‑149 (MQ=255)
                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                                                                                                              <  8:259056/115‑1 (MQ=255)
                                       cgcTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                                                                                                              >  7:259056/1‑115 (MQ=255)
                                                gaccAAGGTCGGGCCCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCtt                                                                                                   <  7:161269/146‑1 (MQ=255)
                                                      gTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCGAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGGTTCTGcag                                                                                             >  5:105000/1‑147 (MQ=255)
                                                      gTTAGTGCTCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTg                                                                                                <  6:105000/146‑1 (MQ=255)
                                                          gTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACa                                                                                          <  4:83608/148‑1 (MQ=255)
                                                          gTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAg                                                                                         <  7:205237/149‑1 (MQ=255)
                                                          gTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCACTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACa                                                                                          <  4:48425/148‑1 (MQ=255)
                                                                  gAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGa                                                                                  <  2:378186/148‑1 (MQ=255)
                                                                  gAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGa                                                                                  <  5:233074/148‑1 (MQ=255)
                                                                  gAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGa                                                                                  <  5:227730/148‑1 (MQ=255)
                                                                              ttGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTg                                                                     >  6:273067/1‑149 (MQ=255)
                                                                                gAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGTGCCTGAGCGAACCGCAGTggg                                                                   >  5:349318/1‑149 (MQ=255)
                                                                                gAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTgg                                                                    <  7:214981/148‑1 (MQ=255)
                                                                                 aGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTGACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTggg                                                                   <  2:257370/148‑1 (MQ=255)
                                                                                 aGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTGACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTggg                                                                   >  1:257370/1‑148 (MQ=255)
                                                                                 aGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTGACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTggg                                                                   >  1:252071/1‑148 (MQ=255)
                                                                                 aGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTGACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTggg                                                                   <  2:252071/148‑1 (MQ=255)
                                                                                   ccATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTAt                                                                                                        >  2:321738/1‑109 (MQ=255)
                                                                                   ccATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTAt                                                                                                        <  1:321738/109‑1 (MQ=255)
                                                                                    cATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGgc                                                               <  7:34563/149‑1 (MQ=255)
                                                                                            gTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAg                                                                       >  4:10420/1‑133 (MQ=255)
                                                                                            gTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAg                                                                       <  3:10420/133‑1 (MQ=255)
                                                                                               gATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCt                                                    >  4:49104/1‑149 (MQ=255)
                                                                                                      gatgatCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGa                                             >  1:253060/1‑149 (MQ=255)
                                                                                                          atCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                                                                                                              >  4:274872/1‑48 (MQ=255)
                                                                                                          atCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTcc                                                                                                                                              <  3:274872/48‑1 (MQ=255)
                                                                                                              gcgTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTgg                                                 <  8:330882/137‑1 (MQ=255)
                                                                                                              gcgTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTgg                                                 >  7:330882/1‑137 (MQ=255)
                                                                                                              gcgTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGaa                                      >  8:207713/1‑148 (MQ=255)
                                                                                                              gcgTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGaa                                      >  8:119466/1‑148 (MQ=255)
                                                                                                                 tttGCCGGAGCTGTTGTTCCGCTACCGGTCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTg                                                                     <  4:211134/114‑1 (MQ=255)
                                                                                                                 tttGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTg                                                                     >  3:211134/1‑114 (MQ=255)
                                                                                                                   ggCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGgc                                 <  8:240837/147‑1 (MQ=255)
                                                                                                                            tgttgtTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCaa                       >  3:69624/1‑149 (MQ=255)
                                                                                                                             gttgttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCtt                                                   <  1:46874/120‑1 (MQ=255)
                                                                                                                             gttgttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCtt                                                   >  2:46874/1‑120 (MQ=255)
                                                                                                                             gttgttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAg                      >  7:359675/1‑149 (MQ=255)
                                                                                                                              ttgttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGc                     >  6:306966/1‑149 (MQ=255)
                                                                                                                               tgttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGcc                    >  7:53029/1‑149 (MQ=255)
                                                                                                                                 ttCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCt                   <  1:72486/148‑1 (MQ=255)
                                                                                                                                   ccGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCCACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTgg                 <  3:159501/148‑1 (MQ=255)
                                                                                                                                   ccGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTgg                 <  7:207713/148‑1 (MQ=255)
                                                                                                                                    cGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGa               <  3:49104/149‑1 (MQ=255)
                                                                                                                                    cGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGa               <  7:119466/149‑1 (MQ=255)
                                                                                                                                    cGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCCGCGAAGAGCGACAAAGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCTAACCGAAGTGGGGGGAGCCGAATACCTTGGGCGACTTCGAAAAGCCGCGCCAGCAAGCCTgg                 >  7:288781/1‑147 (MQ=255)
                                                                                                                                     gCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAg              >  6:33505/1‑149 (MQ=255)
                                                                                                                                         ccGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCTGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGgc           <  4:69624/148‑1 (MQ=255)
                                                                                                                                               aCGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCcgac     >  6:77268/1‑148 (MQ=255)
                                                                                                                                                cGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCcgacga   <  5:77268/149‑1 (MQ=255)
                                                                                                                                                 gTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTACGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCcgacgaa  >  7:245412/1‑148 (MQ=255)
                                                                                                                                                 gTAACTTCCCTGAAACCTTGACCGGAGAAGATCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCcgacg    <  8:245412/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCAGTTGTATGAGGGTTTTATAGGTGACCGCGACCGGCGCTTATGCGAGCAAGTTCGTGCGCTGGAACCGGCGCAGTTGAGCCATGGCCAGTGGATGTTCGATGATCCGCGTTTGCCGGAGCTGTTGTTCCGCTACCGGGCACGTAACTTCCCTGAAACCTTGACCGGCGAAGAGCGGCAACGTTGGTATGGCTTCTGCCAACAGCGCCTGAGCGAACCGCAGTGGGGGGCGCCGAATACCTTGGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAG  >  NC_002947/1554997‑1555290

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