Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 233,832 A→T 11.7% intergenic (+20/+43) PP_0180 → / ← PP_0181 cytochrome c family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,8320AT11.7% 211.0 / 18.9 77intergenic (+20/+43)PP_0180/PP_0181cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base T (4/5);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

GAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAG  >  NC_002947/233684‑233979
                                                                                                                                                    |                                                                                                                                                   
gAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcac                                                                                                                                                    <  2:1369938/150‑1 (MQ=255)
   ctCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                 >  1:1290550/2‑150 (MQ=255)
     cGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGa                                                                                                                                               <  2:2221255/150‑1 (MQ=255)
       aCTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGAcc                                                                                                                                             <  2:2044472/150‑1 (MQ=255)
           gctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCgcgc                                                                                                                                         >  1:1071107/1‑150 (MQ=255)
             tggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCag                                                                                                                                       <  2:56170/150‑2 (MQ=255)
             tggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCaa                                                                                                                                       >  1:471549/1‑150 (MQ=255)
              ggggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCagg                                                                                                                                      <  2:394381/150‑3 (MQ=255)
               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAgg                                                                                                                                     >  1:2171188/1‑150 (MQ=255)
               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGCACCGCGCAAgg                                                                                                                                     >  2:756664/1‑150 (MQ=255)
                 gATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                   >  2:2015565/1‑150 (MQ=255)
                 gATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                   <  1:1862852/150‑1 (MQ=255)
                   ttCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtctc                                                                                                                                 >  2:2345949/1‑150 (MQ=255)
                      aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                 <  1:2214165/131‑1 (MQ=255)
                      aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                 >  2:2214165/1‑131 (MQ=255)
                       tGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGAt                                                                                                                             <  1:87951/150‑1 (MQ=255)
                           gATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTg                                                                                                                         <  2:471549/150‑1 (MQ=255)
                               ccTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCa                                                                                                                     >  1:585582/1‑150 (MQ=255)
                                 tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                   >  1:1941116/1‑150 (MQ=255)
                                      gCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGTGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGa                                                                                                              >  1:2050326/1‑150 (MQ=255)
                                          tCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCt                                                                                                          >  1:816943/1‑150 (MQ=255)
                                             ccgcATGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTATGGCCGCTCGTAGCTGAAAGAAACGAGCAGTATTGCGACCAACTAACTCCAGCCTGCCCCGTACCACACGGACCGCGCAAGGTGTCAGATCTTGCGCAAGCCTGAtccgcaa                                                                                                       <  2:407984/149‑7 (MQ=255)
                                                 aaGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctg                                                                                                   <  1:2015565/150‑1 (MQ=255)
                                                  atcGGTGGGCTTGCGGGCCATGGTGTTCCTTTACGGCCGCTCGTGCCTGGCCCACAAGCGCCCTGCTGCGGTCAACTGACTGGATCCTGTGCTGTCCCACAGGGTCCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                  <  2:273211/148‑1 (MQ=255)
                                                  aGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                  <  2:2171188/150‑1 (MQ=255)
                                                  aGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                  <  2:585582/150‑1 (MQ=255)
                                                   gCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctggc                                                                                                 <  1:1841596/150‑1 (MQ=255)
                                                      gTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACGGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCtgt                                                                                              <  1:1917010/150‑1 (MQ=255)
                                                      gTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCtgt                                                                                              <  2:1761556/150‑1 (MQ=255)
                                                        ggCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTgg                                                                                            >  2:1027838/1‑150 (MQ=255)
                                                             tGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCt                                                                                       <  2:1605362/150‑1 (MQ=255)
                                                              gCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                      <  2:1071107/150‑1 (MQ=255)
                                                                tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                    <  2:1054865/150‑1 (MQ=255)
                                                                tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                    <  2:816943/150‑1 (MQ=255)
                                                                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGTAGTGCTCCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAAt                                                                                >  2:2080236/1‑150 (MQ=255)
                                                                       cGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                             >  1:1339773/1‑150 (MQ=255)
                                                                       cGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                             <  1:265233/150‑1 (MQ=255)
                                                                         tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        <  2:474488/89‑1 (MQ=38)
                                                                         tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        >  1:474488/1‑89 (MQ=38)
                                                                         tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                  <  2:188856/127‑1 (MQ=255)
                                                                         tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                  >  1:188856/1‑127 (MQ=255)
                                                                            tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGTTTTTCTACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                        <  1:151243/150‑1 (MQ=255)
                                                                            tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                        <  2:2012462/150‑1 (MQ=255)
                                                                            tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                        <  2:1941116/150‑1 (MQ=255)
                                                                               tttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCa                                                                     >  2:2066208/1‑150 (MQ=255)
                                                                                ttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        >  1:1825027/1‑82 (MQ=38)
                                                                                ttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        <  2:1825027/82‑1 (MQ=38)
                                                                                   cGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                 <  2:1339773/150‑1 (MQ=255)
                                                                                     gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        <  1:1955958/77‑1 (MQ=38)
                                                                                     gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                        >  2:1955958/1‑77 (MQ=38)
                                                                                     gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                          <  2:1061517/139‑1 (MQ=255)
                                                                                     gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTg                                                                          >  1:1061517/1‑139 (MQ=255)
                                                                                               gttggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAAc                                                     <  1:271288/148‑1 (MQ=255)
                                                                                                  ggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGc                                                  <  1:730741/150‑1 (MQ=255)
                                                                                                     cGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTc                                               >  1:1178009/1‑150 (MQ=255)
                                                                                                               aGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAgg                                     >  2:1287120/1‑150 (MQ=255)
                                                                                                                     gcGGCCAACTGACTGGAGCCTGTGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTc                               <  2:2050326/150‑1 (MQ=255)
                                                                                                                        gCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAgg                            >  1:1163646/1‑150 (MQ=255)
                                                                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                             <  1:2145314/146‑1 (MQ=255)
                                                                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                             >  2:2145314/1‑146 (MQ=255)
                                                                                                                            actgactgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgc                        <  2:2163811/150‑1 (MQ=255)
                                                                                                                             gtgactgGAGCCGGCGCGGTCCCACAGGGACCTCGCAAGGGCTCAGATCTTGCGCAAGCCTGATGGTCATGCTGGCTGTGGTGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgcc                       <  2:209191/149‑1 (MQ=255)
                                                                                                                                  tgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCAt                  >  2:2052479/1‑150 (MQ=255)
                                                                                                                                    gAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                    <  2:2013901/82‑1 (MQ=255)
                                                                                                                                    gAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                    >  1:2013901/1‑82 (MQ=255)
                                                                                                                                        cTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCt            <  1:2066208/150‑1 (MQ=255)
                                                                                                                                         tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGNGCAAGCNTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg           >  1:2002670/1‑150 (MQ=255)
                                                                                                                                         tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg           >  2:1535005/1‑150 (MQ=255)
                                                                                                                                         tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg           >  1:1517769/1‑150 (MQ=255)
                                                                                                                                           cgcgGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGAt         <  1:1535005/150‑1 (MQ=255)
                                                                                                                                            gcgGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATc        >  1:622677/1‑150 (MQ=255)
                                                                                                                                                tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc    <  1:890944/150‑1 (MQ=255)
                                                                                                                                                 cccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa   <  1:756664/150‑1 (MQ=255)
                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGa                                           <  2:235208/109‑1 (MQ=255)
                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGa                                           >  1:235208/1‑109 (MQ=255)
                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg  >  1:485214/1‑150 (MQ=255)
                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg  >  2:1634518/1‑150 (MQ=255)
                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg  >  2:1847443/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAG  >  NC_002947/233684‑233979

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