Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 81,979 A→G 9.8% I97V (ATC→GTC)  PP_0071 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294781,9790AG 9.8% 192.3 / 10.1 72I97V (ATC→GTC) PP_0071hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (3/4);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

TGGTCCACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCG  >  NC_002947/81830‑82127
                                                                                                                                                     |                                                                                                                                                    
tGGTCCACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCa                                                                                                                                                      >  2:857421/1‑150 (MQ=255)
     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                                                                                                                         <  1:690636/150‑1 (MQ=255)
                             gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                         <  2:708434/150‑1 (MQ=255)
                                   gAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATc                                                                                                                                                    <  1:220170/117‑1 (MQ=255)
                                   gAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATc                                                                                                                                                    >  2:220170/1‑117 (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                                                                                                                       >  2:982520/1‑125 (MQ=255)
                                                        gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                       <  1:982520/125‑1 (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                                                                    >  2:630010/1‑89 (MQ=39)
                                                                                                                                               gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                    <  1:630010/89‑1 (MQ=39)
                                                                                                                                               gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                   <  2:1004094/122‑1 (MQ=255)
                                                                                                                                               gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                   >  1:1004094/1‑122 (MQ=255)
                                                                                                                                               gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc               >  1:58548/1‑142 (MQ=255)
                                                                                                                                               gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc               <  2:58548/142‑1 (MQ=255)
                                                                                                                                               gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTaa       >  2:431936/1‑150 (MQ=255)
                                                                                                                                                    cATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCg  >  1:637175/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTCCACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCG  >  NC_002947/81830‑82127

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