Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,162,208 G→C 9.1% A326P (GCT→CCT) ‡ PP_4544 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,162,2080GC 9.1% 224.6 / 10.1 76A326P (GCT→CCT) ‡PP_4544hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/38);  new base C (2/5);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162061‑5162347
                                                                                                                                                   |                                                                                                                                           
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                                                                                                                                    >  1:349716/1‑148 (MQ=255)
         ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                    >  2:154271/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)
                    caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        >  4:98676/1‑149 (MQ=255)
                    caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        <  8:356568/149‑1 (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                                                                                                          >  2:398088/1‑128 (MQ=255)
                                                       tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          <  1:398088/128‑1 (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                                                        >  1:78170/1‑149 (MQ=255)
                                                                                    tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  7:51724/149‑1 (MQ=255)
                                                                                    tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  1:373054/149‑1 (MQ=255)
                                                                                    tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  7:45427/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                                                 >  6:222791/1‑148 (MQ=255)
                                                                                            tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                 <  4:373189/148‑1 (MQ=255)
                                                                                            tGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCg                                                 <  5:70582/148‑1 (MQ=255)
                                                                                                cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACtt                                            <  8:133143/149‑1 (MQ=255)
                                                                                                cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACtt                                            >  1:33340/1‑149 (MQ=255)
                                                                                                cAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACC‑AAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  5:100028/1‑149 (MQ=255)
                                                                                                 aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGGGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  4:65918/149‑1 (MQ=255)
                                                                                                 aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  1:213337/1‑149 (MQ=255)
                                                                                                 aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  6:57267/149‑1 (MQ=255)
                                                                                                 aGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  4:119263/149‑1 (MQ=255)
                                                                                                     cGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGc                                        <  1:154271/148‑1 (MQ=255)
                                                                                                        gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGa                                               >  8:324462/1‑138 (MQ=255)
                                                                                                        gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGa                                               <  7:324462/138‑1 (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                                                                                        <  8:201152/67‑1 (MQ=25)
                                                                                                                                   ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                        >  7:201152/4‑70 (MQ=25)
                                                                                                                                   ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           >  5:219080/4‑115 (MQ=255)
                                                                                                                                   ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTc                                           <  6:219080/112‑1 (MQ=255)
                                                                                                                                     aTCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGtt        >  4:62036/1‑148 (MQ=255)
                                                                                                                                       cAGGCCTGGCTGGCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACAGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTc       >  2:401506/1‑147 (MQ=255)
                                                                                                                                         ggCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         >  2:201730/1‑95 (MQ=255)
                                                                                                                                         ggCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         <  1:201730/95‑1 (MQ=255)
                                                                                                                                          gCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGt  >  7:155780/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
CAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162061‑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: ‑