Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 516,282 A→T 5.6% L238Q (CTG→CAG)  PP_0429 ← short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947516,2820AT 5.6% 217.9 / 3.6 71L238Q (CTG→CAG) PP_0429short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (2/2);  total (36/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 = 7.53e-01

GGAATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCG  >  NC_002947/516133‑516424
                                                                                                                                                     |                                                                                                                                              
ggAATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgca                                                                                                                                                >  1:1342969/1‑150 (MQ=255)
 gAATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcag                                                                                                                                               >  1:1208343/1‑150 (MQ=255)
   aTAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcagca                                                                                                                                             <  2:348308/150‑1 (MQ=255)
   aTAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGgcagca                                                                                                                                             <  2:387056/150‑1 (MQ=255)
        tCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCcgc                                                                                                                                        >  1:626038/1‑150 (MQ=255)
           tCCACAGCACTTGCTCCGTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                     >  2:285603/1‑150 (MQ=255)
             cacaGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCACCGCGGCTGCTGCCGCGCTgg                                                                                                                                   <  2:626038/150‑2 (MQ=255)
                         tCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACgccagcca                                                                                                                       <  2:1208343/150‑1 (MQ=255)
                           cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                     <  1:224306/150‑1 (MQ=255)
                           cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                     >  1:40333/1‑150 (MQ=255)
                           cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCcagccagc                                                                                                                     <  1:946903/150‑1 (MQ=255)
                           cATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGAGGTGAGAACGCcagccagc                                                                                                                     >  1:1064090/1‑150 (MQ=255)
                            aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                <  1:1272580/138‑1 (MQ=255)
                            aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                >  2:1272580/1‑138 (MQ=255)
                            aTGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCg                                                                                                                    <  2:40333/150‑1 (MQ=255)
                                 gAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTgc                                                                                                               <  1:503590/150‑1 (MQ=255)
                                    ggTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCgcc                                                                                                            >  1:429963/1‑150 (MQ=255)
                                       cGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGAGGTGAGAACGCCAGCCAGCGGTTGCgccggc                                                                                                         <  2:1064090/150‑1 (MQ=255)
                                          gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGtt                                                                                                                 >  2:1062619/1‑139 (MQ=255)
                                          gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGtt                                                                                                                 <  1:1062619/139‑1 (MQ=255)
                                          gCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGccggccgg                                                                                                      <  2:579658/150‑1 (MQ=255)
                                                     tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                >  1:1370362/1‑113 (MQ=255)
                                                     tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGaa                                                                                                                                <  2:1370362/113‑1 (MQ=255)
                                                     tCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGt                                                                                           >  2:1228826/1‑150 (MQ=255)
                                                      cAGCCTTGAGCCGCTCTCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGGCGGGCATGATCAGTg                                                                                          <  1:49227/150‑1 (MQ=255)
                                                              aGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                     >  1:315095/1‑99 (MQ=255)
                                                              aGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGt                                                                                                                                     <  2:315095/99‑1 (MQ=255)
                                                                ccGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                <  1:123879/150‑1 (MQ=255)
                                                                 cGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTg                                                                               >  2:61306/1‑150 (MQ=255)
                                                                   cTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTc                                                                             >  1:1109669/1‑150 (MQ=255)
                                                                    tCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTcc                                                                            <  1:1159740/150‑1 (MQ=255)
                                                                       cccAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGt                                                                         <  1:290711/150‑1 (MQ=255)
                                                                        ccAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTg                                                                        >  1:1397112/1‑150 (MQ=255)
                                                                           gTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCAt                                                                                                  <  2:1181941/121‑1 (MQ=255)
                                                                           gTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCAt                                                                                                  >  1:1181941/1‑121 (MQ=255)
                                                                              gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                <  2:604575/136‑1 (MQ=255)
                                                                              gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCtt                                                                                >  1:604575/1‑136 (MQ=255)
                                                                              gCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCa                                                                  <  1:61306/150‑1 (MQ=255)
                                                                               cTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCAc                                                                 <  1:605272/150‑1 (MQ=255)
                                                                                 aTTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACgg                                                               >  1:637334/1‑150 (MQ=255)
                                                                                  ttCCATGCCTGAGGGGCAATCGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCGGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGGGCGGTTCTTGTCCAGTGCCCGCACggc                                                              >  2:158719/1‑149 (MQ=255)
                                                                                  ttCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGgt                                                              >  2:951100/1‑150 (MQ=255)
                                                                                     cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAAc                                                                                                                               <  1:1169744/82‑1 (MQ=255)
                                                                                     cATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAAc                                                                                                                               >  2:1169744/1‑82 (MQ=255)
                                                                                         ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgc                                                       <  2:637334/150‑1 (MQ=255)
                                                                                         ccTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgc                                                       <  2:1342969/150‑1 (MQ=255)
                                                                                          cTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCGAGGGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGAATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTCCAggg                                                      >  2:145367/1‑148 (MQ=255)
                                                                                           tGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                     >  2:787592/1‑150 (MQ=255)
                                                                                             aGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCa                                                   <  2:429963/150‑1 (MQ=255)
                                                                                                         ggcggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCGGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                     >  1:529695/1‑136 (MQ=255)
                                                                                                         ggcggTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCGGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                     <  2:529695/136‑1 (MQ=255)
                                                                                                               ggATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGc                                 >  2:284994/1‑150 (MQ=255)
                                                                                                                     ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCgg                                                                                     <  2:857316/92‑1 (MQ=255)
                                                                                                                     ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCgg                                                                                     >  1:857316/1‑92 (MQ=255)
                                                                                                                     ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGcgcg                           >  2:1159099/1‑150 (MQ=255)
                                                                                                                     ccccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGcgcg                           >  1:314584/1‑150 (MQ=255)
                                                                                                                       ccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGc                         >  2:108156/1‑150 (MQ=255)
                                                                                                                       ccGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGc                         <  2:1397112/150‑1 (MQ=255)
                                                                                                                         ggCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                       <  1:1228826/150‑1 (MQ=255)
                                                                                                                             aGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCa                   <  1:105031/150‑1 (MQ=255)
                                                                                                                              gCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCAt                  >  1:487295/1‑150 (MQ=255)
                                                                                                                                ttGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCg                >  1:229641/1‑150 (MQ=255)
                                                                                                                                   cgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGc             <  1:108156/150‑1 (MQ=255)
                                                                                                                                   cgcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGc             >  1:1198703/1‑150 (MQ=255)
                                                                                                                                    gcACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGcc            >  1:1139855/1‑150 (MQ=255)
                                                                                                                                      aCCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCTGCGGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgc          <  2:550315/150‑1 (MQ=255)
                                                                                                                                      aCCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCcgc          >  1:38927/1‑150 (MQ=255)
                                                                                                                                          gccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCaa      <  1:145367/150‑1 (MQ=255)
                                                                                                                                          gccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCaa      >  2:1039302/1‑150 (MQ=255)
                                                                                                                                           ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                     <  1:972142/102‑1 (MQ=255)
                                                                                                                                           ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAgcgc                                                     >  2:972142/1‑102 (MQ=255)
                                                                                                                                           ccagcGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAAt     <  1:951100/150‑1 (MQ=255)
                                                                                                                                             accgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                       <  2:216515/128‑1 (MQ=255)
                                                                                                                                             accgcgGCTGCTGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCtt                       >  1:216515/3‑130 (MQ=255)
                                                                                                                                              gcgcGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATgcg  >  1:1216808/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGAATAGATCGTCCACAGCACTTGCTCCATGGTGAGGGTCGGGCCGAGTGTACTCAGCCGTGAGCCGCTCGCCCAGTGGCTATTCCATGCCTGAGGGGCAATAGCGGCGGTGGATGGCCCCGGCCAGCTTGCGCACCAGCCAGCGCGGCAGCAGCCGCGGTGCGAACGCCAGCCAGCGGTTGCGCCGGCCGGGCATGATCAGTGCGCGGTTCTTGTCCAGTGCCCGCACGGTGTACAGCGCCACTTCCTCGGGGCTGAGGCAGCGCGGCTTGCCATCGAGCCGCGCAATGCG  >  NC_002947/516133‑516424

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