Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,774,628 C→G 7.7% S604C (TCC→TGC) ‡ PP_1581 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,774,6280CG 7.7% 225.8 / 9.8 78S604C (TCC→TGC) ‡PP_1581BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base G (3/3);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

TCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAGTGTGTCG  >  NC_002947/1774483‑1774773
                                                                                                                                                   |                                                                                                                                                  
tgctggcgGATCTGGGTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGt                                                                                                                                                 <  1:75682/144‑1 (MQ=255)
  cTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATTATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTg                                                                                                                                               <  6:190792/149‑1 (MQ=255)
  cTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTg                                                                                                                                               <  7:89232/149‑1 (MQ=255)
  cTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTg                                                                                                                                               <  4:178870/149‑1 (MQ=255)
   tCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTgg                                                                                                                                              >  8:113211/1‑149 (MQ=255)
     gcgGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGc                                                                                                                                            <  7:113211/149‑1 (MQ=255)
      cgGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACGGAGTTGGGct                                                                                                                                           <  3:253299/149‑2 (MQ=255)
       ggATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCaa                                                                                                                                          <  2:137789/149‑1 (MQ=255)
              gCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTTCTGGAGCCCAACTCCGTTGGGCAAGTCa                                                                                                                                      <  2:142592/146‑1 (MQ=255)
              gCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCAc                                                                                                                                     <  5:235781/147‑1 (MQ=255)
              gCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACGGAGTTGGGctcca                                                                                                                                        <  3:209357/144‑5 (MQ=255)
               cTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTa                                                                                                                                  >  1:230266/1‑149 (MQ=255)
                ttGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTAc                                                                                                                                 >  8:202781/1‑149 (MQ=255)
                      gTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGAtt                                                                                                                           <  3:128572/149‑1 (MQ=255)
                      gTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGAtt                                                                                                                           >  6:390634/1‑149 (MQ=255)
                         ctcCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTGTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATtgt                                                                                                                         <  5:224988/148‑1 (MQ=255)
                          tcCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAAGCCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATtgt                                                                                                                         >  1:96379/1‑148 (MQ=255)
                          tcCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGa                                                                                                                       >  7:304135/1‑149 (MQ=255)
                            cGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGAcc                                                                                                                     >  8:270642/1‑149 (MQ=255)
                            cGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGAcc                                                                                                                     <  2:230266/149‑1 (MQ=255)
                            cGAGCGGTCGACTGTTGACCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGAcc                                                                                                                     >  8:168036/1‑149 (MQ=255)
                               gtggTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTc                                                                                                                  >  2:72565/3‑149 (MQ=255)
                                tggTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCg                                                                                                                 <  1:72565/148‑1 (MQ=255)
                                cGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCg                                                                                                                 >  7:324992/1‑149 (MQ=255)
                                 ggTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                <  8:304135/149‑1 (MQ=255)
                                 ggTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                <  4:111172/149‑1 (MQ=255)
                                   tCGACTGTTGACCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgggg                                                                                                              <  7:168036/149‑1 (MQ=255)
                                        tGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGAtt                                                                                                                           <  6:100059/131‑1 (MQ=255)
                                        tGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGAtt                                                                                                                           >  5:100059/1‑131 (MQ=255)
                                          ttGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATgg                                                                                                       >  3:208173/1‑149 (MQ=255)
                                                tgctgGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGAt                                                                                                 >  4:273883/1‑149 (MQ=255)
                                                  ctgGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCCCACAAAGTGGTGCTGGAGCCCAACGGAGTTGGGctcca                                                                                                                                        >  7:317490/1‑104 (MQ=255)
                                                  ctgGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCCCACAAAGTGGTGCTGGAGCCCAACGGAGTTGGGctcca                                                                                                                                        <  8:317490/108‑5 (MQ=255)
                                                     gAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                >  4:124144/1‑129 (MQ=255)
                                                     gAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                <  3:124144/129‑1 (MQ=255)
                                                     gAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCt                                                                                              <  8:324992/147‑1 (MQ=255)
                                                     gAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTc                                                                                             >  8:333453/1‑148 (MQ=255)
                                                         aacGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                >  7:248239/1‑125 (MQ=255)
                                                         aacGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCgg                                                                                                                <  8:248239/125‑1 (MQ=255)
                                                         aacGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATcgc                                                                                       >  1:48432/1‑149 (MQ=255)
                                                         aacGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATcgc                                                                                       <  7:333453/149‑1 (MQ=255)
                                                             atatCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTAt                                                                                   <  8:195991/149‑1 (MQ=255)
                                                              tatCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTAt                                                                                   >  7:120500/1‑148 (MQ=255)
                                                                  gagaGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACTCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACACGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGAc                                                                                >  1:281265/1‑148 (MQ=255)
                                                                  gagaGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACaa                                                                              >  1:363084/1‑149 (MQ=255)
                                                                   agagGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATAGCGTATGACaac                                                                             >  1:350166/1‑149 (MQ=255)
                                                                      ggAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACaacgaa                                                                          >  6:346659/1‑149 (MQ=255)
                                                                          aaCCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGa                                                                      <  7:270642/149‑1 (MQ=255)
                                                                              ccTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACcg                                                                   >  8:334425/1‑149 (MQ=255)
                                                                                   cGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAg                                                              <  2:96379/149‑1 (MQ=255)
                                                                                    gCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTAACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGc                                                             >  3:132088/1‑149 (MQ=255)
                                                                                      ttGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGCCGATGGGTTGATGCTCATCGCGTATGACAACGAACGGACCGCGAAGCgg                                                           >  6:336091/1‑149 (MQ=255)
                                                                                       ttgtctATGGTATGATGATGGTGGCACCTGGCCTCACAAAGTTGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGc                                                          <  7:279548/146‑1 (MQ=255)
                                                                                        gCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGc                                                          <  2:363084/148‑1 (MQ=255)
                                                                                         ctctATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCgg                                                           >  8:220069/1‑146 (MQ=255)
                                                                                                ggATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGc                                                  <  6:134628/148‑1 (MQ=255)
                                                                                                 gatgatGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGACCGTACCGCGAAGCGGCAGATTCGCAc                                                >  8:357319/1‑149 (MQ=255)
                                                                                                 gatgatGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCAc                                                >  3:208423/1‑149 (MQ=255)
                                                                                                      tgaGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGc                                             >  8:353880/1‑147 (MQ=255)
                                                                                                              tCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCAtt                                   <  8:328861/149‑1 (MQ=255)
                                                                                                              tCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCAtt                                   <  3:273883/149‑1 (MQ=255)
                                                                                                              tCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCAtt                                   <  4:105253/149‑1 (MQ=255)
                                                                                                                   cTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAg                                >  4:99517/1‑148 (MQ=255)
                                                                                                                    tGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAg                                >  8:226421/1‑147 (MQ=255)
                                                                                                                     ggCCTCACATTGTGGTGCTGCAGCCCAACTCCGTTTTGCAAGTCACGTACCCTATTGTGACCTTCGGGTACGATGGGTTGATCCTCATCGCGTTTGTCAACTAACGGTCCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCag                             <  2:350166/149‑1 (MQ=255)
                                                                                                                           acaAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTc                       <  7:202781/149‑1 (MQ=255)
                                                                                                                             aaaGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGt                     >  3:5817/1‑149 (MQ=255)
                                                                                                                                 tggtgCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCTGTGTAAGATGCGCATAGAACCACGAAGCGGCAGATTCCCAATCCCAGGGAAATTGAACAGGAGGTAGTgcc                  >  1:126978/1‑147 (MQ=255)
                                                                                                                                 tggtgCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGcc                 >  3:248745/1‑149 (MQ=255)
                                                                                                                                    tgCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGc              <  7:353880/149‑1 (MQ=255)
                                                                                                                                      cTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCAc            <  7:226421/149‑1 (MQ=255)
                                                                                                                                         ttgCCCAACGGAGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGGGGTCGTGGCCGGCACCAg         >  2:275305/3‑149 (MQ=255)
                                                                                                                                         ttgCCCAACGGAGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAg         >  6:161550/3‑149 (MQ=255)
                                                                                                                                          aGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGGCATTGAGCAGGAGGTCGTGGCCGGCACCAgt        <  8:120500/149‑1 (MQ=255)
                                                                                                                                           gCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAgtg       >  2:309579/1‑149 (MQ=255)
                                                                                                                                           gCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGAGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGAATCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAgtg       >  1:79612/1‑149 (MQ=255)
                                                                                                                                             ccAACTCCGTTGGGCAAGTCACGTACCCGATTGTGAACTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAgtgtg     >  5:102977/1‑149 (MQ=255)
                                                                                                                                                aCTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAGTGTGTCg  <  2:341363/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCCTCGCGGATCTGGCTTGGCCGTTCTCCGAGCGGTCGACTGTTGATCTGCTGGAACAACGATATCGAGAGGAAAACCCTCACGCTTGCTCTATCGGATGATGAGGGGGTCACCTGGCCTCACAAAGTGGTGCTGGAGCCCAACTCCGTTGGGCAAGTCACGTACCCGATTGTGACCTTCGGGGACGATGGGTTGATCCTCATCGCGTATGACAACGAACGGACCGCGAAGCGGCAGATTCGCACTGCCAGGGTCATTGAGCAGGAGGTCGTGGCCGGCACCAGTGTGTCG  >  NC_002947/1774483‑1774773

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: