Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,193,369 A→G 6.4% Q14R (CAG→CGG)  PP_1939 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,193,3690AG 6.4% 180.1 / 2.3 62Q14R (CAG→CGG) PP_1939hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base G (2/2);  total (26/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 = 9.91e-01

GCGAAGGGCCCAGCGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTGGCTTCTTATTA  >  NC_002947/2193229‑2193516
                                                                                                                                             |                                                                                                                                                   
gcgAAGGGCCCAGCGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCa                                                                                                                                             <  2:617/149‑1 (MQ=255)
  gAAGGGCCCAGCGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGc                                                                                                                                           <  7:74379/149‑1 (MQ=255)
      ggCCCAGCGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCg                                                                                                                                       >  6:337454/1‑149 (MQ=255)
             cGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGc                                                                                                                                <  8:56183/149‑1 (MQ=255)
             cGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGc                                                                                                                                <  8:311824/149‑1 (MQ=255)
               ggTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCAc                                                                                                                              <  4:129184/149‑1 (MQ=255)
                  ccATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCAc                                                                                                                              >  4:60269/1‑146 (MQ=255)
                    aTTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCagag                                                                                                                         <  2:207590/149‑1 (MQ=255)
                     ttatAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCGGGTGATCAGCCTCGACCAAGCACCag                                                                                                                           <  8:234188/143‑1 (MQ=255)
                     ttGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCGGGTGATCAGCGTCGACCATGCACCagaag                                                                                                                        >  7:234188/1‑147 (MQ=255)
                     ttGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCGGGTGATCAGCCTCGACCAAGCACCag                                                                                                                           <  8:235949/146‑1 (MQ=255)
                     ttGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAgg                                                                                                                        >  7:93418/1‑149 (MQ=255)
                     ttGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAgg                                                                                                                        >  6:356915/1‑149 (MQ=255)
                     ttGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATAATGTGGGATGGCATCAACATCGCTGAAGTGGCGCGTGTGCGGGTGATCAGCCTCGACCAAGCACCagaag                                                                                                                        >  7:235949/1‑147 (MQ=255)
                       gTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAggg                                                                                                                       <  5:356915/148‑1 (MQ=255)
                         aaaTCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGa                                                                                                                    >  2:84661/1‑149 (MQ=255)
                          aaTCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGTCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCg                                                                                                                     <  1:91644/148‑1 (MQ=255)
                          aaTCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGa                                                                                                                    <  1:317562/148‑1 (MQ=255)
                           aTCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACa                                                                                                                  <  2:332090/149‑1 (MQ=255)
                              gCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGt                                                                                                                <  8:93418/148‑1 (MQ=255)
                              gCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGt                                                                                                                <  3:60269/148‑1 (MQ=255)
                              gCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTg                                                                                                               <  1:84661/149‑1 (MQ=255)
                               cTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGa                                                                                                              <  4:353718/149‑1 (MQ=255)
                               cTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGa                                                                                                              <  7:277749/149‑1 (MQ=255)
                                   ccTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTc                                                                                                          >  1:33264/1‑149 (MQ=255)
                                         gccgccCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGc                                                                                                     >  3:261597/1‑148 (MQ=255)
                                           cgccCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAg                                                                                                  <  4:154956/149‑1 (MQ=255)
                                            gccCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAgg                                                                                                 >  5:120592/1‑149 (MQ=255)
                                               cGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGgcg                                                                                              >  3:186334/1‑149 (MQ=255)
                                                 gCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGcc                                                                                            <  4:261597/149‑1 (MQ=255)
                                                       gggAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATa                                                                                       <  4:135204/148‑1 (MQ=255)
                                                        ggAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATaaa                                                                                     <  2:120469/149‑1 (MQ=255)
                                                         gAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAAt                                                                                    >  4:178309/1‑149 (MQ=255)
                                                         gAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAAt                                                                                    <  6:93127/149‑1 (MQ=255)
                                                                ccTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAAt                                                                             >  7:36988/1‑149 (MQ=255)
                                                                 cTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAAtt                                                                            <  2:4537/149‑1 (MQ=255)
                                                                 cTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAAtt                                                                            <  5:346003/149‑1 (MQ=255)
                                                                      cTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAg                                                                                                             >  8:1198/1‑111 (MQ=255)
                                                                      cTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAg                                                                                                             <  7:1198/111‑1 (MQ=255)
                                                                      cTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGAccc                                                                       <  4:186334/149‑1 (MQ=255)
                                                                        gCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCg                                                                      >  7:30961/1‑148 (MQ=255)
                                                                        gCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCg                                                                      >  4:67370/1‑148 (MQ=255)
                                                                            aCCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCg                                                                                 <  8:106083/133‑1 (MQ=255)
                                                                            aCCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCg                                                                                 >  7:106083/1‑133 (MQ=255)
                                                                              cGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAAt                                                                                        >  7:346341/1‑124 (MQ=255)
                                                                              cGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAAt                                                                                        <  8:346341/124‑1 (MQ=255)
                                                                               gACAAGTTAAATGAGGCGATCATGGGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGCCAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACACGCACAAGAc                                                              >  5:118859/1‑149 (MQ=255)
                                                                                  aaGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGtt                                                           <  5:337454/149‑1 (MQ=255)
                                                                                         aTGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCa                                                    <  3:67370/149‑1 (MQ=255)
                                                                                             ggCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCtt                                                >  3:174609/1‑149 (MQ=255)
                                                                                                gATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTg                                               >  6:68978/1‑147 (MQ=255)
                                                                                                gATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGaa                                             <  3:178309/149‑1 (MQ=255)
                                                                                                     tgtgGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATaa                                        <  5:68978/149‑1 (MQ=255)
                                                                                                           aTCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCaa                                  <  8:30961/149‑1 (MQ=255)
                                                                                                           aTCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCaa                                  <  4:174609/149‑1 (MQ=255)
                                                                                                                aTCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGaa                                             >  2:231074/1‑133 (MQ=255)
                                                                                                                aTCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGaa                                             <  1:231074/133‑1 (MQ=255)
                                                                                                                 tCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAg                            >  4:108247/1‑149 (MQ=255)
                                                                                                                            tGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGccc                  >  3:36380/1‑149 (MQ=255)
                                                                                                                                 tGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTg             >  2:157450/1‑149 (MQ=255)
                                                                                                                                     gcgTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTGGctt         >  8:323669/1‑149 (MQ=255)
                                                                                                                                            cAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTGGCTTCttatta  <  6:120592/149‑1 (MQ=255)
                                                                                                                                             aGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTGGCTTCttatta  >  5:377872/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCGAAGGGCCCAGCGGGTCCATTGTAAATCGCTGCCCTTCAGCCGCCCGGCCTTTGGGAATCGTCCTTACCTGCATACCGACAAGTTAAATGAGGCGATCATGTGGGATCGCATCAACATCGCTGAAGTGGCGCGTGTGCAGGTGATCAGCCTCGACCAAGCACCAGAGGGCGACAGTGAGTTCGACGCTAGGGCGCCGAATAAATTCGTAATTGACCCGCACAAGACGTTCAGCGCAGCTTGAAGATAATCCCAAGACTAGCTGTGCGCCCCGGTGGCTTCTTATTA  >  NC_002947/2193229‑2193516

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: