Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 422,065 A→C 6.0% I628S (ATC→AGC)  PP_0348 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947422,0650AC 6.0% 249.9 / 7.0 83I628S (ATC→AGC) PP_0348hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base C (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

TCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGT  >  NC_002947/421918‑422209
                                                                                                                                                   |                                                                                                                                                 
tCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGa                                                                                                                                                   <  4:276218/148‑1 (MQ=255)
  gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑t                                                                                                                                                 >  3:262004/1‑147 (MQ=255)
  gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑Tc                                                                                                                                                <  4:195187/148‑1 (MQ=255)
  gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCa                                                                                                                                               <  7:73990/149‑1 (MQ=255)
    ttGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAgc                                                                                                                                             <  3:389078/149‑1 (MQ=255)
        cATGCGGCCAATCGACAGCACCAGCAAAGCGCGGATGGGCAGAAGCAGCCTGATGAAATCCAGCAAACGCCAGCTCAGAGCCATGACCGACTCCGCAGTCACTTCCATGGACGTAATGCTGTCGGCGTCGGACTTGGCCA‑TCAGAGccc                                                                                                                                         >  6:226080/1‑149 (MQ=255)
          tGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCgg                                                                                                                                       <  5:99910/149‑1 (MQ=255)
            cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGc                                                                                                                                     >  8:371890/1‑149 (MQ=255)
            cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGC‑TGATCGCCCagg                                                                                                                                      <  5:363951/148‑1 (MQ=255)
             ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGcc                                                                                                                                    >  6:338896/1‑149 (MQ=255)
               ccAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTc                                                                                                                                  >  1:261171/1‑149 (MQ=255)
                cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCg                                                                                                                                 >  2:330802/1‑149 (MQ=255)
                 actCGACAGCACCCGCAACGCGCGGCTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                <  8:300397/147‑1 (MQ=255)
                 aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                >  2:76212/1‑149 (MQ=255)
                 aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGGCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCgg                                                                                                                                <  5:48294/149‑1 (MQ=255)
                        aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                         <  1:213276/149‑1 (MQ=255)
                        aGCACCAGCAACGCGCGGGTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGGTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgta                                                                                                                         <  6:209178/149‑1 (MQ=255)
                          cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                        >  6:170506/1‑148 (MQ=255)
                          cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                        >  5:77799/1‑148 (MQ=255)
                               gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCg                                                                                                                  <  2:261171/149‑1 (MQ=255)
                                     cgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaa                                                                                                            <  7:371890/149‑1 (MQ=255)
                                      gcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                             >  6:357459/1‑147 (MQ=255)
                                        ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                         >  2:279675/1‑149 (MQ=255)
                                        ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                         <  5:338896/149‑1 (MQ=255)
                                        ggTTGGGCAGCAGCAGGATGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCTCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCTA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                         <  5:259167/149‑1 (MQ=255)
                                              gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                   >  8:256229/1‑149 (MQ=255)
                                                agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                                 <  6:39341/149‑1 (MQ=255)
                                                         cTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCggg                                                                                        <  5:170506/149‑1 (MQ=255)
                                                           gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                       >  3:359345/1‑132 (MQ=255)
                                                           gatgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                       <  4:359345/132‑1 (MQ=255)
                                                            atgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGt                                                                                     <  5:357459/149‑1 (MQ=255)
                                                             tgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                    >  7:29160/1‑149 (MQ=255)
                                                               aaaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTgg                                                                                   <  7:528/148‑1 (MQ=255)
                                                                aaTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGt                                                                                  >  8:365945/1‑148 (MQ=255)
                                                                 aTCCAGCAAACGCCTGCTCAGCGCCATGACCGACTCCGCATTCACTTCCGTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCTAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                                 <  5:226080/148‑1 (MQ=255)
                                                                 aTCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTc                                                                                 >  5:378088/1‑148 (MQ=255)
                                                                   ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                           <  7:30933/120‑1 (MQ=255)
                                                                   ccAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGaaaa                                                                                                           >  8:30933/1‑120 (MQ=255)
                                                                    cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                             >  4:329098/1‑149 (MQ=255)
                                                                    cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAAt                                                                             >  7:313245/1‑149 (MQ=255)
                                                                    cAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                              <  1:76212/148‑1 (MQ=255)
                                                                     aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                              <  5:390894/147‑1 (MQ=255)
                                                                     aGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                              >  6:390894/1‑147 (MQ=255)
                                                                      gCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGc                                                                           <  5:219926/149‑1 (MQ=255)
                                                                           cGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGc                                                                       >  1:274771/1‑148 (MQ=255)
                                                                             ccAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCa                                                                    <  2:274771/149‑1 (MQ=255)
                                                                                  tCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                               >  7:88021/1‑149 (MQ=255)
                                                                                    aGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                               <  1:330802/147‑1 (MQ=255)
                                                                                     tcgcCATGACCGACTCCGCATTCACTTCGCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                            <  8:346828/148‑1 (MQ=255)
                                                                                        ccATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTg                                                         <  7:256229/149‑1 (MQ=255)
                                                                                         cATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACTTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGt                                                        >  3:127739/1‑149 (MQ=255)
                                                                                          aTGACCGACTCCGCATTCACTTCCCTGGTCGACATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTc                                                       >  7:298382/1‑149 (MQ=255)
                                                                                           tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                      <  8:313245/149‑1 (MQ=255)
                                                                                           tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                      <  7:365945/149‑1 (MQ=255)
                                                                                              ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                       <  6:147911/81‑1 (MQ=255)
                                                                                              ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATgtagt                                                                                                                       >  5:147911/1‑81 (MQ=255)
                                                                                               cGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGc                                                   >  6:377328/1‑148 (MQ=255)
                                                                                                   tCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCg                                              >  1:318349/1‑149 (MQ=255)
                                                                                                     cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                       >  8:204825/1‑106 (MQ=255)
                                                                                                     cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                       <  7:204825/106‑1 (MQ=255)
                                                                                                         ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCAc                                        <  3:124779/149‑1 (MQ=255)
                                                                                                            aCTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGCCCTCCGCCATCTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGCCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCTGCAACCGCTGCACCGt                                     >  5:78585/1‑149 (MQ=255)
                                                                                                             cTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTGAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                    <  1:279675/149‑1 (MQ=255)
                                                                                                             cTTCCCTGGTCGTCATGCTGTCGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTg                                    <  8:29160/149‑1 (MQ=255)
                                                                                                              ttCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                   >  3:46803/1‑149 (MQ=255)
                                                                                                                cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGGCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTGGGGCGGGGGCAAATGCGTGCTCAGGCGCCCCTTTTCCAGCAACCGCTGCACCGTGGCg                                 >  2:188409/1‑149 (MQ=255)
                                                                                                                cccTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                  <  6:164764/148‑1 (MQ=255)
                                                                                                                cccTGGGCGTCATGCTGACGGCGGCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGc                                  <  8:273413/148‑1 (MQ=255)
                                                                                                                    ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGAGTCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                                >  1:47358/1‑147 (MQ=255)
                                                                                                                     gtcgtcATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGGCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAg                             <  6:77799/148‑1 (MQ=255)
                                                                                                                            tGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAAt                     >  5:8515/1‑149 (MQ=255)
                                                                                                                              cTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATgg                   <  6:8515/149‑1 (MQ=255)
                                                                                                                               tGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGc                  <  3:235763/149‑1 (MQ=255)
                                                                                                                                  cggcgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACGTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAgg               >  6:167553/1‑149 (MQ=255)
                                                                                                                                     cgtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcgc            >  2:329589/1‑149 (MQ=255)
                                                                                                                                     ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                   >  6:299022/5‑126 (MQ=255)
                                                                                                                                     ccgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                   <  5:299022/122‑1 (MQ=255)
                                                                                                                                      gtcggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCt           <  2:318349/149‑1 (MQ=255)
                                                                                                                                      cgaggcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGAGCTCAGGCGCCCCTTGTCCCGCAACAGCTGCACCGTGGCGACAGGCAGCATTGGCAGGCGCt           >  4:238413/4‑149 (MQ=255)
                                                                                                                                         ggcCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTgg         <  4:46803/148‑1 (MQ=255)
                                                                                                                                          gcCCGGGCGC‑TGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc       >  7:219290/1‑149 (MQ=255)
                                                                                                                                            cTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta     >  6:205116/1‑149 (MQ=255)
                                                                                                                                             tGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtac    >  3:377633/1‑149 (MQ=255)
                                                                                                                                               ggCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACTTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtacgt  <  4:127739/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGA‑TCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGT  >  NC_002947/421918‑422209

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