Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,162,204 G→C 9.2% W324C (TGG→TGC) ‡ PP_4544 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,162,2040GC 9.2% 227.9 / 13.0 76W324C (TGG→TGC) ‡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 = 1.43e-01

CCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162057‑5162347
                                                                                                                                                   |                                                                                                                                               
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                                                                                                                                    >  2:154271/1‑148 (MQ=255)
             ggCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGa                                                                                                                                    <  1:149603/148‑1 (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                                                                                                                                                >  4:334163/1‑127 (MQ=255)
                      accaACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggc                                                                                                                                                <  3:334163/127‑1 (MQ=255)
                        caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        >  4:98676/1‑149 (MQ=255)
                        caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        <  8:356568/149‑1 (MQ=255)
                          aCGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGt                                                                                                                      <  5:317365/149‑1 (MQ=255)
                          aCGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGt                                                                                                                      <  3:98676/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                                                                                                                 >  3:264650/1‑136 (MQ=255)
                                            aaGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctg                                                                                                                 <  4:264650/136‑1 (MQ=255)
                                                     gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgc                                                                                           >  6:70582/1‑149 (MQ=255)
                                                     gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcg                                                                                            <  1:195500/148‑1 (MQ=255)
                                                     gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcg                                                                                            <  7:66973/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:373054/149‑1 (MQ=255)
                                                                                        tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        >  1:78170/1‑149 (MQ=255)
                                                                                        tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCAc                                                        <  7:51724/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                                            >  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                                           >  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                                                                                        >  7:201152/4‑70 (MQ=25)
                                                                                                                                       ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAt                                                                                        <  8:201152/67‑1 (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)
                                                                                                                                                   |                                                                                                                                               
CCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162057‑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: ‑