Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,169,775 T→A 6.1% W449R (TGG→AGG)  PP_4547 → glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,169,7750TA 6.1% 198.4 / 2.9 66W449R (TGG→AGG) PP_4547glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base A (2/2);  total (33/33)
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.68e-01

CACTGACTGGCTGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAG  >  NC_002947/5169627‑5169924
                                                                                                                                                    |                                                                                                                                                      
cACTGACTGGCTGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTg                                                                                                                                                       <  1:437633/150‑1 (MQ=255)
   tgactgGCTGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTa                                                                                                                                                    >  2:416991/1‑150 (MQ=255)
    gactgGCTGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTAc                                                                                                                                                   <  2:1394824/150‑1 (MQ=255)
        ggatGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTga                                                                                                                                               >  2:574459/4‑149 (MQ=255)
            gACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCa                                                                                                                                           <  2:398930/150‑1 (MQ=255)
               cGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgtc                                                                                                                                        <  2:1296002/150‑3 (MQ=255)
               cGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                        >  1:377407/1‑150 (MQ=255)
                   cTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCCGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGa                                                                                                                                    >  2:15494/1‑150 (MQ=255)
                            cTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTaaa                                                                                                                           >  1:263352/1‑150 (MQ=255)
                             tGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCt                                                                                                                                                 <  2:1359678/127‑1 (MQ=255)
                             tGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCt                                                                                                                                                 >  1:1359678/1‑127 (MQ=255)
                                  cATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACaa                                                                                                                    <  2:931720/150‑1 (MQ=255)
                                  cATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCCGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACaa                                                                                                                    <  1:15494/150‑1 (MQ=255)
                                    ttCCGATTGGGCCCGTGAAGGATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAgg                                                                                                                  >  1:440118/1‑150 (MQ=255)
                                    ttCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAgg                                                                                                                  <  2:624747/150‑1 (MQ=255)
                                     tacGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCGGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGa                                                                                                                 >  1:289615/3‑150 (MQ=255)
                                         aTTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAacac                                                                                                             <  2:263352/150‑1 (MQ=255)
                                                ccGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGaa                                                                                                      >  1:237745/1‑150 (MQ=255)
                                                   tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                     >  1:10135/1‑148 (MQ=255)
                                                   tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                     <  2:10135/148‑1 (MQ=255)
                                                   tGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACgc                                                                                                   >  2:878975/1‑150 (MQ=255)
                                                         aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                          <  2:122000/122‑1 (MQ=255)
                                                         aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAt                                                                                                                          >  1:122000/1‑122 (MQ=255)
                                                         aTTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCAtt                                                                                             <  2:237745/150‑1 (MQ=255)
                                                           tGGGCTGAGACTGGCTGGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGa                                                                                           <  2:658757/150‑1 (MQ=255)
                                                           tGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTAc                                                                                                                                                   <  1:126006/95‑1 (MQ=255)
                                                           tGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTAc                                                                                                                                                   >  2:126006/1‑95 (MQ=255)
                                                             ggCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGaac                                                                                         >  1:942088/1‑150 (MQ=255)
                                                              gCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTAc                                                                                                                                                   <  2:214394/92‑1 (MQ=255)
                                                              gCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTAc                                                                                                                                                   >  1:214394/1‑92 (MQ=255)
                                                                gAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGaacaac                                                                                      <  1:1039064/150‑1 (MQ=255)
                                                                         cTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                        >  2:1110437/1‑92 (MQ=255)
                                                                         cTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGc                                                                                                                                        <  1:1110437/92‑1 (MQ=255)
                                                                               gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgt                                                                                                                                         <  1:884188/85‑2 (MQ=37)
                                                                               gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGCCTGGTGCAGGTACCAGCGCCAgt                                                                                                                                         >  2:884188/1‑84 (MQ=25)
                                                                               gtgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAgcgc                                                                       >  2:1337845/1‑150 (MQ=255)
                                                                                tgtACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAgcgcg                                                                      >  2:1178279/1‑150 (MQ=255)
                                                                                     cTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGcc                                                                                                                                       <  2:529483/81‑1 (MQ=255)
                                                                                     cTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGcc                                                                                                                                       >  1:529483/1‑81 (MQ=255)
                                                                                      tGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGt                                                                <  2:440118/150‑1 (MQ=255)
                                                                                      tGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCAGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAATAAACTACG‑ACATCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGt                                                                >  1:1156085/1‑150 (MQ=255)
                                                                                       gAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTc                                                               <  1:1178279/150‑1 (MQ=255)
                                                                                           gTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCctac                                                           <  1:416991/150‑1 (MQ=255)
                                                                                               aGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                     >  2:884152/1‑104 (MQ=255)
                                                                                               aGCGCGCCGAGTACCGCCAGTTCATGGCAGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAAc                                                                                                     <  1:884152/104‑1 (MQ=255)
                                                                                                 cgcgcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAgc                                                     <  2:377407/150‑1 (MQ=255)
                                                                                                    gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGcc                                                  >  1:1197863/1‑150 (MQ=255)
                                                                                                    gcCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGcc                                                  >  1:1197156/1‑150 (MQ=255)
                                                                                                             cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                         >  1:606344/1‑150 (MQ=255)
                                                                                                             cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                         >  2:876118/1‑150 (MQ=255)
                                                                                                             cGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAAc                                         <  1:981082/150‑1 (MQ=255)
                                                                                                              gCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACg                                        >  2:887143/1‑150 (MQ=255)
                                                                                                              gCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACg                                        >  1:1075286/1‑150 (MQ=255)
                                                                                                               ccAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGa                                       >  1:782670/1‑150 (MQ=255)
                                                                                                               ccAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCAACAGGCCTGAACCCCTAAAT‑AAGGAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGa                                       <  1:878975/150‑1 (MQ=255)
                                                                                                                cAGCTCATGGCCGAAGGCAGCGAACAGCACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACCCCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGAc                                      <  2:145404/150‑1 (MQ=255)
                                                                                                                  gTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCAGGCGCAGCGCGCGTCGTCCTACTACAGCGCCACCCTCCACGACca                                    <  1:1176574/150‑1 (MQ=255)
                                                                                                                        ggCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGc                              <  2:956658/150‑1 (MQ=255)
                                                                                                                        ggCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGc                              <  2:782670/150‑1 (MQ=255)
                                                                                                                           cGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTa                           <  2:289615/150‑1 (MQ=255)
                                                                                                                                   gCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCa                                           <  2:505230/126‑1 (MQ=255)
                                                                                                                                   gCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCa                                           >  1:505230/1‑126 (MQ=255)
                                                                                                                                    cGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACAGCGTACCCGCAACt                  >  2:369229/1‑150 (MQ=255)
                                                                                                                                     gAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTc                 >  2:1326321/1‑150 (MQ=255)
                                                                                                                                         aGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCAACCTCAACGACCCCCCTGCGTACCCGCAACTCAAgg             <  2:1156085/150‑1 (MQ=255)
                                                                                                                                           gccTGGTGCAGGTACCAGCGCCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGc           >  1:1353054/3‑150 (MQ=255)
                                                                                                                                                    tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg  <  1:1326321/150‑1 (MQ=255)
                                                                                                                                                    tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg  >  1:64770/1‑150 (MQ=255)
                                                                                                                                                    tGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAg  <  1:887143/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CACTGACTGGCTGACCGCACTCGACGCCCTGCAACATTCCGATTGGGCCCGTGAAGCATTGGGCGAAGACTTCCTTGGCGTGTACCTGAAGGTCAAGCGCGCCGAGTACCGCCAGTTCATGGCCGAAGTCAGCGAACAGGACTGGCGCTGGTACCTGCACCAGGCCTGAGCCCCTAAAT‑AACAAGGAACACCCATGAACGCAGCATTGAACAACGGCCCGGCGCAGCGCGCGCCGTCCTACTACAGCGCCACCCTCAACGACCACACTGCGTACCCGCAACTCAAGGGCACGGTGCAG  >  NC_002947/5169627‑5169924

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