Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,162,202 T→A 7.9% W324R (TGG→AGG) ‡ PP_4544 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,162,2020TA 7.9% 218.7 / 7.6 76W324R (TGG→AGG) ‡PP_4544hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/38);  new base A (2/4);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

CATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162054‑5162347
                                                                                                                                                    |                                                                                                                                                 
cATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCct                                                                                                                                                   >  5:136261/1‑149 (MQ=255)
 atacTGCAGGCGTTCGGCCACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGTCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctg                                                                                                                                                  <  8:382068/146‑1 (MQ=255)
   ccTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggc                                                                                                                                                >  6:322478/1‑149 (MQ=255)
    cTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggct                                                                                                                                               >  7:152948/1‑149 (MQ=255)
       cAGTCGTTCTGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACATCACCGACGGCGAAGATCAGGCctggctggc                                                                                                                                            <  7:323594/149‑1 (MQ=255)
          tCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctggctgg                                                                                                                                         <  6:96201/149‑1 (MQ=255)
             ttCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCt                                                                                                                                      >  8:66973/1‑149 (MQ=255)
              tCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCt                                                                                                                                      >  2:195500/1‑148 (MQ=255)
                ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaa                                                                                                                                   <  5:322478/149‑1 (MQ=255)
                ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                    <  1:149603/148‑1 (MQ=255)
                ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                    >  2:154271/1‑148 (MQ=255)
                ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                    >  1:349716/1‑148 (MQ=255)
                 gCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaa                                                                                                                                  <  5:2295/149‑1 (MQ=255)
                  cAACACCACCAACTAGTTCAATCTGGTGGAAGGCTTTTTCCATATAGAATTCACTGGCTTGGCCATTCCGTTTTCGGCATATGGCCAGTACGTGTAGAAAACCGACAGCAGCGACGGCTAAGATCAGGCCTGCCGGCCAGGCCTGatct                                                                                                                                 <  8:276698/149‑4 (MQ=255)
                  cAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaaa                                                                                                                                 <  2:349716/149‑1 (MQ=255)
                    acacCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTGatcttc                                                                                                                               <  4:306013/149‑6 (MQ=255)
                      accaccAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTGatcttcgc                                                                                                                             <  8:18004/149‑8 (MQ=255)
                         accaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggc                                                                                                                                                <  3:334163/127‑1 (MQ=255)
                         accaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggc                                                                                                                                                >  4:334163/1‑127 (MQ=255)
                           caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        <  8:356568/149‑1 (MQ=255)
                           caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        >  4:98676/1‑149 (MQ=255)
                             aCGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGt                                                                                                                      <  3:98676/149‑1 (MQ=255)
                             aCGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGt                                                                                                                      <  5:317365/149‑1 (MQ=255)
                                  ttCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctg                                                                                                                 >  2:342615/1‑149 (MQ=255)
                                      aTCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagc                                                                                                             <  6:136261/149‑1 (MQ=255)
                                          ggtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctgg                                                                                                         >  6:408212/1‑149 (MQ=255)
                                          ggtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          >  2:373054/1‑148 (MQ=255)
                                           gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          <  8:152948/147‑1 (MQ=255)
                                           gtggAACGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctgga                                                                                                        <  8:60177/149‑1 (MQ=255)
                                               aaGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctg                                                                                                                 <  4:264650/136‑1 (MQ=255)
                                               aaGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctg                                                                                                                 >  3:264650/1‑136 (MQ=255)
                                                        gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgc                                                                                           >  6:70582/1‑149 (MQ=255)
                                                        gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcg                                                                                            <  7:66973/148‑1 (MQ=255)
                                                        gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcg                                                                                            <  1:195500/148‑1 (MQ=255)
                                                          cAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgcgt                                                                                         >  3:301120/1‑148 (MQ=255)
                                                           aGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                        <  5:408212/149‑1 (MQ=255)
                                                              tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          <  1:398088/128‑1 (MQ=255)
                                                              tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          >  2:398088/1‑128 (MQ=255)
                                                                   ttCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCg                                                                                >  7:133143/1‑149 (MQ=255)
                                                                            ttGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAATAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGc                                                                                 <  4:204940/139‑1 (MQ=255)
                                                                            ttGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGc                                                                                 >  3:204940/1‑139 (MQ=255)
                                                                            ttGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggt                                                                       >  3:373189/1‑149 (MQ=255)
                                                                              ggCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtgg                                                                     >  1:257783/1‑149 (MQ=255)
                                                                              cgCCCTTCCGATCTCGGCCTATGGCCAGTACGTGAAGAACACCGACACCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCAGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                      <  4:282530/147‑1 (MQ=255)
                                                                                ccATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGGGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtggt                                                                    >  3:65918/1‑148 (MQ=255)
                                                                                ccATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtggt                                                                    >  3:119263/1‑148 (MQ=255)
                                                                                 cATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGGCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtggt                                                                    >  4:29872/1‑147 (MQ=255)
                                                                                   ttCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGc                                                                <  2:257783/149‑1 (MQ=255)
                                                                                           tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         <  1:279872/148‑1 (MQ=255)
                                                                                           tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  7:45427/149‑1 (MQ=255)
                                                                                           tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  1:373054/149‑1 (MQ=255)
                                                                                           tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        >  1:78170/1‑149 (MQ=255)
                                                                                           tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  7:51724/149‑1 (MQ=255)
                                                                                             ggCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCg                                                      <  1:235917/149‑1 (MQ=255)
                                                                                                cTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGAtt                                                   >  5:22748/1‑149 (MQ=255)
                                                                                                   tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCgg                                                <  1:342615/149‑1 (MQ=255)
                                                                                                   tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                 <  4:373189/148‑1 (MQ=255)
                                                                                                   tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                 <  5:70582/148‑1 (MQ=255)
                                                                                                   tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                 >  6:222791/1‑148 (MQ=255)
                                                                                                       cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACtt                                            >  1:33340/1‑149 (MQ=255)
                                                                                                       cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACtt                                            <  8:133143/149‑1 (MQ=255)
                                                                                                       cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACC‑AAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  5:100028/1‑149 (MQ=255)
                                                                                                        aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGGGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  4:65918/149‑1 (MQ=255)
                                                                                                        aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  6:57267/149‑1 (MQ=255)
                                                                                                        aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  1:213337/1‑149 (MQ=255)
                                                                                                        aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  4:119263/149‑1 (MQ=255)
                                                                                                            cGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGc                                        <  1:154271/148‑1 (MQ=255)
                                                                                                               gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGa                                               <  7:324462/138‑1 (MQ=255)
                                                                                                               gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGa                                               >  8:324462/1‑138 (MQ=255)
                                                                                                                  gaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGc                                  >  2:350798/1‑148 (MQ=255)
                                                                                                                    acacCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTac                               <  5:222791/149‑1 (MQ=255)
                                                                                                                      accgaccgCCCCCTCGGCGCAGATCAGGCCGGGCGGGCTGTCCTGTAATCCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTGAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTacac                             <  8:320815/149‑1 (MQ=255)
                                                                                                                               aCCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcg                    <  2:213337/149‑1 (MQ=255)
                                                                                                                                   aCGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGt                <  2:33340/149‑1 (MQ=255)
                                                                                                                                        gAAGATCAGGCCTGGCGGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACaa           <  5:391257/149‑1 (MQ=255)
                                                                                                                                          ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                        >  7:201152/4‑70 (MQ=25)
                                                                                                                                          ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                        <  8:201152/67‑1 (MQ=25)
                                                                                                                                          ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  6:219080/112‑1 (MQ=255)
                                                                                                                                          ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  5:219080/4‑115 (MQ=255)
                                                                                                                                            aTCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGtt        >  4:62036/1‑148 (MQ=255)
                                                                                                                                              cAGGCCTGGCTGGCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACAGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTc       >  2:401506/1‑147 (MQ=255)
                                                                                                                                                ggCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         <  1:201730/95‑1 (MQ=255)
                                                                                                                                                ggCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         >  2:201730/1‑95 (MQ=255)
                                                                                                                                                 gCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGt  >  7:155780/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162054‑5162347

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