Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,431,272 C→A 6.3% Q348K (CAG→AAG)  PP_3925 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,431,2720CA 6.3% 253.9 / 5.9 79Q348K (CAG→AAG) PP_3925hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/33);  new base A (2/3);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

AACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4431125‑4431418
                                                                                                                                                   |                                                                                                                                                  
aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                    <  4:311178/148‑1 (MQ=255)
aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                    >  3:311178/1‑148 (MQ=255)
 aCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAg                                                                                                                                                  <  8:194560/149‑1 (MQ=255)
   ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGt                                                                                                                                                >  8:220993/1‑149 (MQ=255)
      gATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCa                                                                                                                                              <  3:170372/148‑1 (MQ=255)
       aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtt                                                                                                                                            <  1:243805/149‑1 (MQ=255)
        ttACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                           <  8:17066/149‑1 (MQ=255)
            tactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                                       <  8:61754/149‑5 (MQ=255)
               taTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                    >  7:148974/1‑149 (MQ=255)
                 ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                    >  1:333626/1‑147 (MQ=255)
                 ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCa                                                                                                                                  <  8:271931/149‑1 (MQ=255)
                       tGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                            <  1:261994/149‑1 (MQ=255)
                         gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTg                                                                                                                           <  1:13060/148‑1 (MQ=255)
                           cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAg                                                                                                                         >  6:334167/1‑148 (MQ=255)
                            tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                                   >  1:48022/1‑137 (MQ=255)
                            tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                                   <  2:48022/137‑1 (MQ=255)
                                 cGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                                  >  2:250163/1‑149 (MQ=255)
                                  gCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                                  >  4:258563/1‑148 (MQ=255)
                                     tGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCac                                                                                                               >  8:249882/1‑148 (MQ=255)
                                       gcaggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacg                                                                                                              <  4:203202/147‑1 (MQ=255)
                                          ggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgac                                                                                                         <  8:106776/149‑1 (MQ=255)
                                           gcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                        <  1:256759/149‑1 (MQ=255)
                                            caTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                        >  2:205572/1‑148 (MQ=255)
                                                gTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                                   <  5:238678/149‑1 (MQ=255)
                                                 ttCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                                   >  2:159507/1‑148 (MQ=255)
                                                   cagcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACa                                                                                                 <  7:220993/148‑1 (MQ=255)
                                                     gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                                 <  3:119931/98‑1 (MQ=255)
                                                     gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                                 >  4:119931/1‑98 (MQ=255)
                                                      cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAgg                                                                                              >  7:95287/1‑148 (MQ=255)
                                                       agATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCTTTTCGCACAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCAAGCATCTGAGCACAGCCCCCGACGACTGATAACAAGGa                                                                                             >  3:237004/1‑148 (MQ=255)
                                                            tactacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                                       >  8:296169/1‑149 (MQ=255)
                                                              ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGc                                                                                      <  8:95287/148‑1 (MQ=255)
                                                                 ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         >  4:90965/1‑142 (MQ=255)
                                                                 ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                         <  3:90965/142‑1 (MQ=255)
                                                                  cGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTc                                                                                 <  8:148974/149‑1 (MQ=255)
                                                                       ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGAGCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                            <  7:105934/149‑1 (MQ=255)
                                                                       ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                            <  1:250163/149‑1 (MQ=255)
                                                                        tacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAACCAGCATGTCCAAGAc                                                                           >  4:29252/1‑149 (MQ=255)
                                                                          ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     <  8:212156/121‑1 (MQ=255)
                                                                          ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGAcc                                                                          >  6:21001/1‑148 (MQ=255)
                                                                          ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCAGTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                     >  7:212156/1‑121 (MQ=255)
                                                                            aTGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCAc                                                                       <  5:71965/149‑1 (MQ=255)
                                                                                ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                                <  1:207631/88‑1 (MQ=255)
                                                                                ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                                >  2:207631/1‑88 (MQ=255)
                                                                                ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGt                                                                   >  8:305583/1‑149 (MQ=255)
                                                                                 cAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGtt                                                                  <  5:21001/149‑1 (MQ=255)
                                                                                     cccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  2:268988/1‑149 (MQ=255)
                                                                                      ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAcc                                                             >  8:98853/1‑149 (MQ=255)
                                                                                      ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  5:228918/1‑148 (MQ=255)
                                                                                      ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  8:94097/1‑148 (MQ=255)
                                                                                        aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTc                                                           >  3:247072/1‑149 (MQ=255)
                                                                                          gACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         <  1:205572/149‑1 (MQ=255)
                                                                                              aGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     <  1:268988/149‑1 (MQ=255)
                                                                                               gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGCCGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                     >  6:182000/1‑148 (MQ=255)
                                                                                               gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              <  4:161219/139‑1 (MQ=255)
                                                                                               gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                              >  3:161219/1‑139 (MQ=255)
                                                                                                 gTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCaa                                                   <  6:228918/148‑1 (MQ=255)
                                                                                                       cccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGcc                                            >  6:320147/1‑149 (MQ=255)
                                                                                                        ccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCt                                           >  3:217661/1‑149 (MQ=255)
                                                                                                          aGTTCATTCACATGAACCGCTTGCTTGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCttt                                         <  5:334167/149‑1 (MQ=255)
                                                                                                               attcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTc                                    <  1:159507/149‑1 (MQ=255)
                                                                                                                ttgacaTGAACGGCTTGCTGGAACAATTGACCTTGCAGGTCATTTCCAAGTCCTGCCTCTGAGCACAGCCCAGGACGACGGATAACAAGGAAATGAGCATGTGCAAGACGCTCCTGTTCGTCCTCGACGGGAAGATTGCCTTTGTTTcc                                   <  6:187312/146‑1 (MQ=255)
                                                                                                                     aTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgcgc                              >  3:134793/1‑149 (MQ=255)
                                                                                                                        aaCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCa                            <  7:249882/148‑1 (MQ=255)
                                                                                                                         aCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGc                          >  8:233403/1‑149 (MQ=255)
                                                                                                                           cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCg                        >  4:33991/1‑149 (MQ=255)
                                                                                                                             cTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg                      <  5:320147/149‑1 (MQ=255)
                                                                                                                                gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGCGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCACGACGGCAAGAATGCCTAAGATTCCGGCGCCAGCCGTGGc                    >  7:296205/1‑148 (MQ=255)
                                                                                                                                  tGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATc                 >  3:195323/1‑149 (MQ=255)
                                                                                                                                    ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg               >  6:255332/5‑149 (MQ=255)
                                                                                                                                    ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg               >  1:143470/5‑149 (MQ=255)
                                                                                                                                    gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         >  2:16002/1‑107 (MQ=255)
                                                                                                                                    gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                         <  1:16002/107‑1 (MQ=255)
                                                                                                                                        aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg           <  4:223562/149‑1 (MQ=255)
                                                                                                                                        aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg           >  4:257488/1‑149 (MQ=255)
                                                                                                                                         aaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAgg          >  6:32449/1‑149 (MQ=255)
                                                                                                                                           tGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGc         >  4:153813/1‑148 (MQ=255)
                                                                                                                                            gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGcc        >  8:265609/1‑148 (MQ=255)
                                                                                                                                            gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCa       >  2:55271/1‑149 (MQ=255)
                                                                                                                                             aCCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCAt      >  4:17478/1‑149 (MQ=255)
                                                                                                                                                  gCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCCTCGCCGATGCCacggct  <  4:247072/148‑6 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCATCGCC  >  NC_002947/4431125‑4431418

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