Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 590,361 G→A 6.3% G162D (GGC→GAC)  PP_0500 → dTDP‑4‑rhamnose reductase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947590,3610GA 6.3% 189.4 / 2.0 63G162D (GGC→GAC) PP_0500dTDP‑4‑rhamnose reductase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base A (2/2);  total (27/36)
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.66e-01

GCGAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCTGTGGGGCACCTAC  >  NC_002947/590220‑590509
                                                                                                                                             |                                                                                                                                                    
gCGAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTc                                                                                                                                               <  2:10000/149‑1 (MQ=255)
 cGAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTcc                                                                                                                                              <  5:188062/149‑1 (MQ=255)
  gAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCt                                                                                                                                             <  8:150681/149‑1 (MQ=255)
  gAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCt                                                                                                                                             <  6:129455/149‑1 (MQ=255)
    aaaGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGa                                                                                                                                           <  6:287738/149‑1 (MQ=255)
      aGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGAcc                                                                                                                                         <  2:244066/149‑1 (MQ=255)
       ggACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGAccc                                                                                                                                        <  5:20009/149‑1 (MQ=255)
       ggACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGAccc                                                                                                                                        <  5:24567/149‑1 (MQ=255)
              ccAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGa                                                                                                                                 <  6:137297/149‑1 (MQ=255)
                aGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGa                                                                                                                                           >  4:259447/1‑137 (MQ=255)
                aGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGa                                                                                                                                           <  3:259447/137‑1 (MQ=255)
                    ccGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAAc                                                                                                                            <  6:220788/148‑1 (MQ=255)
                           gCTTGCGTGGCCAGGGCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAg                                                                                                                    <  8:707/149‑1 (MQ=255)
                             ttGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCt                                                                                                                                             <  4:218353/122‑1 (MQ=255)
                             ttGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCt                                                                                                                                             >  3:218353/1‑122 (MQ=255)
                                  tGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGa                                                                                                                     <  2:233351/141‑1 (MQ=255)
                                  tGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGa                                                                                                                     >  1:233351/1‑141 (MQ=255)
                                  tGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGACCGCTTCCTGACCCGTGCCGAGCAACCACAGGAgctgctgc                                                                                                             >  6:161766/1‑149 (MQ=255)
                                  tGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGACCGCTTCCTGACCCGTGCCGAGCAACCACAGGAgctgctgc                                                                                                             >  6:165934/1‑149 (MQ=255)
                                      cAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGc                                                                                                         <  8:100232/149‑1 (MQ=255)
                                           ccTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCGTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCgat                                                                                                     >  1:85310/1‑148 (MQ=255)
                                               tGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATc                                                                                                 <  3:39904/148‑1 (MQ=255)
                                               tGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGACCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATc                                                                                                 <  5:165934/148‑1 (MQ=255)
                                               tGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGACCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATc                                                                                                 <  5:161766/148‑1 (MQ=255)
                                                   gcATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGt                                                                                            >  5:242032/1‑149 (MQ=255)
                                                   gcATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGt                                                                                            <  6:225027/149‑1 (MQ=255)
                                                    cATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTg                                                                                           <  6:141811/149‑1 (MQ=255)
                                                        gAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCaa                                                                                       <  2:185209/149‑1 (MQ=255)
                                                        gAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCCATGATCGACGTGGCaa                                                                                       >  2:218185/1‑149 (MQ=255)
                                                         aGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAAc                                                                                      >  1:159559/1‑149 (MQ=255)
                                                         aGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAAc                                                                                      >  5:150729/1‑149 (MQ=255)
                                                         aGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAAc                                                                                      >  4:284895/1‑149 (MQ=255)
                                                                   gcgtgcgGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCCATGGCGCATTCGGACGCTACCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCgg                                                                            >  6:297233/1‑149 (MQ=255)
                                                                       gcgGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCga                                                                        <  2:159559/149‑1 (MQ=255)
                                                                          gCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCgatga                                                                     >  7:230159/1‑149 (MQ=255)
                                                                           ccTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCg                                                                             <  6:243471/140‑1 (MQ=255)
                                                                           ccTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCg                                                                             >  5:243471/1‑140 (MQ=255)
                                                                            cTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTAGTTGCTGGACGAGAGCATCGATGGCGCATGGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGacgatg                                                                   <  1:200935/149‑1 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTggg                                                                                                                                                      >  8:271801/1‑63 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTggg                                                                                                                                                      <  7:271801/63‑1 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGc                                                                               <  6:197014/134‑1 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGc                                                                               >  5:197014/1‑134 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCgatgat                                                                    <  2:135328/145‑1 (MQ=255)
                                                                               cccGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCgatgat                                                                    >  1:135328/1‑145 (MQ=255)
                                                                                    agcaTGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGa                                                           <  4:207015/149‑1 (MQ=255)
                                                                                     gcaTGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCCATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAg                                                          <  1:218185/149‑1 (MQ=255)
                                                                                          tgtTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATc                                                     <  6:242032/149‑1 (MQ=255)
                                                                                           gtTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAAcc                                                                                     <  6:208254/116‑1 (MQ=255)
                                                                                           gtTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGt                                                         <  4:201355/144‑1 (MQ=255)
                                                                                           gtTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGt                                                         >  3:201355/1‑144 (MQ=255)
                                                                                              gctgcGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGt                                                 >  8:219708/1‑149 (MQ=255)
                                                                                              gctgcGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGt                                                 >  5:99219/1‑149 (MQ=255)
                                                                                               ctgcGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGt                                                 <  3:284895/148‑1 (MQ=255)
                                                                                                tgcGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCg                                               >  8:138238/1‑149 (MQ=255)
                                                                                                           tGGTTGCTGGACGAGAGCATAGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAAACCGACGCCGGTCGATGATGCCGCACGAGTGATCCTGTCGGTGCTCAAGc                                     >  7:173606/1‑148 (MQ=255)
                                                                                                                  tGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCgg                                              >  2:94337/1‑132 (MQ=255)
                                                                                                                  tGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCgg                                              <  1:94337/132‑1 (MQ=255)
                                                                                                                    gACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTg                           >  6:167388/1‑149 (MQ=255)
                                                                                                                       gagaGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCaa                        <  7:138238/149‑1 (MQ=255)
                                                                                                                         gagCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATg                      >  1:510/1‑149 (MQ=255)
                                                                                                                             aTCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATgc                     >  4:201019/1‑146 (MQ=255)
                                                                                                                             aTCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATgc                     <  3:201019/146‑1 (MQ=255)
                                                                                                                              tCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGcgccg                 >  5:14607/1‑149 (MQ=255)
                                                                                                                                gATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCt               >  4:170647/1‑149 (MQ=255)
                                                                                                                                    gcgcATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCtgtg            <  8:230159/148‑1 (MQ=255)
                                                                                                                                           gggCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCTGTGGGGCAcc     >  8:104021/1‑148 (MQ=255)
                                                                                                                                             gCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCTGTGGGGCACCTAc  <  5:167388/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCGAAAAGGACGAGCCAGTGCCGCTGGGCTTGCGTGGCCAGGCCCTATGGCGCATCGAGCAGAGCGTGCGTGCGGCCTGCCCGCAGCATGTGTTGCTGCGCTTTGGCTGGTTGCTGGACGAGAGCATCGATGGCGCATTGGGCCGCTTCCTGACCCGTGCCGAGCAACCACAGGAGCTGCTGCTGGCCGATGATCGACGTGGCAACCCGACGCCGGTCGATGATGCCGCCCGAGTGATCCTGTCGGTGCTCAAGCAACTGGATTGCAATGCGCCGCTGTGGGGCACCTAC  >  NC_002947/590220‑590509

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