Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,815,544 A→G 10.0% intergenic (+28/+63) PP_5090 → / ← PP_5091 cell division protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,815,5440AG10.0% 155.5 / 7.2 60intergenic (+28/+63)PP_5090/PP_5091cell division protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (3/3);  total (30/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 = 8.75e-01

CTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAACCCGCCGG  >  NC_002947/5815396‑5815692
                                                                                                                                                    |                                                                                                                                                    
cTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGa                                                                                                                                                      <  6:134704/149‑1 (MQ=255)
    tACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACACAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                                  <  1:335001/149‑1 (MQ=255)
      ccGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCAAg                                                                                                                                                  >  3:37999/1‑147 (MQ=255)
      ccGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGTCGCATCGCTGCGGCAAg                                                                                                                                                  <  4:37999/147‑1 (MQ=255)
            gCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTCGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCAGCAGCGATGCGACAAGGCcttctt                                                                                                                                          >  1:372449/1‑149 (MQ=255)
                gTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTtata                                                                                                                                      <  2:339414/149‑1 (MQ=255)
                      ccGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCTGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACa                                                                                                                                <  3:313615/149‑1 (MQ=255)
                        gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                                  <  5:205516/129‑1 (MQ=255)
                        gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAg                                                                                                                                                  >  6:205516/1‑129 (MQ=255)
                        gTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                                              <  5:325906/149‑1 (MQ=255)
                            aGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTg                                                                                                                          >  5:229021/1‑149 (MQ=255)
                             gCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGc                                                                                                                         >  6:159399/1‑149 (MQ=255)
                                 ccGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGt                                                                                                                     <  6:75210/149‑1 (MQ=255)
                                    cgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                                              >  7:325585/1‑137 (MQ=255)
                                    cgAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGc                                                                                                                              <  8:325585/137‑1 (MQ=255)
                                     gAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                                 <  2:183232/149‑1 (MQ=255)
                                      aGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                                 >  4:152472/1‑148 (MQ=255)
                                      aGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAg                                                                                                                 <  3:152472/148‑1 (MQ=255)
                                          gTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCt                                                                                                            <  4:59382/149‑1 (MQ=255)
                                                aGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAag                                                                                                      <  6:317971/149‑1 (MQ=255)
                                                      acaGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGccc                                                                                                >  5:213309/1‑149 (MQ=255)
                                                             cagcTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCAcc                                                                                         >  5:80103/1‑149 (MQ=255)
                                                                   gctggGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgc                                                                                   <  6:80103/149‑1 (MQ=255)
                                                                         gCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccgcc                                                                             <  6:229021/149‑1 (MQ=255)
                                                                         gCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAgccccgcc                                                                             >  6:305342/1‑149 (MQ=255)
                                                                            gggTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGccccgccccg                                                                          >  2:207483/1‑149 (MQ=255)
                                                                             ggTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGccccgccccg                                                                          >  2:354664/1‑148 (MQ=255)
                                                                                    cAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCtt                                                                  >  1:358150/1‑149 (MQ=255)
                                                                                     aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                                 >  8:83533/1‑149 (MQ=255)
                                                                                     aaCCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                                 <  1:367301/149‑1 (MQ=255)
                                                                                      cgcTGTTGCAGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTc                                                                     <  5:204705/142‑1 (MQ=255)
                                                                                      cgcTGTTGCAGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTc                                                                     >  6:204705/3‑144 (MQ=255)
                                                                                        cTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCg                                                               >  5:426231/1‑148 (MQ=255)
                                                                                           gtgctgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGa                                                             <  4:110546/146‑1 (MQ=255)
                                                                                              ctgcAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACAcc                                                         <  5:305342/148‑1 (MQ=255)
                                                                                                 cAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                                      <  1:207483/148‑1 (MQ=255)
                                                                                                  aaCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGa                                                                                                                   <  2:38396/86‑1 (MQ=255)
                                                                                                  aaCAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGa                                                                                                                   >  1:38396/1‑86 (MQ=255)
                                                                                                       cGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTc                                                >  8:305267/1‑148 (MQ=255)
                                                                                                          cAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACg                                            >  5:60521/1‑149 (MQ=255)
                                                                                                             aCCCGTCAGTAACGAAGAAGGCCTTGCAGCAGCGATGCCACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCACCGCCCCGCATTCCTTTCGGACCCCTGGCTGTTCAACGTc                                          >  4:342868/1‑148 (MQ=255)
                                                                                                              cccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTcc                                         <  6:213309/148‑1 (MQ=255)
                                                                                                              cccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                                          >  7:300149/1‑147 (MQ=255)
                                                                                                               ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                                          >  5:16685/1‑146 (MQ=255)
                                                                                                               ccGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCAGTGAAGGGGCTGTAAAGCAGCCCCGCCCCCGTCAGCCCCGCCCCGCATTACTTTCGGACACCTGGCTGTTCAACGta                                          >  5:77523/1‑145 (MQ=255)
                                                                                                                        aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCCCCGTCAGCCCCGCCCCGCATTCATTTCGGCCAACAGGCGGTTCAAAGGCCAGAAGGCATaa                              >  8:385707/1‑148 (MQ=255)
                                                                                                                        aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCCGCCCCCCCCCGCCTTCCTTTCGGACACCTGGCTGTTCCACGTCCAGAAGTCCTAc                              >  3:26589/1‑149 (MQ=255)
                                                                                                                        aCGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATAc                              <  2:358150/149‑1 (MQ=255)
                                                                                                                         cGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACa                             <  6:426231/149‑1 (MQ=255)
                                                                                                                            agaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCa                          >  1:6083/1‑149 (MQ=255)
                                                                                                                             gaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCACGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                                 >  6:257346/1‑109 (MQ=255)
                                                                                                                             gaagGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCACGCCACCGTCAGCCCCGCCCCGCATTCCttt                                                                 <  5:257346/109‑1 (MQ=255)
                                                                                                                               agGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCAcc                        <  7:305267/148‑1 (MQ=255)
                                                                                                                                  ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                                          >  3:399739/1‑127 (MQ=255)
                                                                                                                                  ccTTGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTc                                          <  4:399739/127‑1 (MQ=255)
                                                                                                                                  ccTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCcca                    <  6:16685/149‑1 (MQ=255)
                                                                                                                                    ttGTCGCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGACCAGAAGTCATACAGCACCCccac                   >  7:111884/1‑148 (MQ=255)
                                                                                                                                    ttGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                                      >  6:61680/1‑113 (MQ=255)
                                                                                                                                    ttGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTgg                                                      <  5:61680/113‑1 (MQ=255)
                                                                                                                                         gCATCGCTGCGGCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGACCAGAAGTCATACAGCACCCCCACCAGGaa             <  8:111884/149‑1 (MQ=255)
                                                                                                                                               aTGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGACCCGCATTCCTATCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAAccc       <  2:372449/149‑1 (MQ=255)
                                                                                                                                                    aCAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAACCCGCCgg  <  6:60521/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGTACCGCGTGCTGGTTGGGCCGTTCAGCAACCGCGAGCAGTTGACAGGGGCACAGAAGCAGCTGGCTGGGGCCGGGTTCAGCAACCTGTTGCTGCAACAGCGGCAGACCCGTCAGTAACGAAGAAGGCCTTGCCGCAGCGATGCGACAAGGCCTTCTTTATATGGACAGCATTGCCCGTGAAGGGGCTGTAAAGCAGCCCCGCCACCGTCAGCCCCGCCCCGCATTCCTTTCGGACACCTGGCTGTTCAACGTCCAGAAGTCATACAGCACCCCCACCAGGAACAACCCGCCGG  >  NC_002947/5815396‑5815692

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