Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,577,644 T→A 6.0% L20Q (CTG→CAG)  PP_3158 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,577,6440TA 6.0% 210.5 / 3.0 67L20Q (CTG→CAG) PP_3158hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base A (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

CTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTG  >  NC_002947/3577498‑3577789
                                                                                                                                                     |                                                                                                                                                 
ctcCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTgg                                                                                                                                                 >  7:326427/1‑149 (MQ=255)
 tcCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTgg                                                                                                                                                 >  5:61465/1‑148 (MQ=255)
   cAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                              <  6:169346/149‑1 (MQ=255)
    aGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                              >  4:278657/1‑148 (MQ=255)
    aGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                              <  3:278657/148‑1 (MQ=255)
     gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                            <  8:326427/149‑1 (MQ=255)
     gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                            <  8:109335/149‑1 (MQ=255)
     gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                            <  2:166208/149‑1 (MQ=255)
      ccACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                              <  1:388858/146‑1 (MQ=255)
      ccACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                              >  2:388858/1‑146 (MQ=255)
       cACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCt                                                                                                                                          <  1:186146/149‑1 (MQ=255)
          cGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCtgtg                                                                                                                                       >  1:109937/1‑149 (MQ=255)
             aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                    >  7:179686/1‑149 (MQ=255)
             aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                    >  6:20488/1‑149 (MQ=255)
                 tttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTg                                                                                                                                 >  7:105352/1‑148 (MQ=255)
                  ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGt                                                                                                                                <  8:179686/148‑1 (MQ=255)
                  ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTa                                                                                                                               >  6:160556/1‑149 (MQ=255)
                   tCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTAc                                                                                                                              >  3:244080/1‑149 (MQ=255)
                     ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCTACCCTGTTCCGATGGTTGCAAGGTACTCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                             <  8:295010/148‑1 (MQ=255)
                     ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                             >  8:270317/1‑148 (MQ=255)
                      cATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGCGCACAGGCATGCCAGGGCACCTGtcagta                                                                                                                           <  6:219781/149‑6 (MQ=38)
                       aTTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACgcag                                                                                                                          <  1:266494/149‑1 (MQ=255)
                             gggCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGcagcagcag                                                                                                                    <  2:109937/149‑1 (MQ=255)
                                ctctATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGc                                                                                                                 >  4:114196/1‑149 (MQ=255)
                                   tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGcc                                                                                                                <  6:256792/147‑1 (MQ=255)
                                   tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAAGTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                              >  1:87417/1‑149 (MQ=255)
                                     tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCg                                                                                                               <  2:271396/146‑1 (MQ=255)
                                     tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                              <  7:39925/147‑1 (MQ=255)
                                      gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                   >  3:137270/1‑140 (MQ=18)
                                      gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                   <  4:137270/144‑5 (MQ=18)
                                           gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATgg                                                                                                       <  4:244080/148‑1 (MQ=255)
                                           gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                      <  6:61465/149‑1 (MQ=255)
                                                gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                      <  4:185055/144‑1 (MQ=255)
                                                gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                      >  3:185055/1‑144 (MQ=255)
                                                      cTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcg                                                                                           <  7:270317/149‑1 (MQ=255)
                                                           tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCGTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                      >  7:218030/1‑149 (MQ=255)
                                                           tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                      >  4:183299/1‑149 (MQ=255)
                                                             ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                    >  6:242441/1‑149 (MQ=255)
                                                             ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacg                                                                                     <  5:160556/148‑1 (MQ=255)
                                                                gcgcAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacggca                                                                                  >  1:379792/1‑148 (MQ=255)
                                                                    aGCCGCCCAATATCGACCCGTAGCCGAAGGTAGCAAGCTACGCGAATTTCCTACCCCGCCTGCTACTG‑‑‑ACAGTTGCCCTGGCATGCCTGTGCGCCTTTACGCATCAGCAGGTCCGTGACATGGTCAACCAGCTCGGCAACGTCAAAgcc                                                                             <  7:166119/149‑1 (MQ=255)
                                                                            aaTATCGACCCGGTGCCGATGGTTGCAAGTTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCt                                                                     >  6:215466/1‑149 (MQ=255)
                                                                               atCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGtt                                                                  >  6:374055/1‑149 (MQ=255)
                                                                                  gACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCa                                                               <  8:105352/149‑1 (MQ=255)
                                                                                    cccGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGa                                                             <  2:379792/149‑1 (MQ=255)
                                                                                          gCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCg                                                       <  7:72046/149‑1 (MQ=255)
                                                                                                   tGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                    <  2:281141/111‑1 (MQ=255)
                                                                                                     cAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTg                                            <  3:114196/149‑1 (MQ=255)
                                                                                                      aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                           <  5:215466/149‑1 (MQ=255)
                                                                                                      aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                           >  2:102131/1‑149 (MQ=255)
                                                                                                            cgcgAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGc                                     >  1:1461/1‑149 (MQ=255)
                                                                                                                aaTTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTc                                 <  7:120340/149‑1 (MQ=255)
                                                                                                                       tccCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                          >  6:15830/1‑149 (MQ=255)
                                                                                                                         ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCTCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGTCCGCCTGGCCACGCTGTCGATgcgc                          <  8:138593/147‑1 (MQ=255)
                                                                                                                         ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                          <  8:218030/147‑1 (MQ=255)
                                                                                                                         ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGcgcg                         >  4:365268/1‑148 (MQ=255)
                                                                                                                               cTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCt                  >  1:93426/1‑149 (MQ=255)
                                                                                                                                  ccaCTG‑‑‑ACCCGGGCCCTTGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGGGACATGGGCAACCACCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa               <  5:144422/147‑1 (MQ=255)
                                                                                                                                   tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa               >  2:128631/1‑148 (MQ=255)
                                                                                                                                   tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGACTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACAGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa               >  4:93652/1‑148 (MQ=255)
                                                                                                                                        cgtaCAGGCGCACAGGCATGCCAGGGCACCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAg            >  7:43058/4‑149 (MQ=255)
                                                                                                                                           aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                          <  4:69522/84‑1 (MQ=255)
                                                                                                                                           aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                          >  3:69522/1‑84 (MQ=255)
                                                                                                                                           aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGc           <  2:1461/147‑1 (MQ=255)
                                                                                                                                            cAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGGCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTg        <  8:376579/149‑1 (MQ=255)
                                                                                                                                                  cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTAGctc  >  7:140944/1‑148 (MQ=255)
                                                                                                                                                  cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACctg  <  5:15830/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTG  >  NC_002947/3577498‑3577789

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