Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 220,936 C→G 6.2% A89G (GCG→GGG)  PP_5431 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947220,9360CG 6.2% 205.6 / 4.9 65A89G (GCG→GGG) PP_5431hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

CCAGGAATGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTATTTAAATT  >  NC_002947/220790‑221084
                                                                                                                                                  |                                                                                                                                                     
ccAGGAATGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCa                                                                                                                                                    <  1:873251/150‑1 (MQ=255)
  aGGAATGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCa                                                                                                                                                    >  1:409513/1‑148 (MQ=255)
  aGGAATGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCa                                                                                                                                                    <  2:409513/148‑1 (MQ=255)
     aaTGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCa                                                                                                                                               >  2:77139/1‑150 (MQ=255)
      aTGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAggtcgg                                                                                                                                              <  2:514311/150‑5 (MQ=255)
          tttGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAgctg                                                                                                                                          >  2:895583/1‑150 (MQ=255)
                gaggatgTCGTGAGCAACAGCCGCCGCCATCTCAACCAAACACACCTGAACACGACCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAAc                                                                                                                                    <  2:881935/147‑1 (MQ=255)
                     tgTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTc                                                                                                                               <  2:1177158/150‑1 (MQ=255)
                         gTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACg                                                                                                                           >  1:1034028/1‑150 (MQ=255)
                         gTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACg                                                                                                                           >  2:166121/1‑150 (MQ=255)
                            gagaACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGgcg                                                                                                                        >  2:907412/1‑150 (MQ=255)
                                        aaCCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTc                                                                                                            <  1:166121/150‑1 (MQ=255)
                                           cAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTcagc                                                                                                         <  1:917442/150‑1 (MQ=255)
                                           cAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTcagc                                                                                                         <  1:895583/150‑1 (MQ=255)
                                               tCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTcagcttcc                                                                                                     <  1:44421/150‑5 (MQ=255)
                                                   acaaaGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCaa                                                                                                 >  2:907098/1‑150 (MQ=255)
                                                     aaaGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCGCCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATa                                                                                               >  1:785499/1‑150 (MQ=255)
                                                       aGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAAc                                                                                             <  2:1034028/150‑1 (MQ=255)
                                                              tGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTg                                                                                      >  1:1014521/1‑150 (MQ=255)
                                                              tGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTg                                                                                      >  2:1169830/1‑150 (MQ=255)
                                                                gACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTTACCGCCAAGGGCGTATTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACCGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGAc                                                                                    >  2:348821/1‑150 (MQ=255)
                                                                 aCATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACg                                                                                   >  2:457074/1‑150 (MQ=255)
                                                                  cATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACga                                                                                  >  1:193619/1‑150 (MQ=255)
                                                                   aTGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTc                                                                                                                                       <  1:566937/96‑1 (MQ=255)
                                                                   aTGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTc                                                                                                                                       >  2:566937/1‑96 (MQ=255)
                                                                    tGATCAGGGTGCTCGAGGACCTGATCGACGTGTTGACCAACAAGCGCGTGTTCAACATGACCGACCTGCCCTCCGGCGCTCATGACAAGCATATCAATCATTGCTGGGCGAGGAAGGCGATCAGCTGCCTCAATAACCTGATTGACgtag                                                                                <  2:193619/150‑1 (MQ=255)
                                                                     gATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTcag                                                                                                          <  2:47711/123‑1 (MQ=255)
                                                                     gATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTcag                                                                                                          >  1:47711/1‑123 (MQ=255)
                                                                          gcgTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAgc                                                                                                                                            >  2:452555/1‑84 (MQ=255)
                                                                          gcgTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAgc                                                                                                                                            <  1:452555/84‑1 (MQ=255)
                                                                          gcgTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCACTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGc                                                                          <  2:423870/150‑1 (MQ=255)
                                                                              gCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCGCCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAggg                                                                      <  2:785499/150‑1 (MQ=255)
                                                                              gCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAggg                                                                      >  1:987445/1‑150 (MQ=255)
                                                                                tCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAgggcg                                                                    >  2:584163/1‑150 (MQ=255)
                                                                                  gAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACAGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAgggcggg                                                                  >  1:987174/1‑150 (MQ=255)
                                                                                    ggACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCa                                                                                                                             >  2:1250587/1‑89 (MQ=255)
                                                                                    ggACCGGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCa                                                                                                                             <  1:1250587/89‑1 (MQ=255)
                                                                                        cTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGAta                                                            >  2:1018447/1‑149 (MQ=255)
                                                                                         tGATCGACGTGTTCACAGCCAAGAGCGAGTTCAAGATCACCGACCTGCCCGGCGGCGCGCAGGCCAAGCAGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCt                                                          >  2:531284/1‑150 (MQ=255)
                                                                                          gATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAgg                                                                       <  2:96583/137‑1 (MQ=255)
                                                                                          gATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAgg                                                                       >  1:96583/1‑137 (MQ=255)
                                                                                           aTCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGa                                                        <  1:1169830/150‑1 (MQ=255)
                                                                                              gACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTg                                                     <  2:987445/150‑1 (MQ=255)
                                                                                              gACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCGCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTg                                                     <  2:987174/150‑1 (MQ=255)
                                                                                                    ttGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGa                                                        >  2:320918/1‑141 (MQ=255)
                                                                                                    ttGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGa                                                        <  1:320918/141‑1 (MQ=255)
                                                                                                       aCCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCt                                            >  2:452022/1‑150 (MQ=255)
                                                                                                           ccAAGGGCGTGTTAAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCcg                                        >  2:1149437/1‑150 (MQ=255)
                                                                                                            cAAGGGCGTGTTCAAGATCACCGACTTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCcgc                                       >  1:726685/1‑150 (MQ=255)
                                                                                                                   gtgtTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGTATCCTCAATAACCTGATTCACGACGAAGGACATCGTATGCTGATTTGACTGCAGGACGCTCCCGCGACAGTg                                >  1:287072/1‑150 (MQ=255)
                                                                                                                      ttCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGa                                                <  2:139433/131‑1 (MQ=255)
                                                                                                                      ttCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGa                                                >  1:139433/1‑131 (MQ=255)
                                                                                                                         aaGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTg                          >  1:816920/1‑150 (MQ=255)
                                                                                                                         aaGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTg                          <  1:907412/150‑1 (MQ=255)
                                                                                                                               aCCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTc                    <  1:77139/150‑1 (MQ=255)
                                                                                                                                ccGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTAtt                     >  2:951479/1‑148 (MQ=255)
                                                                                                                                ccGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTAtt                     <  1:951479/148‑1 (MQ=255)
                                                                                                                                    ccTGCGCGCCGGCGGGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTaa               >  2:759211/1‑150 (MQ=255)
                                                                                                                                    ccTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCATCAGCCTCAATAAACTGATTGACGTCGAAATCAAGGAGTCTCGTTGCGTGTTGCTGGACGCTCCCGCGACAGGGATATTGTTATTCCTTaa               <  1:348821/150‑1 (MQ=255)
                                                                                                                                     cTGCGCGCCGGCGGGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAg                                                                        <  2:1248152/93‑1 (MQ=255)
                                                                                                                                     cTGCGCGCCGGCGGGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAg                                                                        >  1:1248152/1‑93 (MQ=255)
                                                                                                                                     cTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTc                                                                                                            <  2:174980/57‑1 (MQ=255)
                                                                                                                                     cTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTc                                                                                                            >  1:174980/1‑57 (MQ=255)
                                                                                                                                       gcgcgccgGCGGGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATgg            <  1:759211/147‑1 (MQ=255)
                                                                                                                                          cgccgGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTAt         <  2:816920/150‑1 (MQ=255)
                                                                                                                                               gcgcgcAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCACAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTATTTaaa    <  2:756154/150‑1 (MQ=255)
                                                                                                                                                cgcgcAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTATTTAAat   <  2:726685/150‑1 (MQ=255)
                                                                                                                                                 gcgcAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTATTTAAatt  >  1:691130/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAGGAATGGTTTGCCGATGATGTCGTGGAGAACAGCCTCAACCAGCTCAAACAAAGCGACCTGGACATGATCCGCGTGCTCGAGGACCTGATCGACGTGTTGACCGCCAAGGGCGTGTTCAAGATCACCGACCTGCCCGCCGGCGCGCAGGCCAAGCTGCTCAACCGTTCCACGGCGCGCAAGGCGCTCAGCAGCCTCAATAACCTGATTGACGAAGAAGAGCAGGGCGGGCTGATCTGACTGCAGGACGCTCCCGCGACAGTGATATTGTTATTCCTTAATGGTATTTAAATT  >  NC_002947/220790‑221084

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