Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,849 A→T 12.3% K120* (AAG→TAG) ‡ PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,8490AT12.3% 208.4 / 16.4 73K120* (AAG→TAG) ‡PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (5/4);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CAGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTC  >  NC_002947/1899701‑1899996
                                                                                                                                                    |                                                                                                                                                   
cagccagcCGCCAGTGCGCACCAGCATTGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCcaacaa                                                                                                                                                    >  1:393419/1‑150 (MQ=255)
        cgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCa                                                                                                                                            <  1:1098323/150‑1 (MQ=255)
         ggcAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCTGAGCGTCTGCGCGGCATCATCGACCATCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCag                                                                                                                                           <  2:34363/148‑2 (MQ=255)
          ccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAg                                                                                                                                          >  2:326732/1‑150 (MQ=255)
          ccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAg                                                                                                                                          <  1:1063816/150‑1 (MQ=255)
              tGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCt                                                                                                                                      <  2:1317050/150‑1 (MQ=255)
                   aCCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATc                                                                                                                                 >  2:1001426/1‑150 (MQ=255)
                     cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCg                                                                                                                               >  1:876242/1‑150 (MQ=255)
                     cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCg                                                                                                                               >  2:674514/1‑150 (MQ=255)
                      aGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGc                                                                                                                                              >  2:507720/1‑134 (MQ=255)
                      aGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGc                                                                                                                                              <  1:507720/134‑1 (MQ=255)
                             aCCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCAATTCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGcc                                                                                                                       >  1:1170071/1‑150 (MQ=255)
                                  gCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                  <  2:1273547/150‑1 (MQ=255)
                                    cATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCTTTGCGCTGTCCGTGTCCGTGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGa                                                                                                                <  2:1404970/150‑1 (MQ=255)
                                     aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGtt                                                                                                                                                 >  1:695500/1‑116 (MQ=255)
                                     aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGtt                                                                                                                                                 <  2:695500/116‑1 (MQ=255)
                                         cTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCa                                                                                                           >  1:327865/1‑150 (MQ=255)
                                              aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcg                                                                                                      >  2:1176374/1‑150 (MQ=255)
                                                       cGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGc                                                                                             >  1:205187/1‑150 (MQ=255)
                                                       cGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGc                                                                                             >  1:359922/1‑150 (MQ=255)
                                                         tCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCAt                                                                                           >  2:354425/1‑150 (MQ=255)
                                                          cTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                          <  1:944644/150‑1 (MQ=255)
                                                            gcgcGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgc                                                                                        <  1:571050/150‑1 (MQ=255)
                                                                 gCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcg                                                                                   >  2:734512/1‑150 (MQ=255)
                                                                      atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAg                                                                                                                                          <  2:589714/90‑1 (MQ=255)
                                                                      atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAg                                                                                                                                          >  1:589714/1‑90 (MQ=255)
                                                                      atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCACGCCTGTGCATGtccagccccagc                                                                              >  2:267575/1‑150 (MQ=255)
                                                                       tcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                             >  2:1083508/1‑118 (MQ=255)
                                                                       tcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                             <  1:1083508/118‑1 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                           <  1:354425/150‑1 (MQ=255)
                                                                           ccAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                         <  2:359922/150‑1 (MQ=255)
                                                                            cAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAgg                                                                                                                                         >  2:822940/1‑85 (MQ=255)
                                                                            cAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAgg                                                                                                                                         <  1:822940/85‑1 (MQ=255)
                                                                              gCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACATGGTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAATCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTa                                                                      >  1:99659/1‑150 (MQ=255)
                                                                              gCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTa                                                                      >  2:705221/1‑150 (MQ=255)
                                                                                 cGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACaa                                                                   <  1:326732/150‑1 (MQ=255)
                                                                                        gATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCg                                                            <  2:1327358/150‑1 (MQ=255)
                                                                                         aTGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCg                                                                                                                               >  2:1296172/1‑82 (MQ=255)
                                                                                         aTGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCg                                                                                                                               <  1:1296172/82‑1 (MQ=255)
                                                                                         aTGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGc                                                           <  1:105511/150‑1 (MQ=255)
                                                                                            ccGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCCGTCCCTGGTCAACAGGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                        <  1:674514/150‑1 (MQ=255)
                                                                                            ccGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                        <  1:620982/150‑1 (MQ=255)
                                                                                            ccGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCGACAAGCTGTCGCTTg                                                        <  1:1000854/150‑1 (MQ=255)
                                                                                              gTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGgc                                                      >  2:1275963/1‑150 (MQ=255)
                                                                                                      gTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                              <  1:1001426/150‑1 (MQ=255)
                                                                                                      gTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                              <  2:393419/150‑1 (MQ=255)
                                                                                                       tACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTTACTTGGCGCCGATGGc                                             <  2:710008/150‑1 (MQ=255)
                                                                                                       tACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGc                                             >  1:791207/1‑150 (MQ=255)
                                                                                                          gtgtCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCa                                          >  1:1146779/1‑150 (MQ=255)
                                                                                                             tCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                       >  2:724425/1‑150 (MQ=255)
                                                                                                             tCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                       >  2:189215/1‑150 (MQ=255)
                                                                                                              ccGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                             >  2:1083174/1‑79 (MQ=255)
                                                                                                              ccGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                             <  1:1083174/79‑1 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                  <  2:905601/64‑1 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                  >  1:905601/1‑64 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                            <  2:1146779/150‑1 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                            <  2:35602/150‑1 (MQ=255)
                                                                                                                         tcacTGGATTCCCATTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAg                           <  2:791207/147‑1 (MQ=255)
                                                                                                                             tGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACt                       >  2:23504/1‑150 (MQ=255)
                                                                                                                              ggAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACtt                      <  2:327865/150‑1 (MQ=255)
                                                                                                                                 ggCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGa                   <  1:1176374/150‑1 (MQ=255)
                                                                                                                                   cccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGc                 <  2:876242/150‑1 (MQ=255)
                                                                                                                                     cAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACaa                                                                   <  1:642472/98‑1 (MQ=255)
                                                                                                                                     cAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACaa                                                                   >  2:642472/1‑98 (MQ=255)
                                                                                                                                         ccttGGCCAACTTGTTGGCCAGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                              >  2:1109169/4‑115 (MQ=255)
                                                                                                                                         ccttGGCCAACTTGTTGGCCAGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATgg                                              <  1:1109169/112‑1 (MQ=255)
                                                                                                                                         ccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGa                         <  1:547635/133‑1 (MQ=255)
                                                                                                                                         ccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGa                         >  2:547635/4‑136 (MQ=255)
                                                                                                                                         ccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCg           >  2:764077/4‑150 (MQ=255)
                                                                                                                                           ttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc         >  2:43015/2‑150 (MQ=255)
                                                                                                                                           ttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc         >  1:922419/2‑150 (MQ=255)
                                                                                                                                           cTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCCCCTGGCAGTTACCAGGACTTCGAGCACGCCGAc         >  2:496778/1‑150 (MQ=255)
                                                                                                                                             ggCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                      <  2:198082/119‑1 (MQ=255)
                                                                                                                                             ggCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                      >  1:198082/1‑119 (MQ=255)
                                                                                                                                                caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt    >  1:384840/1‑150 (MQ=255)
                                                                                                                                                caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt    >  1:444017/1‑150 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAg                           >  1:1353047/1‑126 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAg                           <  2:1353047/126‑1 (MQ=255)
                                                                                                                                                  acaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTc  >  1:1016867/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTC  >  NC_002947/1899701‑1899996

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