Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 858,325 C→G 9.9% A168P (GCA→CCA)  phr ← deoxyribodipyrimidine photolyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947858,3250CG 9.9% 287.7 / 17.7 99A168P (GCA→CCA) phrdeoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base C (37/51);  new base G (5/5);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

AGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGC  >  NC_002947/858177‑858466
                                                                                                                                                    |                                                                                                                                                
aGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGcc                                                                                                                                                 >  1:274935/1‑150 (MQ=255)
  gAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAg                                                                                                                                               <  2:545963/150‑1 (MQ=255)
        tGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAgggcggg                                                                                                                                         <  1:1028724/150‑1 (MQ=255)
         gACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcgggt                                                                                                                                        <  2:2060631/150‑2 (MQ=255)
              gCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCt                                                                                                                                  >  1:1383297/1‑150 (MQ=255)
                gggACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGc                                                                                                                               >  2:55517/1‑150 (MQ=255)
                 ggACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGCCAGGCTGCg                                                                                                                              <  1:1313279/150‑1 (MQ=255)
                 ggACTGTGCCTCAGCTTCACCGGATCGCTAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCg                                                                                                                              <  1:1209162/150‑1 (MQ=255)
                   aCTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGt                                                                                                                            <  2:313799/150‑1 (MQ=255)
                     tgtgCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGc                                                                                                                          <  2:1348768/150‑1 (MQ=255)
                     tgtgCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGc                                                                                                                          >  1:1439300/1‑150 (MQ=255)
                        gCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTCTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGTGGGCAGGCTGCGGTGCAGa                                                                                                                       <  1:1719073/150‑1 (MQ=255)
                         ccTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGAt                                                                                                                      <  1:55517/150‑1 (MQ=255)
                          cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGAta                                                                                                                     >  2:47161/1‑149 (MQ=255)
                          cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATg                                                                                                                     >  1:134040/1‑150 (MQ=255)
                          cTCAGATTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATg                                                                                                                     >  1:215245/1‑150 (MQ=255)
                           tCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGc                                                                                                                    >  1:1503706/1‑150 (MQ=255)
                           tCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGACAGGGCGCGGAGGCCGCGGTTCATATGc                                                                                                                    >  2:137699/1‑150 (MQ=255)
                               cTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCaa                                                                                                                >  1:819353/1‑150 (MQ=255)
                                ttCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGTAGATGCTCAAg                                                                                                               >  1:400351/1‑150 (MQ=255)
                                 tCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGc                                                                                                             <  1:93341/150‑1 (MQ=255)
                                  cACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCAGATGTTCAGCGGGGCCTGGCTCAACACCCGCCGTTCCAGTGCGGGCAGGCTGCGGTGCAGATGCTCAAGGc                                                                                                             <  1:47161/150‑1 (MQ=255)
                                   aCCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCa                                                                                                            <  2:274935/150‑1 (MQ=255)
                                   aCCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcggg                                                                                                                                         <  1:1106576/123‑1 (MQ=255)
                                   aCCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcggg                                                                                                                                         >  2:1106576/1‑123 (MQ=255)
                                    ccGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATTTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAg                                                                                                           <  2:2399833/150‑1 (MQ=255)
                                    ccGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAg                                                                                                           <  1:2001948/150‑1 (MQ=255)
                                    ccGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAg                                                                                                           <  1:2002423/150‑1 (MQ=255)
                                    ccGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAg                                                                                                           <  2:1057603/150‑1 (MQ=255)
                                      ggCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGCCAgtt                                                                                                         >  1:769824/1‑148 (MQ=255)
                                      ggCTGGCCAGGGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCGCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGt                                                                                                         <  2:862330/150‑1 (MQ=255)
                                            ccAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTcttctt                                                                                                   >  2:2380967/1‑150 (MQ=255)
                                             cAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTg                                                                                                  >  1:717890/1‑150 (MQ=255)
                                             cAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcggg                                                                                                                                         >  2:1332668/1‑113 (MQ=255)
                                             cAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcggg                                                                                                                                         <  1:1332668/113‑1 (MQ=255)
                                               gTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGaa                                                                                               <  1:911897/150‑1 (MQ=255)
                                                tGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAAc                                                                                              <  1:1238133/150‑1 (MQ=255)
                                                   tCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCAGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTgg                                                                                           <  2:215245/150‑1 (MQ=255)
                                                     gcgcAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCt                                                                                         <  2:1383297/150‑1 (MQ=255)
                                                       gcAATGCGGGCGTGGTTTTTCCGACCCACTCTACATCTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGCGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGa                                                                                       <  1:1139936/150‑1 (MQ=255)
                                                             cGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGTAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCt                                                                                 <  2:400351/150‑1 (MQ=255)
                                                             cGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCt                                                                                 <  1:1087829/150‑1 (MQ=255)
                                                              gggCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTg                                                                                <  2:2368189/150‑1 (MQ=255)
                                                                gCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGa                                                                              >  1:1028476/1‑150 (MQ=255)
                                                                 cGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                             <  1:2380967/150‑1 (MQ=255)
                                                                     gTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTaa                                                                         <  1:175587/150‑1 (MQ=255)
                                                                      tttttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcgg                                                                                                                                          <  2:235989/87‑1 (MQ=255)
                                                                      tttttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCTGGCAcggcgg                                                                                                                                          >  1:235989/1‑87 (MQ=255)
                                                                          tCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGcc                                                                    <  2:434783/150‑1 (MQ=255)
                                                                            cGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCAc                                                                  <  1:2326256/150‑1 (MQ=255)
                                                                                   ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                           >  2:965665/1‑150 (MQ=255)
                                                                                   ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                           >  1:2322827/1‑150 (MQ=255)
                                                                                        cATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTa                                                      >  2:551648/1‑150 (MQ=255)
                                                                                        cATGTTGCGGCACTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGAGCTCAACTGTATTCTGGGATTTGCAGGAGTGGGGGTGCAGATGCGCAATGCAGGTCTTCTTCAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTa                                                      <  1:1151294/150‑1 (MQ=255)
                                                                                          tGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACg                                                    <  2:1745018/150‑1 (MQ=255)
                                                                                           gTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACgg                                                   <  1:490775/150‑1 (MQ=255)
                                                                                            ttGCGGGATTGGGTCGCTGCTGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGt                                                  <  1:1889031/150‑1 (MQ=255)
                                                                                              gCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGc                                                <  2:717890/150‑1 (MQ=255)
                                                                                              gCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGc                                                <  1:137699/150‑1 (MQ=255)
                                                                                                gggATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTg                                                                                <  2:2173619/116‑1 (MQ=255)
                                                                                                gggATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTg                                                                                >  1:2173619/1‑116 (MQ=255)
                                                                                                  gATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgg                                            >  2:1114128/1‑150 (MQ=255)
                                                                                                    ttGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTTAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                          <  2:1503706/150‑1 (MQ=255)
                                                                                                       ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggacgga                                       >  2:608051/1‑150 (MQ=255)
                                                                                                        gTCGCTGCCGATGTTCAGCGTTGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                      >  1:195393/1‑150 (MQ=255)
                                                                                                        gTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaa                                      >  2:2231303/1‑150 (MQ=255)
                                                                                                         tCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaag                                     >  1:2399283/1‑150 (MQ=255)
                                                                                                           gCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTATGGTGCCCGGACGgaagaa                                   <  1:1145405/150‑1 (MQ=255)
                                                                                                             tGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACa                                 >  2:2040091/1‑150 (MQ=255)
                                                                                                               ccGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGc                               >  2:235417/1‑150 (MQ=255)
                                                                                                               ccGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGc                               <  1:2231303/150‑1 (MQ=255)
                                                                                                                 gATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagtga                                                         >  1:1267446/1‑122 (MQ=255)
                                                                                                                 gATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagtga                                                         <  2:1267446/122‑1 (MQ=255)
                                                                                                                 gATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTg                             <  1:965665/150‑1 (MQ=255)
                                                                                                                    gTTCAGCGGAGCATGGGGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCCGTGACACCCAGTTGCCGGATGTTTTTGCAACTGGCTGAATACCTAGACGTACCCGCCCCTGTCTCGGTGGACGGTCCCCGGACGGAAGAACAGCTGGTc                          <  1:551648/150‑1 (MQ=255)
                                                                                                                     ttCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                             <  2:2129680/114‑1 (MQ=255)
                                                                                                                     ttCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                             >  1:2129680/1‑114 (MQ=255)
                                                                                                                     ttCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTcc                         >  1:2157899/1‑149 (MQ=255)
                                                                                                                      tCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGa                        >  2:1653276/1‑150 (MQ=255)
                                                                                                                      tCAGAGGGGGCTGGCCCTTTACACGCCGTACCAGGGCGGGGAGTCTTCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGa                        <  2:2322827/150‑1 (MQ=255)
                                                                                                                       cAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAg                       <  1:235417/150‑1 (MQ=255)
                                                                                                                       cAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAg                       <  1:2040091/150‑1 (MQ=255)
                                                                                                                         gCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGt                     <  2:819353/150‑1 (MQ=255)
                                                                                                                         gCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTATTCTTGAACTGGCTGAACACCTGGaa                                                                             >  1:393149/1‑94 (MQ=255)
                                                                                                                         gCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTATTCTTGAACTGGCTGAACACCTGGaa                                                                             <  2:393149/94‑1 (MQ=255)
                                                                                                                          cGGGGCCCGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTg                    <  1:1207492/150‑1 (MQ=255)
                                                                                                                                 tGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTgc             >  2:2239711/1‑150 (MQ=255)
                                                                                                                                 tGGCGCTTCACCCGCCGTGCCAGGGAGGGTAGGATGCAGTGGAGATGGACGAGACGGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTgc             >  1:909679/1‑150 (MQ=255)
                                                                                                                                      cTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTg        >  2:1893253/1‑150 (MQ=255)
                                                                                                                                        ctgcccGCCCTGGCACGGCGGGTAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGa      >  1:133790/4‑150 (MQ=255)
                                                                                                                                        ctgcccGCCCTGGCACGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                   <  1:311639/118‑1 (MQ=255)
                                                                                                                                        ctgcccGCCCTGGCACGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                   >  2:311639/4‑121 (MQ=255)
                                                                                                                                         tgcccGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGAt     >  2:2103166/3‑150 (MQ=255)
                                                                                                                                          aCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                             <  2:534464/93‑1 (MQ=255)
                                                                                                                                          aCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                             >  1:534464/1‑93 (MQ=255)
                                                                                                                                          aCCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATg    >  1:2361947/1‑150 (MQ=255)
                                                                                                                                           cccGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                             <  2:594153/76‑1 (MQ=255)
                                                                                                                                           cccGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                             >  1:594153/1‑76 (MQ=255)
                                                                                                                                            cctccgTGCCCGGGCGGGCCGGCTACGGTGCAATTGCTCAAGGCAGGTTTTCTTGAACTGGCTGAACAGCTGTAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGc  <  2:446670/147‑1 (MQ=255)
                                                                                                                                            ccgccgTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTAGGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGc  >  1:446670/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGTTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGGC  >  NC_002947/858177‑858466

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