Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,004,815 C→T 7.8% R466C (CGC→TGC) ‡ PP_2627 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,004,8150CT 7.8% 181.9 / 4.3 64R466C (CGC→TGC) ‡PP_2627membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base T (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CTCAGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGT  >  NC_002947/3004669‑3004957
                                                                                                                                                  |                                                                                                                                               
cTCAGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAgcgc                                                                                                                                               >  6:133842/1‑149 (MQ=255)
cTCAGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAgcgc                                                                                                                                               <  6:301051/149‑1 (MQ=255)
cTCAGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAgcgc                                                                                                                                               <  2:181563/149‑1 (MQ=255)
   aGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCAc                                                                                                                                            <  6:49980/149‑1 (MQ=255)
          cacGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGTGGCGCTTGCtgc                                                                                                                                     <  6:241669/149‑2 (MQ=255)
           acgacgAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCtgtg                                                                                                                                    <  2:65110/149‑1 (MQ=255)
           acgacgAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCtgtg                                                                                                                                    >  4:161115/1‑149 (MQ=255)
             gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtggg                                                                                                                                  >  7:5746/1‑149 (MQ=255)
             gacgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtggg                                                                                                                                  >  6:205985/1‑149 (MQ=255)
              acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 >  5:115414/1‑149 (MQ=255)
              acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 <  4:109358/149‑1 (MQ=255)
              acgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggt                                                                                                                                 <  5:133842/149‑1 (MQ=255)
               cgaAGCTGCCCGTGCGTGTAGGTAAGGCACAGGGGAAGGCGTTTGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGTGCTGGCGGGTGACAATGTTCTGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                <  2:35439/149‑1 (MQ=255)
               cgaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtg                                                                                                                                >  4:203798/1‑149 (MQ=255)
                gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                               >  1:69553/1‑149 (MQ=255)
                gaAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTgtgggtgg                                                                                                                               <  3:161115/149‑1 (MQ=255)
                       ccccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTg                                                                                                                        <  3:203798/148‑1 (MQ=255)
                        cccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgt                                                                                                                      <  2:143049/149‑1 (MQ=255)
                        cccTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGgt                                                                                                                      <  5:205985/149‑1 (MQ=255)
                          cTGCTTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGTAAGCGCCACTTGCTGTGGGTGGGCTCTGgtgt                                                                                                                    <  5:232314/149‑1 (MQ=255)
                                 aaGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATATTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATc                                                                                                             >  8:43716/1‑149 (MQ=255)
                                       aGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGc                                                                                                        <  2:69553/148‑1 (MQ=255)
                                          cacaGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGtt                                                                                                     >  7:60984/1‑148 (MQ=255)
                                               gggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTAt                                                                                               >  5:321245/1‑149 (MQ=255)
                                                ggcaGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATc                                                                                              <  4:250337/146‑1 (MQ=255)
                                                ggAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATc                                                                                              <  6:115414/149‑1 (MQ=255)
                                                   aGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGc                                                                                            <  8:60984/148‑1 (MQ=255)
                                                      cGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGCAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGtt                                                                                        >  1:234460/1‑149 (MQ=255)
                                                           gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGc                                                                                      <  7:56419/146‑1 (MQ=255)
                                                           gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGc                                                                                      >  8:56419/1‑146 (MQ=255)
                                                           gTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAg                                                                                   >  2:94366/1‑149 (MQ=255)
                                                              tATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTg                                                                                >  1:78492/1‑149 (MQ=255)
                                                                tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  4:234502/1‑148 (MQ=255)
                                                                tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCt                                                                                 <  6:321245/146‑1 (MQ=255)
                                                                tttCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  4:137979/1‑148 (MQ=255)
                                                                    aTCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCg                                                                          >  4:130061/1‑149 (MQ=255)
                                                                     tCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGt                                                                         >  5:336384/1‑149 (MQ=255)
                                                                           aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               >  1:322921/1‑137 (MQ=255)
                                                                           aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTgg                                                                               <  2:322921/137‑1 (MQ=255)
                                                                           aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGcc                                                                   <  2:78492/149‑1 (MQ=255)
                                                                           aagcaTTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGcc                                                                   >  1:328750/1‑149 (MQ=255)
                                                                                ttCAGGCGGGTTATCTATCAGGAAGCAGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACa                                                               >  1:178395/1‑148 (MQ=255)
                                                                                   aGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCaa                                                           >  7:96686/1‑149 (MQ=255)
                                                                                    ggCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                          <  3:130061/149‑1 (MQ=255)
                                                                                    ggCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAAt                                                          >  7:97287/1‑149 (MQ=255)
                                                                                      cGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATgg                                                        <  6:336384/149‑1 (MQ=255)
                                                                                       gggTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgc                                                       <  1:94366/149‑1 (MQ=255)
                                                                                        ggTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGgcg                                                      <  8:5746/149‑1 (MQ=255)
                                                                                                atcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCAAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCt                                              >  8:295842/1‑149 (MQ=255)
                                                                                                 tcAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTGGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAg                                                 <  2:234460/145‑1 (MQ=255)
                                                                                                  cAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            >  6:156056/1‑149 (MQ=255)
                                                                                                   aGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            <  8:97287/148‑1 (MQ=255)
                                                                                                   aGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGc                                            >  2:269968/1‑148 (MQ=255)
                                                                                                     gAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAg                                         >  8:222561/1‑149 (MQ=255)
                                                                                                                   gggTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCaaaa                           <  2:262479/149‑1 (MQ=255)
                                                                                                                   gggTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCaaaa                           <  1:269968/149‑1 (MQ=255)
                                                                                                                        aCATTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                      <  3:234502/149‑1 (MQ=255)
                                                                                                                          aaTGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAg                      >  5:64650/1‑147 (MQ=255)
                                                                                                                             gTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTa                 >  4:71456/1‑149 (MQ=255)
                                                                                                                                    ttGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGt          >  8:88470/1‑149 (MQ=255)
                                                                                                                                        gcgcAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc       >  8:90622/1‑148 (MQ=255)
                                                                                                                                        ccacagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCAc       >  1:227172/4‑149 (MQ=255)
                                                                                                                                          gcagcaAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGaa    <  2:328750/149‑1 (MQ=255)
                                                                                                                                          acagcaAGTGGCGCTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACg      >  6:222252/2‑147 (MQ=255)
                                                                                                                                            agcaAGTGGCGCTTGCTGCGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  2:53207/149‑1 (MQ=255)
                                                                                                                                            agcaAGTGGCGCTTGCTGCGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  2:48347/149‑1 (MQ=255)
                                                                                                                                            agcaAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGt  <  6:64650/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTCAGCCCTACACGACGAAGCTGCCCCTGCGTGAAGGTAAGGCACAGGGGAAGGCGTTGGTCTATTTCATCCAGCAAGCATTCAGGCGGGTTATCTATCAGGAAGCCGGGCTGGCGGGTGACAATGTTCGGGTTGCGCGCAGCAAGCGCCACTTGCTGTGGGTGGGCTCTGGTGTAGGCATCCTGGCGTTCAGTATCGCCGTTGCCAGCTGGCAGCGTTACGCCGACATCAATGGCGCCAAGGCTGCCAGTGTGCTCGCCAAAAGCCAGGAGTACAGCGGTCACGAAGT  >  NC_002947/3004669‑3004957

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: