Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,330 A→C 7.1% D286A (GAC→GCC)  gltA → citrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,3300AC 7.1% 208.6 / 6.4 70D286A (GAC→GCC) gltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (32/33);  new base C (2/3);  total (34/36)
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.69e-01

CACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAGAACGAC  >  NC_002947/4740184‑4740478
                                                                                                                                                    |                                                                                                                                                    
caccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGAc                                                                                                                                                     >  5:355087/1‑148 (MQ=255)
caccGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACa                                                                                                                                                    <  3:54510/149‑1 (MQ=255)
        gTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAATATCGACAAATTCATc                                                                                                                                            <  4:282674/149‑1 (MQ=255)
         gcttgCAGGCTCGTCGGGTGCCACCCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTTGGGCCGGTGCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATCAAATCGGCTATGTCTCGAACATCGACAAATTCATCg                                                                                                                                           <  8:296866/145‑1 (MQ=255)
            ggCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgcc                                                                                                                                         >  6:247082/1‑148 (MQ=255)
             gCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGTTCGAGACATCGCCGATTTCATcca                                                                                                                                          <  4:348323/146‑3 (MQ=37)
                ggCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgccaaggc                                                                                                                                    <  2:133240/149‑1 (MQ=255)
                  cTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGTTCGAGACATCGCCGATTTCATcca                                                                                                                                          <  8:345695/141‑3 (MQ=37)
                  cTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGTTCGAGACATCGCCGATTTCATcca                                                                                                                                          >  7:345695/1‑139 (MQ=37)
                   tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGTTCGAGACATCGCCGATTTCATccagcatg                                                                                                                                     <  2:148152/145‑8 (MQ=25)
                   tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAATTTGTCGATGTTCGAGACATCGCCGATTTCATccagcatg                                                                                                                                     >  1:148152/1‑138 (MQ=25)
                   tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAAt                                                                                                                                                 <  3:418845/133‑1 (MQ=255)
                   tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAAt                                                                                                                                                 >  4:418845/1‑133 (MQ=255)
                   tcgtcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAAGTCATCgccaaggcc                                                                                                                                   >  8:295858/1‑147 (MQ=255)
                     gtcgTGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgccaaggccaag                                                                                                                                >  3:231877/1‑148 (MQ=255)
                      tcgGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGa                                                                                                                              >  2:1205/1‑149 (MQ=255)
                             ccAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATc                                                                                                                                            >  7:282576/1‑128 (MQ=255)
                             ccAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATc                                                                                                                                            <  8:282576/128‑1 (MQ=255)
                                aCCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGAc                                                                                                                    <  2:265743/149‑1 (MQ=255)
                                        gCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTc                                                                                                              >  1:229661/1‑147 (MQ=255)
                                        gCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCaa                                                                                                            <  4:390568/149‑1 (MQ=255)
                                            gTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACaaa                                                                                                                                                  >  5:132170/1‑107 (MQ=255)
                                            gTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACaaa                                                                                                                                                  <  6:132170/107‑1 (MQ=255)
                                            gTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgc                                                                                                                                          <  6:427988/115‑1 (MQ=255)
                                            gTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgc                                                                                                                                          >  5:427988/1‑115 (MQ=255)
                                             tATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGa                                                                                                       <  7:391829/149‑1 (MQ=255)
                                             tATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGa                                                                                                       <  7:371307/149‑1 (MQ=255)
                                               tCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgccaaggccaa                                                                                                                                 <  8:6564/121‑1 (MQ=255)
                                               tCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCgccaaggccaa                                                                                                                                 >  7:6564/1‑121 (MQ=255)
                                               tCGCAGACGTTATCGAAGAACTGTGGGTCCAGGCCAAATTCAGTGCGAACGAAGACGTACTGACCATGCTGGATGTAATCGGCGATGTCTCGAACATCGACAAATTCATCGACAAGGCAAAGGACAAGAACGACCCCTTCTAGCTGAtc                                                                                                     >  4:70823/1‑148 (MQ=255)
                                                cgccgcCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGACAAGGCCAAGGACAAGAACGACCCGTTCAAGCGGATgg                                                                                                    >  8:290221/1‑149 (MQ=255)
                                                   cgccgGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCCAGCTGATGGGCt                                                                                                 >  6:202496/1‑149 (MQ=255)
                                                    gccgGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCtt                                                                                                <  1:98589/149‑1 (MQ=255)
                                                      cgGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCg                                                                                              <  2:290653/149‑1 (MQ=255)
                                                        gTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGt                                                                                            >  1:317769/1‑149 (MQ=255)
                                                         tATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGt                                                                                            >  1:420081/1‑148 (MQ=255)
                                                           tCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCa                                                                                          <  2:342488/148‑1 (MQ=255)
                                                              ccgcACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCGAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACcgc                                                                                      >  7:141271/1‑149 (MQ=255)
                                                                gcACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCgt                                                                                    <  2:229661/149‑1 (MQ=255)
                                                                   gccgtgggccgggccCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCgtg                                                                                   >  6:431244/6‑148 (MQ=255)
                                                                       ggggCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTGAAGCTGATGGGCTTCGGTCACCGCGTGTAc                                                                                >  1:116012/1‑146 (MQ=255)
                                                                       ggggCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAg                                                                             <  5:15440/149‑1 (MQ=255)
                                                                         ggcccggcccACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCCAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGaa                                                                           >  7:102704/1‑149 (MQ=255)
                                                                          gcccggcccACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAAc                                                                          >  6:431162/1‑149 (MQ=255)
                                                                                   cACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGAccc                                                                  >  8:197359/1‑149 (MQ=255)
                                                                                             cgaacgaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCAcc                                                        >  6:36195/1‑149 (MQ=255)
                                                                                                 cgaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtga                                                    <  1:71911/149‑1 (MQ=255)
                                                                                                  gaaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtgat                                                   >  6:376579/1‑149 (MQ=255)
                                                                                                   aaGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGtgatg                                                  <  5:247082/149‑1 (MQ=255)
                                                                                                      ccGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAg                                               <  7:197359/149‑1 (MQ=255)
                                                                                                          aCTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGa                                           >  8:364174/1‑149 (MQ=255)
                                                                                                           cTGACCATGCTCGTTGAAATCCGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTTTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGAc                                          <  4:238617/149‑1 (MQ=255)
                                                                                                              aCCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCt                                        <  5:431162/148‑1 (MQ=255)
                                                                                                               ccATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGc                                      <  5:431244/149‑1 (MQ=255)
                                                                                                                cATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcg                                     >  3:33868/1‑149 (MQ=255)
                                                                                                                 aTGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcga                                    <  7:364174/149‑1 (MQ=255)
                                                                                                                    cTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGcgacg                                  <  2:317769/148‑1 (MQ=255)
                                                                                                                     tGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGaa                                <  2:116012/149‑1 (MQ=255)
                                                                                                                       gATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCCACGaa                                <  7:392771/147‑1 (MQ=255)
                                                                                                                          gAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCt                           >  3:213457/1‑149 (MQ=255)
                                                                                                                           aaaTCGGCGTTGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTg                          >  6:411107/1‑149 (MQ=255)
                                                                                                                                  cGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGc                   <  2:420081/149‑1 (MQ=255)
                                                                                                                                  cGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGc                   >  8:219437/1‑149 (MQ=255)
                                                                                                                                        ctcGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGCGCTTCGGTCACCGCGTGTCCAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTCCGACGAAGTCCTGCGCGAGCTGGGc              >  7:280407/1‑148 (MQ=255)
                                                                                                                                           gAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCa             <  1:1205/146‑1 (MQ=255)
                                                                                                                                           gAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATc           <  5:376579/148‑1 (MQ=255)
                                                                                                                                            aaCATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCaa         <  4:33868/149‑1 (MQ=255)
                                                                                                                                               cATCGACAAATTCATCGCCAAGGGCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCCCCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAg        >  4:60781/1‑148 (MQ=255)
                                                                                                                                               cATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAg        >  8:164736/1‑148 (MQ=255)
                                                                                                                                                aTCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAGaa      <  3:70823/149‑1 (MQ=255)
                                                                                                                                                    aCAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAGAACGAc  <  4:231877/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACCGTTCGTCTGGCAGGCTCGTCGGGTGCCAACCCGTTCGCCTGTATCGCCGCCGGTATCGCCGCACTGTGGGGCCCGGCCCACGGCGGTGCGAACGAAGCCGTACTGACCATGCTGGATGAAATCGGCGATGTCTCGAACATCGACAAATTCATCGCCAAGGCCAAGGACAAGAACGACCCGTTCAAGCTGATGGGCTTCGGTCACCGCGTGTACAAGAACCGCGACCCGCGCGCCACCGTGATGAAGCAGACCTGCGACGAAGTCCTGCGCGAGCTGGGCATCAAGAACGAC  >  NC_002947/4740184‑4740478

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: