Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 873,421 G→A 5.3% G18S (GGC→AGC)  PP_5457 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947873,4210GA 5.3% 277.7 / 3.5 93G18S (GGC→AGC) PP_5457hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base A (3/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569
                                                                                                                                                     |                                                                                                                                                    
ccGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTggcagg                                                                                                                                                     >  2:659298/1‑150 (MQ=255)
 cGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGCCCGACACCTGATGCAGGGAAGCTACCA‑ACCGGAAAGCCTGCTTGAGGCGCAAATTCTggcaggg                                                                                                                                                    >  1:585928/1‑149 (MQ=255)
     gATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGAGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGc                                                                                                                                                >  2:1484821/1‑150 (MQ=255)
       ttGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCACGTCGACGACACCTATTGC‑TAGTCCCCGGTATAGCTGTGAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCtt                                                                                                                                              <  1:976503/150‑1 (MQ=255)
       gtGCCAGGTGTGCGGGGAGCTTGTGGTGTGCATTCTGCAGGTGCATGGTGACGACTGGATGCATGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTTGCAGGCATGCtt                                                                                                                                              <  2:585928/149‑1 (MQ=255)
          ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGcc                                                                                                                                           >  1:244108/1‑150 (MQ=255)
          ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGcc                                                                                                                                           >  1:1704076/1‑150 (MQ=255)
           cAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                          <  2:446797/150‑1 (MQ=255)
               tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTCTCTACGG‑AACGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGGGa                                                                                                                                      <  2:544171/150‑1 (MQ=255)
               tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGa                                                                                                                                      <  1:1830501/150‑1 (MQ=255)
               tgtgCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAAGCATGCCTGCCAGcat                                                                                                                                      <  2:1519805/150‑3 (MQ=255)
                    tGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                 >  2:535819/1‑150 (MQ=255)
                     gTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGAt                                                                                                                                >  2:920705/1‑150 (MQ=255)
                      tCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATc                                                                                                                               <  1:494092/150‑1 (MQ=255)
                         gCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                            >  2:316299/1‑150 (MQ=255)
                          ccTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGAACCTTTATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATg                                                                                                                           <  1:118648/150‑1 (MQ=255)
                          ccTGTGGTGTGCATTCTGCAGGTGCAGGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTTAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATg                                                                                                                           >  2:803597/1‑150 (MQ=255)
                            tgtgGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAATCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTa                                                                                                                         <  2:1312422/150‑1 (MQ=255)
                             gtggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTac                                                                                                                        <  2:665999/150‑1 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACg                                                                                                                     >  1:1146053/1‑150 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACg                                                                                                                     <  1:920705/150‑1 (MQ=255)
                                gtgtGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACg                                                                                                                     >  1:223080/1‑150 (MQ=255)
                                  gtgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTg                                                                                                                   >  1:1642088/1‑150 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGTAATGATGCTTGCCAGCGAGGGGATCGATGTACACGTgg                                                                                                                  >  1:575886/1‑150 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                 >  2:673052/1‑135 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                 <  1:673052/135‑1 (MQ=255)
                                   tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTgg                                                                                                                  >  1:89204/1‑150 (MQ=255)
                                      aTTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCAAGCGAGGGGAGCGATGTACACGTGGTCg                                                                                                               >  2:1545900/1‑150 (MQ=255)
                                              ggTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGAc                                                                                                       <  2:1090164/150‑1 (MQ=255)
                                               gTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGAcc                                                                                                      >  2:1456621/1‑150 (MQ=255)
                                                tGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCt                                                                                                     <  2:1146053/150‑1 (MQ=255)
                                                       tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑AGCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGCGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              >  2:1351356/1‑150 (MQ=255)
                                                       tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              <  1:748510/150‑1 (MQ=255)
                                                       tGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGt                                                                                              <  2:1872158/150‑1 (MQ=255)
                                                           gAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTtccg                                                                                          >  1:736150/1‑150 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGCCACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGa                                                                                                        <  2:925868/126‑1 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGa                                                                                                        >  1:925868/1‑126 (MQ=255)
                                                                     ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCt                                                                                >  2:1375073/1‑150 (MQ=255)
                                                                         aGTCCGCCGCGAGGACGACACCTGATGCAGAGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGccc                                                                            >  2:1221402/1‑150 (MQ=255)
                                                                          gTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGNAAAGNCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGACGTGACCTTGTCGGTGCCGCCGGCGAGCTGCACt                                                                           >  1:211936/1‑150 (MQ=255)
                                                                           tCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCtg                                                                          >  1:1683119/1‑149 (MQ=255)
                                                                              gccgcGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAg                                                                                  <  1:111092/139‑1 (MQ=255)
                                                                              gccgcGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAg                                                                                  >  2:111092/1‑139 (MQ=255)
                                                                               ccgcGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAg                                                                      >  1:31757/1‑150 (MQ=255)
                                                                               ccgcGAGGACGACACCTGATGCAGAGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAg                                                                      <  1:1221402/150‑1 (MQ=255)
                                                                                cgcgAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     <  1:535819/150‑1 (MQ=255)
                                                                                cgcgAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAgg                                                                     <  2:1761818/150‑1 (MQ=255)
                                                                                    aGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCct                                                                 >  1:1117955/1‑150 (MQ=255)
                                                                                     ggACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCctg                                                                >  2:1751402/1‑150 (MQ=255)
                                                                                       acgacACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCctgct                                                              <  2:31757/150‑1 (MQ=255)
                                                                                        cgacACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCctgctg                                                             <  1:855180/150‑1 (MQ=255)
                                                                                         gacACCTGATGCAGCGTATCAACGA‑ACCGGATAGCCCGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgg                                                            <  2:575886/150‑1 (MQ=255)
                                                                                           cacCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGGCCAGCGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTgggg                                                          <  1:18283/150‑1 (MQ=255)
                                                                                            acCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGc                                                         >  2:1653142/1‑150 (MQ=255)
                                                                                            acCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTGGCAGGGCCTGCTGGGGc                                                         >  1:549568/1‑150 (MQ=255)
                                                                                             ccTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTATAAGTGGTCTGCCGTTACCTGGGCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCt                                                        >  1:697752/1‑150 (MQ=255)
                                                                                              cTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCtc                                                       >  2:1627628/1‑150 (MQ=255)
                                                                                                   gCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTcgtcgt                                                  <  1:316299/150‑1 (MQ=255)
                                                                                                   gCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTCCTGGGGCtcgtcct                                                  >  2:205334/1‑148 (MQ=255)
                                                                                                     ggcttagcTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                >  2:505716/8‑150 (MQ=255)
                                                                                                     aGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                >  1:718761/1‑150 (MQ=255)
                                                                                                     aGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                <  1:1653142/150‑1 (MQ=255)
                                                                                                     aGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                <  1:1751402/150‑1 (MQ=255)
                                                                                                     aGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTgg                                                <  2:500211/150‑1 (MQ=255)
                                                                                                         tATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtga                                            <  1:187882/150‑1 (MQ=255)
                                                                                                         tATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCGTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtga                                            <  1:1545900/150‑1 (MQ=255)
                                                                                                           tCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCtgatg                                          >  1:195646/1‑150 (MQ=255)
                                                                                                               cGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCa                                      >  2:1443556/1‑150 (MQ=255)
                                                                                                               cGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCa                                      <  2:89204/150‑1 (MQ=255)
                                                                                                                gA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGc                                                                              >  2:365264/1‑109 (MQ=255)
                                                                                                                gA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGc                                                                              <  1:365264/109‑1 (MQ=255)
                                                                                                                 a‑aCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAgg                                    >  1:1012035/1‑150 (MQ=255)
                                                                                                                 a‑aCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCCCGTCGTGGCTGATGAGCAgg                                    <  1:963649/150‑1 (MQ=255)
                                                                                                                   aCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   <  1:505716/150‑1 (MQ=255)
                                                                                                                    ccGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGTGCGGTGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGcc                                  <  1:416548/150‑1 (MQ=255)
                                                                                                                     cGGAAAGCCTGATTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCg                                 >  2:1419995/1‑150 (MQ=255)
                                                                                                                          cgatctcTTGAGGCGCAATTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCACCGGCTGCGGATCCTCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATgc                            <  2:1005308/144‑1 (MQ=255)
                                                                                                                              tGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCAc                        >  1:6435/1‑150 (MQ=255)
                                                                                                                                 ttGAGGCGCAAATAGCGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATgcgagatcg                     >  1:1005308/1‑144 (MQ=255)
                                                                                                                                  tGAGGCGCAAAGGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCa                    <  2:170375/150‑1 (MQ=255)
                                                                                                                                      gcgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTg                                                                                                   <  1:227767/67‑1 (MQ=255)
                                                                                                                                      gcgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTg                                                                                                   >  2:227767/1‑67 (MQ=255)
                                                                                                                                       cgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGa               <  1:535318/150‑1 (MQ=255)
                                                                                                                                       cgcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGa               <  1:535329/150‑1 (MQ=255)
                                                                                                                                        gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCa                         <  2:1224732/139‑1 (MQ=255)
                                                                                                                                        gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCa                         >  1:1224732/1‑139 (MQ=255)
                                                                                                                                        ccAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGAt              >  1:544830/2‑150 (MQ=255)
                                                                                                                                         cAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGCTGATCAGGCCGAATGCGCACGGCAGCTGATc             <  1:1484821/150‑1 (MQ=255)
                                                                                                                                          cctcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   >  1:311152/5‑127 (MQ=255)
                                                                                                                                          cctcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGc                                   <  2:311152/123‑1 (MQ=255)
                                                                                                                                           ctcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGa           >  2:1525446/4‑150 (MQ=255)
                                                                                                                                              gCTGGCAGGCGTGCTTGCCCGCGCGGGGATCGGAGCGCACGTGGTCGCCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGc        <  2:1451220/150‑1 (MQ=255)
                                                                                                                                              gCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGc        <  2:1475198/150‑1 (MQ=255)
                                                                                                                                              gCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGc        >  2:1327620/1‑150 (MQ=255)
                                                                                                                                                tGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGt      >  2:482568/1‑150 (MQ=255)
                                                                                                                                                 ggcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTa     >  1:1091037/1‑150 (MQ=255)
                                                                                                                                                    aggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACaa  >  2:663578/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGAGGATTGCCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGA‑ACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTGCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAA  >  NC_002947/873273‑873569

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: ‑