Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,386,384 A→C 8.3% R240R (CGT→CGG)  PP_2094 ← nitrate‑binding protein NasS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,386,3840AC 8.3% 173.8 / 7.4 60R240R (CGT→CGG) PP_2094nitrate‑binding protein NasS
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base C (3/2);  total (31/29)
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.82e-01

CTCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAACGCAGA  >  NC_002947/2386237‑2386532
                                                                                                                                                   |                                                                                                                                                    
cTCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCAc                                                                                                                                                     <  6:133895/149‑1 (MQ=255)
cTCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCAc                                                                                                                                                     <  6:234494/149‑1 (MQ=255)
 taaTAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    >  5:93150/3‑149 (MQ=255)
 tCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    >  2:356787/1‑149 (MQ=255)
 tCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    >  2:147907/1‑149 (MQ=255)
 tCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    >  3:215533/1‑149 (MQ=255)
  cATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACgg                                                                                                                                                   >  5:52304/1‑149 (MQ=255)
   atatTGGGCGTGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGGTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTAGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    <  7:91628/147‑1 (MQ=255)
   aTAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACgg                                                                                                                                                   <  1:53788/148‑1 (MQ=255)
   aTAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACg                                                                                                                                                    <  5:364619/147‑1 (MQ=255)
         gCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGt                                                                                                                                            >  2:291539/1‑149 (MQ=255)
           cGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTg                                                                                                                                          <  2:302043/149‑1 (MQ=255)
            gAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTgg                                                                                                                                         <  1:356787/149‑1 (MQ=255)
            gAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTgg                                                                                                                                         <  2:315884/149‑1 (MQ=255)
              gCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTgggg                                                                                                                                       <  2:201869/149‑1 (MQ=255)
              gCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTGGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCCCGTGCGCTgatcaa                                                                                                                                        <  8:340442/148‑6 (MQ=255)
               cAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCCCGTGCGCTgatcaagg                                                                                                                                      <  5:325902/149‑8 (MQ=255)
               cAGGCGGGGGGCTATGCTGGCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGGTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGAGCGCCTTGATCAGCGCACGGGCGCTGTTGGGGt                                                                                                                                      <  1:313174/149‑1 (MQ=255)
                 ggCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGGTGTTGGGGTAt                                                                                                                                    <  6:248886/149‑1 (MQ=255)
                         gCTCTGCAGTCCCGGGGGGTGGCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGc                                                                                                                            <  6:259498/149‑1 (MQ=255)
                          cTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGc                                                                                                                            >  3:298179/1‑148 (MQ=255)
                                  tCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGAcgcg                                                                                                                   <  1:147907/149‑1 (MQ=255)
                                   ccAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCgg                                                                                                                  <  4:215533/149‑1 (MQ=255)
                                    cAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCgg                                                                                                                  <  4:298179/148‑1 (MQ=255)
                                         gggTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAg                                                                                                            >  3:320253/1‑149 (MQ=255)
                                             gtCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGc                                                                                                          >  6:128850/1‑147 (MQ=255)
                                               ccAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAg                                                                                                       >  8:170329/1‑148 (MQ=255)
                                                cAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGa                                                                                                     <  6:93150/149‑1 (MQ=255)
                                                 aGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGAc                                                                                                    >  8:71756/1‑149 (MQ=255)
                                                   gTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCt                                                                                                  >  2:29806/1‑149 (MQ=255)
                                                    tAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCtt                                                                                                 >  4:314240/1‑149 (MQ=255)
                                                        cgcTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTcc                                                                                             <  7:71756/149‑1 (MQ=255)
                                                          cTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCgg                                                                                           <  7:170329/149‑1 (MQ=255)
                                                             gcgcTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCggg                                                                                          <  3:314240/147‑1 (MQ=255)
                                                               gcTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGAt                                                                                      >  2:181603/1‑149 (MQ=255)
                                                                 tGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCtt                                                                                                 >  7:154943/1‑136 (MQ=255)
                                                                 tGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCtt                                                                                                 <  8:154943/136‑1 (MQ=255)
                                                                 tGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATcc                                                                                    >  3:189353/1‑149 (MQ=255)
                                                                      aGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGc                                                                                 >  3:224609/1‑147 (MQ=255)
                                                                                     cggcggTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCtggt                                                                >  2:47118/1‑149 (MQ=255)
                                                                                             ttCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAgg                                                         <  1:291539/148‑1 (MQ=255)
                                                                                                    cTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGcc                                                 <  4:189353/149‑1 (MQ=255)
                                                                                                       tGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTg                                              >  4:222230/1‑149 (MQ=255)
                                                                                                         cTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGc                                            >  8:51257/1‑149 (MQ=255)
                                                                                                          tCGATGAATCGGCTGGCGGTGAGTATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGcc                                           >  5:180234/1‑149 (MQ=255)
                                                                                                             aTGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCttt                                        >  1:238314/1‑149 (MQ=255)
                                                                                                                 aTCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGCCTCCTGGCCAGGGTGAAGCCTTGGCCTTTgga                                     >  5:307105/1‑147 (MQ=255)
                                                                                                                    ggCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCAcc                                 >  1:201073/1‑149 (MQ=255)
                                                                                                                       tGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACcgcg                              >  2:301710/1‑149 (MQ=255)
                                                                                                                                 gATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGcc                     >  3:87368/1‑148 (MQ=255)
                                                                                                                                        cccaacAGCGCCCGTGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCt              >  1:379291/6‑148 (MQ=255)
                                                                                                                                        cccaacAGCGCCCGTGCGCTGATGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCt              >  7:279323/6‑148 (MQ=255)
                                                                                                                                         ccaacAGCGCCCGTGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCa            >  5:307500/5‑149 (MQ=255)
                                                                                                                                          gATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCa            <  2:238314/148‑1 (MQ=255)
                                                                                                                                          gATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCa            <  2:371270/148‑1 (MQ=255)
                                                                                                                                           aTCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGAGCGCCTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCAc           <  4:224609/148‑1 (MQ=255)
                                                                                                                                             cAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCaa        >  5:87591/1‑149 (MQ=255)
                                                                                                                                              aGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAAc       <  4:320253/149‑1 (MQ=255)
                                                                                                                                                 gcACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAACGCa    <  2:201073/149‑1 (MQ=255)
                                                                                                                                                  cACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAACg      >  8:37768/1‑146 (MQ=255)
                                                                                                                                                   aCGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAACGCAGa  <  1:47118/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTCATAATGGCCGAGCAGGCGGGGGGCTATGCTGTCCAGCGGGGTGTCCAGGTAAGCGCTTGCGCTGAGCAGTTGCGCGGTGCTGCGGCGGTTTTCCGGGCTTTGCTCGATGAATCGGCTGGCGGTGAGGATCGCCTTGATCAGCGCACGGGCGCTGTTGGGGTATTGCTCGGCAAAGACGCGGGCGCAGGCCAGGACCTTTTCCGGGTGATCCGGCCAGATCGACTGGCTGGTGGCCAGGGTGAAGCCTTGGCCTTTGGCCACCGCGTCGGCGGCCCAGGGCTCACCAACGCAGA  >  NC_002947/2386237‑2386532

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