Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,882,554 A→G 5.7% *103W (TAG→TGG)  PP_2537 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,882,5540AG 5.7% 217.4 / 2.3 70*103W (TAG→TGG) PP_2537hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base G (2/2);  total (33/37)
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.73e-01

CCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTGGCCGC  >  NC_002947/2882407‑2882702
                                                                                                                                                   |                                                                                                                                                     
ccGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGc                                                                                                                                                     >  2:1570972/1‑150 (MQ=255)
        aCCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTGGCTAGGCttc                                                                                                                                             <  1:1519477/150‑2 (MQ=255)
          ccGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGCCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGcc                                                                                                                                          >  2:568575/1‑150 (MQ=255)
              aaGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTggg                                                                                                                                      <  2:547178/150‑1 (MQ=255)
                ggCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCAAGCCTggcta                                                                                                                                    <  1:1476363/150‑4 (MQ=255)
                       gACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACCCTCAGCTGCTGCAGCTACGTCTCGCGCGCTAAGTCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGaa                                                                                                                             <  1:200128/150‑1 (MQ=255)
                          gggCCTGACGAACGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACtt                                                                                                                          >  2:808570/2‑150 (MQ=255)
                          aGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACtt                                                                                                                          <  1:1570972/150‑1 (MQ=255)
                           ggCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTg                                                                                                                         <  1:1580167/150‑1 (MQ=255)
                                    acCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTg                                                                                                                <  1:1289123/150‑1 (MQ=255)
                                         cAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCgg                                                                                                           <  2:938031/150‑1 (MQ=255)
                                          aGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGt                                                                                                          >  1:167044/1‑150 (MQ=255)
                                           ggTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGtt                                                                                                         <  2:1138082/150‑1 (MQ=255)
                                           ggTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGCCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGtt                                                                                                         <  1:568575/150‑1 (MQ=255)
                                            gTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCt                                                                                                                  >  2:35721/1‑140 (MQ=255)
                                            gTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCt                                                                                                                  <  1:35721/140‑1 (MQ=255)
                                               atgatgaAGGCCACGTACTACATGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTg                                                                                                                         <  2:693685/130‑1 (MQ=255)
                                               atgatgaAGGCCACGTACTACATGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTg                                                                                                                         >  1:693685/1‑130 (MQ=255)
                                                    gaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTc                                                                                                >  2:1362073/1‑150 (MQ=255)
                                                     ccggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCt                                                                                               >  1:1286289/3‑150 (MQ=255)
                                                     aaGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCt                                                                                               <  2:1565121/150‑1 (MQ=255)
                                                         ccACGTACTACCTGGCCGGGATGGCAAAGTCACCGCTCGACGGCCCGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTaa                                                                                           <  2:1503479/150‑1 (MQ=255)
                                                          cACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAac                                                                                          <  1:1362073/150‑1 (MQ=255)
                                                           aCGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTggg                                                                                                                                      <  2:720332/105‑1 (MQ=255)
                                                           aCGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTggg                                                                                                                                      >  1:720332/1‑105 (MQ=255)
                                                                 tacCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAg                                                                                   >  1:614431/1‑150 (MQ=255)
                                                                    cTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCgg                                                                                                           <  1:239696/123‑1 (MQ=255)
                                                                    cTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCgg                                                                                                           >  2:239696/1‑123 (MQ=255)
                                                                    cTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAg                                                                                >  2:1188942/1‑150 (MQ=255)
                                                                        ccGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGaaa                                                                            <  2:1286289/150‑1 (MQ=255)
                                                                            gATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATTCGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATg                                                                        <  1:1500133/150‑1 (MQ=255)
                                                                              tGGCCGAGTAACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCa                                                                      <  1:34723/150‑1 (MQ=255)
                                                                                  cAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGa                                                                  >  1:1278873/1‑150 (MQ=255)
                                                                                      tCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTg                                                              <  2:1278873/150‑1 (MQ=255)
                                                                                               gACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTGTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGt                                                     >  2:334095/1‑150 (MQ=255)
                                                                                                aCGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCAtt                                                                    >  1:1645477/1‑134 (MQ=255)
                                                                                                aCGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCAtt                                                                    <  2:1645477/134‑1 (MQ=255)
                                                                                                     gacctgctgCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGaa                                               >  1:56414/4‑150 (MQ=255)
                                                                                                     cAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACtt                                                                                                                          <  2:1633809/75‑1 (MQ=255)
                                                                                                     cAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACtt                                                                                                                          >  1:1633809/1‑75 (MQ=255)
                                                                                                         tgctgcAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCtt                                           <  2:614431/150‑1 (MQ=255)
                                                                                                           ctgcAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGtt                                                                                                         >  2:1638189/1‑86 (MQ=255)
                                                                                                           ctgcAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGtt                                                                                                         <  1:1638189/86‑1 (MQ=255)
                                                                                                             gcAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACAt                                                                         <  2:1338270/116‑1 (MQ=255)
                                                                                                             gcAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACAt                                                                         >  1:1338270/1‑116 (MQ=255)
                                                                                                                 aTACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGc                                   <  1:808570/150‑1 (MQ=255)
                                                                                                                 aTACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGc                                   >  2:335577/1‑150 (MQ=255)
                                                                                                                     gTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTg                               >  1:404315/1‑150 (MQ=255)
                                                                                                                      tCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGc                              >  1:1297097/1‑150 (MQ=255)
                                                                                                                                  aGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCtaa                  <  2:1297097/150‑3 (MQ=255)
                                                                                                                                  aGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGa                  <  2:404315/150‑1 (MQ=255)
                                                                                                                                    gCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATg                <  2:413046/150‑1 (MQ=255)
                                                                                                                                       ggggCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTc             >  2:1362561/1‑150 (MQ=255)
                                                                                                                                        gggCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTcc            >  1:961168/1‑150 (MQ=255)
                                                                                                                                           gcaaGCCTGGCTAGGCTTCGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTg                                                    <  1:1391949/106‑1 (MQ=255)
                                                                                                                                           gcaaGCCTGGCTAGGCTTCGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTg                                                    >  2:1391949/3‑108 (MQ=255)
                                                                                                                                            gAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCtcgt                                      >  2:161072/1‑120 (MQ=255)
                                                                                                                                            gAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCtcgt                                      <  1:161072/120‑1 (MQ=255)
                                                                                                                                              aGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTgg      <  2:167044/150‑1 (MQ=255)
                                                                                                                                               gCCTGGCTAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTGgc     >  2:1642456/1‑150 (MQ=255)
                                                                                                                                               gCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTaa                                                                                                     <  1:114107/54‑1 (MQ=255)
                                                                                                                                               gCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTaa                                                                                                     >  2:114107/1‑54 (MQ=255)
                                                                                                                                               gCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAaca                                                                                         <  1:742829/66‑1 (MQ=255)
                                                                                                                                               gCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAaca                                                                                         >  2:742829/1‑66 (MQ=255)
                                                                                                                                                 cTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGaa                                               <  2:931261/106‑1 (MQ=255)
                                                                                                                                                 cTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGaa                                               >  1:931261/1‑106 (MQ=255)
                                                                                                                                                 cTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCtt                                           >  1:156817/1‑110 (MQ=255)
                                                                                                                                                 cTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCtt                                           <  2:156817/110‑1 (MQ=255)
                                                                                                                                                 cTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTGgccg   >  1:384874/1‑150 (MQ=255)
                                                                                                                                                  tAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTGgccgc  >  2:756278/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGACAGAAGGAAACATGCATTGAGCTGGGACAAGGTGGGGAAGCTTCTCGTCGCCGTGCTGATCGGCCTGATGGTCCTGTTGGCCGC  >  NC_002947/2882407‑2882702

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: