Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,349,933 G→A 5.4% L129L (CTG→CTA PP_2945 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,349,9330GA 5.4% 214.6 / 2.2 74L129L (CTG→CTAPP_2945two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base A (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GTTAGCCTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAGCGGCGTA  >  NC_002947/3349789‑3350081
                                                                                                                                                |                                                                                                                                                    
gTTAGCCTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGCCCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGc                                                                                                                                                 <  1:1858995/150‑1 (MQ=255)
gTTAGCCTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGc                                                                                                                                                 <  1:162151/150‑1 (MQ=255)
      cTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCTGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGa                                                                                                                                           >  2:1831872/1‑150 (MQ=255)
      cTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGa                                                                                                                                           >  1:490751/1‑150 (MQ=255)
       tttGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGAt                                                                                                                                          >  1:480296/1‑150 (MQ=255)
         tGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCCGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCg                                                                                                                                        >  2:377968/1‑150 (MQ=255)
             gACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCtt                                                                                                                                     >  1:1097786/1‑149 (MQ=255)
             gACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCtt                                                                                                                                     <  2:1097786/149‑1 (MQ=255)
             gACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTg                                                                                                                                    <  2:1447645/150‑1 (MQ=255)
             gACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTg                                                                                                                                    >  2:1476485/1‑150 (MQ=255)
                aTGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTg                                                                                                                                    >  2:772756/1‑147 (MQ=255)
                aTGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTg                                                                                                                                    <  1:772756/147‑1 (MQ=255)
                aTGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTg                                                                                                                                 >  1:10646/1‑150 (MQ=255)
                aTGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTg                                                                                                                                 >  1:109495/1‑150 (MQ=255)
                 tGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGa                                                                                                                                <  1:1476485/150‑1 (MQ=255)
                           aGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGaaaaa                                                                                                                      >  2:1470927/1‑150 (MQ=255)
                                 atgggatggCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGtt                                                                                                                <  2:1562877/150‑1 (MQ=255)
                                 atgggatggCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGtt                                                                                                                >  1:1137929/1‑150 (MQ=255)
                                     gatggCGCTCAGACCATTGAGCACCTGTGGGACAAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGac                                                                                                            >  1:628170/1‑150 (MQ=255)
                                       tggCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGacaa                                                                                                          >  2:531622/1‑150 (MQ=255)
                                           gcTCAGACCATTGAGCACCTGTGGGGATCGGGCCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATc                                                                                                      >  1:1791199/1‑150 (MQ=255)
                                           gcTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATc                                                                                                      <  1:725027/150‑1 (MQ=255)
                                              cggggcATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAg                                                                                                   <  1:9995/145‑1 (MQ=255)
                                              cAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAg                                                                                                   <  1:1361908/150‑1 (MQ=255)
                                              cAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCAACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAg                                                                                                   <  2:481517/150‑1 (MQ=255)
                                                gACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGt                                                                                                 >  2:802335/1‑150 (MQ=255)
                                                   cATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTAGGATGAGTTGCTCGACCGGCTGCAGGCGCACGGGGGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCa                                                                                              >  1:890949/1‑150 (MQ=255)
                                                    aTTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGCCCTGCTGCAGGCCCAGCATCGCTGGCTGATGCTGAAAAAACCGTTCGACAACATCGAGGTCCAg                                                                                             <  1:1470927/150‑1 (MQ=255)
                                                         gCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAAt                                                                                        >  1:1812817/1‑150 (MQ=255)
                                                           aCCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGtt                                                                                                                >  1:500589/1‑124 (MQ=255)
                                                           aCCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGtt                                                                                                                <  2:500589/124‑1 (MQ=255)
                                                                    aaGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCt                                                                             >  2:1353882/1‑150 (MQ=255)
                                                                     agagGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCtt                                                                            <  1:1394249/150‑1 (MQ=255)
                                                                     agagGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCtt                                                                            >  2:1146200/1‑150 (MQ=255)
                                                                        ggACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGCTGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGGCCAGCAAATGGCCAGCACCTTGCt                                                                         >  1:192590/1‑150 (MQ=255)
                                                                            ccGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTACAGGCCCAc                                                                                                                                             <  2:913404/78‑1 (MQ=255)
                                                                            ccGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTACAGGCCCAc                                                                                                                                             >  1:913404/1‑78 (MQ=255)
                                                                            ccGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTACAGGCCCAc                                                                                                                                             >  1:913549/1‑78 (MQ=255)
                                                                            ccGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTACAGGCCCAc                                                                                                                                             <  2:913549/78‑1 (MQ=255)
                                                                             cGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGt                                                                                                                 <  1:546108/105‑1 (MQ=255)
                                                                             cGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGt                                                                                                                 >  2:546108/1‑105 (MQ=255)
                                                                                    gCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGa                                                             <  2:10646/150‑1 (MQ=255)
                                                                                      aGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCTTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGaaa                                                           >  1:1456640/1‑150 (MQ=255)
                                                                                         tAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGa                                                        <  1:531622/150‑1 (MQ=255)
                                                                                          aGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGAc                                                       <  2:779850/150‑1 (MQ=255)
                                                                                              gtcTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCaa                                                   <  2:480296/150‑1 (MQ=255)
                                                                                                 tGCACTGCCTACTCTGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGa                                                <  1:1831872/150‑1 (MQ=255)
                                                                                                  gCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAg                                               <  1:377968/150‑1 (MQ=255)
                                                                                                     cTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCAc                                                                               <  1:192369/115‑1 (MQ=255)
                                                                                                     cTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCAc                                                                               >  2:192369/1‑115 (MQ=255)
                                                                                                        ccTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGc                                         >  2:652477/1‑150 (MQ=255)
                                                                                                         cTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCt                                        <  1:802335/150‑1 (MQ=255)
                                                                                                            cTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGAc                                                       >  1:1503327/1‑132 (MQ=255)
                                                                                                            cTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGAc                                                       <  2:1503327/132‑1 (MQ=255)
                                                                                                               ggACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAAt                                  >  1:669118/1‑150 (MQ=255)
                                                                                                               ggACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAAt                                  >  2:1863916/1‑150 (MQ=255)
                                                                                                               ggACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAAt                                  <  1:1146200/150‑1 (MQ=255)
                                                                                                                     ttCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCa                            <  2:669118/150‑1 (MQ=255)
                                                                                                                        cTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCt                         >  1:1770642/1‑150 (MQ=255)
                                                                                                                         tGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAg                                                                                                   >  1:1769583/1‑75 (MQ=255)
                                                                                                                         tGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAg                                                                                                   <  2:1769583/75‑1 (MQ=255)
                                                                                                                              tGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACa                   >  2:1599054/1‑150 (MQ=255)
                                                                                                                                 gTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGc                                                                                     >  2:956380/1‑81 (MQ=255)
                                                                                                                                 gTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGc                                                                                     <  1:956380/81‑1 (MQ=255)
                                                                                                                                 gTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCg                                          >  2:1330610/1‑124 (MQ=255)
                                                                                                                                 gTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCg                                          <  1:1330610/124‑1 (MQ=255)
                                                                                                                                      tCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAACCGTTCGACAACATCGCGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCGGCTCGAACAGCGGGgaga          >  2:624706/1‑150 (MQ=255)
                                                                                                                                         aCCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAgc        <  2:192590/150‑1 (MQ=255)
                                                                                                                                          ccGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGAc                                                                      >  2:1465141/1‑87 (MQ=255)
                                                                                                                                          ccGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGAc                                                                      <  1:1465141/87‑1 (MQ=255)
                                                                                                                                            ggCTGAAGGACCAAGACCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAgcggc     <  2:109495/150‑1 (MQ=255)
                                                                                                                                              cTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAgcgg      >  1:735832/1‑147 (MQ=255)
                                                                                                                                              cTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCGTGCGGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAgcgg      <  2:735832/147‑1 (MQ=255)
                                                                                                                                               tGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAGCGGCGTa  <  2:1770642/150‑1 (MQ=255)
                                                                                                                                               cgCAGGCCCACGATCGCTTGCTGATCCTAAAAAAGCCGTTCGATAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAGCGGCGTa  <  2:628170/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTTAGCCTTTGTCGACATGCGAATGCCAGGTGGATGGGATGGCGCTCAGACCATTGAGCACCTGTGGGAAGAGGACCCGCTGTTGCAGGTAGTCGTCTGCACTGCCTACTCGGACTATTCCTGGGATGAGTTGCTCGACCGGCTGCAGGCCCACGATCGCTTGCTGATCCTGAAAAAACCGTTCGACAACATCGAGGTCCAGCAAATGGCCAGCACCTTGCTGACCAAGTGGGAAATGACCCAACGAGCCTCGCTGAAAATGCACCAGCTCGAACAGCGGGTAGAGCGGCGTA  >  NC_002947/3349789‑3350081

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: