Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,646,266 A→T 5.1% D540V (GAC→GTC)  PP_4954 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,646,2660AT 5.1% 239.1 / 2.7 78D540V (GAC→GTC) PP_4954membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/33);  new base T (2/2);  total (43/35)
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.74e-01

CTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGC  >  NC_002947/5646124‑5646410
                                                                                                                                              |                                                                                                                                                
cTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTgg                                                                                                                                              >  6:288877/1‑147 (MQ=255)
   gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                         <  8:32207/149‑1 (MQ=255)
   gCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCGGTTCCTGCCTGACTGGCAggg                                                                                                                                         >  6:109805/1‑149 (MQ=255)
    ccACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcaggaa                                                                                                                                        <  5:263118/149‑3 (MQ=255)
      aCTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCAGTcaggcaggaaca                                                                                                                                      <  2:144942/149‑5 (MQ=255)
          ccGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTCTTCCTGCCTGACTGGCAGGGCCGACgg                                                                                                                                  <  1:15763/149‑1 (MQ=255)
               aTCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                             >  2:240399/1‑149 (MQ=255)
               aTCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                             >  3:6011/1‑149 (MQ=255)
                   cccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAg                                                                                                                         >  8:208831/1‑149 (MQ=255)
                   cccTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAg                                                                                                                         <  3:120489/149‑1 (MQ=255)
                        aTGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCt                                                                                                                    <  2:239256/149‑1 (MQ=255)
                           gTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGc                                                                                                                 <  2:250613/149‑1 (MQ=255)
                           gTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGc                                                                                                                 <  1:227662/149‑1 (MQ=255)
                                 ttCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATAc                                                                                                           >  2:29042/1‑149 (MQ=255)
                                    tGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCAGGCCagggtatt                                                                                                        >  1:55335/1‑142 (MQ=255)
                                      cTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATAccc                                                                                                         <  1:240399/146‑1 (MQ=255)
                                       ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                         >  4:115259/1‑113 (MQ=255)
                                       ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAggg                                                                                                                                         <  3:115259/113‑1 (MQ=255)
                                       ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATCCCCTGGc                                                                                                     >  6:187208/1‑149 (MQ=255)
                                       ttCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCGGAACAAGGCGCTGGCCAATACCCTgg                                                                                                      >  5:315589/1‑148 (MQ=255)
                                          aaCCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTg                                                                                                  >  2:222556/1‑149 (MQ=255)
                                           aCCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCTCCTGTTCGATACCCTGGGCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgc                                                                                                 <  7:223966/149‑1 (MQ=255)
                                            ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                             >  3:232719/1‑120 (MQ=255)
                                            ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGaa                                                                                                                             <  4:232719/120‑1 (MQ=255)
                                            ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                       <  6:90728/142‑1 (MQ=255)
                                            ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                       >  5:90728/1‑142 (MQ=255)
                                            ccAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgcg                                                                                                >  6:206218/1‑149 (MQ=255)
                                             cAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgcgc                                                                                               >  5:277097/1‑149 (MQ=255)
                                               gATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTgcgc                                                                                               >  6:222729/1‑147 (MQ=255)
                                                 tCGGCGCTGGCTACGGCCTGTTCCGGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGc                                                                                           <  6:19001/149‑1 (MQ=255)
                                                 tCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGc                                                                                           <  6:277097/149‑1 (MQ=255)
                                                    gCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACgg                                                                                                                                  >  5:278681/1‑107 (MQ=255)
                                                    gCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACgg                                                                                                                                  <  6:278681/107‑1 (MQ=255)
                                                     cGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                       >  6:268022/1‑133 (MQ=255)
                                                     cGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTg                                                                                                       <  5:268022/133‑1 (MQ=255)
                                                       aTGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTAc                                                                                     >  8:111086/1‑149 (MQ=255)
                                                          gCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTg                                                                                  >  6:303083/1‑149 (MQ=255)
                                                           cTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                 >  1:54120/1‑149 (MQ=255)
                                                           cTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgc                                                                                 >  7:5227/1‑149 (MQ=255)
                                                              cGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgcgc                                                                               <  2:54120/148‑1 (MQ=255)
                                                                  cTGTTCCTGCCGCGCCTGTTGGATACCCTGGTCGGCAGCCGTATTGCAATTCTTGCAGTGTTCCTGTTCCTGCCGTACTGGATGGGCCTTCGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGAt                                                                          <  1:95605/149‑1 (MQ=255)
                                                                       ccTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      <  1:222556/148‑1 (MQ=255)
                                                                       ccTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      <  5:222729/148‑1 (MQ=255)
                                                                        cTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      >  5:141088/1‑147 (MQ=255)
                                                                        cTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGACTGACTGGCAGGGCCGACGGCTGAACAAGGAGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTTCCTGCGCCAGATCATg                                                                      >  5:38039/1‑147 (MQ=255)
                                                                         tGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTCATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATgcag                                                                   >  3:321327/1‑149 (MQ=255)
                                                                           ccGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATgcagca                                                                 <  5:206218/149‑1 (MQ=255)
                                                                             gcgcCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTTTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGt                                                               <  2:271669/149‑1 (MQ=255)
                                                                               gcCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTAc                                                                                     <  5:120516/125‑1 (MQ=255)
                                                                               gcCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTAc                                                                                     >  6:120516/1‑125 (MQ=255)
                                                                                       gATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgcgc                                                                               <  4:291994/123‑1 (MQ=255)
                                                                                       gATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTgcgc                                                                               >  3:291994/1‑123 (MQ=255)
                                                                                       gATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACgg                                                     >  8:263215/1‑149 (MQ=255)
                                                                                           cccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCaa                                                  <  2:99313/148‑1 (MQ=255)
                                                                                            ccTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAgc                                                >  4:306615/1‑149 (MQ=255)
                                                                                                gTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGcga                                            <  4:6011/149‑1 (MQ=255)
                                                                                                    gCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCgacgac                                        >  6:188578/1‑149 (MQ=255)
                                                                                                    gCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCgacgac                                        <  6:291495/149‑1 (MQ=255)
                                                                                                        ccTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGGCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTgg                                    >  4:319463/1‑149 (MQ=255)
                                                                                                               gCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCg                             >  4:159488/1‑149 (MQ=255)
                                                                                                                ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      >  3:29320/1‑107 (MQ=255)
                                                                                                                ccATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATg                                                                      <  4:29320/107‑1 (MQ=255)
                                                                                                                         cAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCgg                   >  4:86616/1‑149 (MQ=255)
                                                                                                                           gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg                 <  7:111086/149‑1 (MQ=255)
                                                                                                                           gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg                 >  5:75250/1‑149 (MQ=255)
                                                                                                                           gtgtTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGgcg                 >  1:130685/1‑149 (MQ=255)
                                                                                                                              ttcctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCaa              >  1:65932/1‑149 (MQ=255)
                                                                                                                                cctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCgg                   <  5:215433/142‑1 (MQ=255)
                                                                                                                                cctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCgg                   >  6:215433/1‑142 (MQ=255)
                                                                                                                                 ctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGc           >  8:204257/1‑149 (MQ=255)
                                                                                                                                 ctgttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGc           >  1:6239/1‑149 (MQ=255)
                                                                                                                                    ttcctgCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca        <  5:303083/149‑1 (MQ=255)
                                                                                                                                     ccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca        >  3:3412/2‑148 (MQ=255)
                                                                                                                                     ccctgccagTCAGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCca        >  7:196573/2‑148 (MQ=255)
                                                                                                                                          cctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAAc    <  5:288877/147‑1 (MQ=255)
                                                                                                                                          cctgACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGc  <  8:5227/149‑1 (MQ=255)
                                                                                                                                             gACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTgg                        >  3:97293/1‑124 (MQ=255)
                                                                                                                                             gACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTgg                        <  4:97293/124‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
CTGGCCACTGCCGCGATCACCCTGATGGTGCTGTTCTGCTTCAACCAGATCGGCGATGGCTACGGCCTGTTCCTGCCGCGCCTGTTCGATACCCTGGTCGGCAGCCTGATTGCCATTCTTGCAGTGTTCCTGTTCCTGCCTGACTGGCAGGGCCGACGGCTGAACAAGGCGCTGGCCAATACCCTGGCCTGCGCCAGCGTGTACCTGCGCCAGATCATGCAGCAGTATGCCCACGGCAAGCGCGACGACCTGGCCTACCGCCTGGCCCGGCGCAATGCCCACAACGC  >  NC_002947/5646124‑5646410

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