Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 81,984 G→C 10.0% E98D (GAG→GAC PP_0071 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294781,9840GC10.0% 213.5 / 13.3 70E98D (GAG→GACPP_0071hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (3/4);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

CACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTA  >  NC_002947/81835‑82133
                                                                                                                                                     |                                                                                                                                                     
cACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAg                                                                                                                                                       >  1:708434/1‑150 (MQ=255)
    tcgCCCTATACCTGCCCAGTGCTTACGAAACCGGCTTATCAGCTTTGACGTGTCCTCAGAATCAAGACCTTTGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGa                                                                                                                                                   <  1:444682/148‑1 (MQ=255)
       cccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGc                                                                                                                                                >  1:779612/1‑150 (MQ=255)
        ccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGcc                                                                                                                                               >  2:159094/1‑150 (MQ=255)
        ccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGcc                                                                                                                                               >  2:678786/1‑150 (MQ=255)
            taCCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCg                                                                                                                                           <  1:159094/150‑1 (MQ=255)
               cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCgagt                                                                                                                                        <  2:228810/150‑4 (MQ=255)
               cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCg                                                                                                                                        >  2:107594/1‑150 (MQ=255)
                   gcAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCAc                                                                                                                                    >  2:1073612/1‑150 (MQ=255)
                        gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                               <  2:708434/150‑1 (MQ=255)
                        gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                               <  1:690636/150‑1 (MQ=255)
                              gAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTcc                                                                                                                         >  2:460017/1‑150 (MQ=255)
                                   cGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGc                                                                                                                    >  1:435646/1‑150 (MQ=255)
                                         cAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcg                                                                                                              >  2:746647/1‑150 (MQ=255)
                                          aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCacg                                                                                                                      <  2:323062/141‑3 (MQ=255)
                                          aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                      >  1:323062/1‑141 (MQ=255)
                                              cTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCgctgct                                                                                                         <  1:1073612/150‑1 (MQ=255)
                                               tttGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTa                                                                                                        >  2:208938/1‑150 (MQ=255)
                                                 tGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATc                                                                                                      <  1:857421/150‑1 (MQ=255)
                                                   gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                             <  1:982520/125‑1 (MQ=255)
                                                   gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                             >  2:982520/1‑125 (MQ=255)
                                                     cTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAgg                                                                                                  <  1:107594/150‑1 (MQ=255)
                                                      tACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGc                                                                                                 <  1:242382/150‑1 (MQ=255)
                                                             gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                          <  2:779612/150‑1 (MQ=255)
                                                             gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                          <  1:108382/150‑1 (MQ=255)
                                                             gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGCCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                          >  1:521799/1‑150 (MQ=255)
                                                              tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                         <  2:457561/150‑1 (MQ=255)
                                                              tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                         >  1:899187/1‑150 (MQ=255)
                                                              tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGCAATAGCGGGCTTGCCCGCGCCTGCACTCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                         >  1:444914/1‑150 (MQ=255)
                                                                  gaagGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTAcc                                                                                     >  1:219797/1‑150 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                                  <  2:1132386/86‑1 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                                  >  1:1132386/1‑86 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCCAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCtgt                                                                                  <  2:345896/150‑1 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGACGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCtgt                                                                                  >  2:424364/1‑150 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGACGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATa                                                                                                                 <  1:904166/119‑1 (MQ=255)
                                                                     ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGACGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATa                                                                                                                 >  2:904166/1‑119 (MQ=255)
                                                                            gCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                      >  2:70014/1‑107 (MQ=255)
                                                                            gCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                      <  1:70014/107‑1 (MQ=255)
                                                                             cATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                          <  1:460017/150‑1 (MQ=255)
                                                                             cATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                          <  2:899187/150‑1 (MQ=255)
                                                                                      gTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCCTCGAGTCGACGCCATAGCCGTCAGGCCCAGCCTTCAAGGCAGAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCATGCCcg                                                                 >  1:537080/1‑150 (MQ=255)
                                                                                               gCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGCGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCg                                                        <  1:1150073/150‑1 (MQ=255)
                                                                                               gCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCg                                                        <  1:678786/150‑1 (MQ=255)
                                                                                                 aTGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                      >  2:549708/1‑150 (MQ=255)
                                                                                                  tGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCa                                                     >  1:428513/1‑150 (MQ=255)
                                                                                                   gCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCaa                                                    >  2:47372/1‑150 (MQ=255)
                                                                                                    ggtttGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCaaa                                                   <  1:369720/147‑1 (MQ=255)
                                                                                                     gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTGGCCGCCATCGCCGTCAAGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTCACTGTGGCAGCGGGCGTGCCCGCGAAAACCGGCAAAg                                                  >  2:369720/1‑150 (MQ=255)
                                                                                                     gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                      <  1:1061094/82‑1 (MQ=255)
                                                                                                     gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                      >  2:1061094/1‑82 (MQ=255)
                                                                                                     gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAg                                                  >  1:245581/1‑150 (MQ=255)
                                                                                                      aTTGAGGTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGc                                                 >  2:1130884/1‑150 (MQ=255)
                                                                                                          agaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGt                                             <  2:219797/150‑1 (MQ=255)
                                                                                                           gaTGGCAGAGGGTGTGATCGAGGCGATGCTCGAGCGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                            <  2:427868/150‑1 (MQ=255)
                                                                                                           gaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                            >  2:956609/1‑150 (MQ=255)
                                                                                                                         tgATCGAGGCGATCCTCCATGGCATCGAGTCGACGCCAGCGCCGTCACGCCCAGCCTCCAAGGTAGATTACTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGt                              <  1:208938/150‑1 (MQ=37)
                                                                                                                          gATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTg                             <  2:444914/150‑1 (MQ=255)
                                                                                                                              gAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACCCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGctt                         >  2:1055099/1‑150 (MQ=255)
                                                                                                                                ggCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGcttct                       <  2:428513/150‑1 (MQ=255)
                                                                                                                                   gATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttcc                    >  2:406183/1‑149 (MQ=255)
                                                                                                                                    aTACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTcgc                   <  1:746647/150‑1 (MQ=255)
                                                                                                                                     tACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTGCTTcgcg                  >  1:486617/1‑150 (MQ=255)
                                                                                                                                          gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                          <  1:630010/89‑1 (MQ=39)
                                                                                                                                          gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                          >  2:630010/1‑89 (MQ=39)
                                                                                                                                          gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                         <  2:1004094/122‑1 (MQ=255)
                                                                                                                                          gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                         >  1:1004094/1‑122 (MQ=255)
                                                                                                                                          gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc                     <  2:58548/142‑1 (MQ=255)
                                                                                                                                          gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc                     >  1:58548/1‑142 (MQ=255)
                                                                                                                                          gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTaa             >  2:431936/1‑150 (MQ=255)
                                                                                                                                               cATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCg        >  1:637175/1‑150 (MQ=255)
                                                                                                                                                   gAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTcc    <  1:431936/150‑1 (MQ=255)
                                                                                                                                                     gTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTACACCCGCTCCTa  >  1:872540/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTA  >  NC_002947/81835‑82133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: