Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,228 A→T 11.0% intergenic (+56/+75) prfC → / ← potF‑I peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,2280AT11.0% 205.5 / 13.6 76intergenic (+56/+75)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base T (4/4);  total (36/40)
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.54e-01

CTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGC  >  NC_002947/1013085‑1013375
                                                                                                                                                |                                                                                                                                                   
ctacctGGCGCCGACCCGGGTCAACCTGGCGCTGGTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                                                <  5:64976/149‑1 (MQ=255)
ctacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                                                <  8:147755/149‑1 (MQ=255)
ctacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                                                >  7:212472/1‑149 (MQ=255)
ctacctGGCGCCGACCCGGGTCAACCTCGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                                                <  4:175291/149‑1 (MQ=255)
ccacctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAc                                                                                                                                                <  8:42047/147‑1 (MQ=255)
    ctGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAcagac                                                                                                                                            <  4:348757/149‑5 (MQ=255)
     tGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGAc                                                                                                                                           >  7:208074/1‑149 (MQ=255)
            gACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAgg                                                                                                                                    <  6:289496/149‑1 (MQ=255)
            gACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAgg                                                                                                                                    <  6:167040/149‑1 (MQ=255)
             acgcGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                     >  4:6110/4‑131 (MQ=255)
             aCCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCTCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAgg                                                                                                                                    >  1:258817/1‑149 (MQ=255)
             aCCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGa                                                                                                                                                     <  3:6110/131‑1 (MQ=255)
               ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAcc                                                                                                                                               >  5:107957/1‑135 (MQ=255)
               ccGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCAcc                                                                                                                                               <  6:107957/135‑1 (MQ=255)
                 gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTTCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTctgg                                                                                                                               <  3:255168/149‑4 (MQ=37)
                 gggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtg                                                                                                                               >  8:75973/1‑149 (MQ=255)
                  ggTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                              >  4:262637/1‑149 (MQ=255)
                      aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGGCATGTGATa                                                                                                                          >  4:282349/1‑149 (MQ=255)
                       aCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                         >  3:43334/1‑149 (MQ=255)
                       aCCTGGCGCTGATGGAACAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATAGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                         >  3:38413/1‑149 (MQ=255)
                         cTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCt                                                                                                                       <  2:119041/149‑1 (MQ=255)
                            gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGTGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                    <  8:50575/149‑1 (MQ=25)
                            gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                     <  8:212472/148‑1 (MQ=255)
                             cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                   <  1:114241/149‑1 (MQ=255)
                             cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                   <  7:75973/149‑1 (MQ=255)
                             cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                   <  4:43334/149‑1 (MQ=255)
                             cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTgg                                                                                                                   <  4:38413/149‑1 (MQ=255)
                              gcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTggg                                                                                                                  >  5:316152/1‑149 (MQ=255)
                               cTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGa                                                                                                                 <  3:328225/149‑1 (MQ=21)
                                 gATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGGACGCCCGCACCCACACGAATCCCCACACCGGGAGGTCATGTGAAGCCGGTGGgcgc                                                                                                               >  4:154753/1‑149 (MQ=38)
                                  aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAg                                                                                                                <  2:35566/147‑1 (MQ=255)
                                    ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                    >  6:325073/1‑141 (MQ=21)
                                    ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg                                                                                                                    <  5:325073/141‑1 (MQ=21)
                                    ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                         >  7:189792/1‑136 (MQ=255)
                                    ggAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc                                                                                                                         <  8:189792/136‑1 (MQ=255)
                                      aaGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCCCCCCCCCACCCACAAGAAACACCAACCCAGGCTGTGAAGTTATAACTGTGGGAGCggggg                                                                                                          >  1:141774/1‑147 (MQ=255)
                                           cgcTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGccc                                                                                                     <  6:316152/149‑1 (MQ=255)
                                                   cGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                              <  6:195455/117‑1 (MQ=255)
                                                   cGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                              >  7:16276/1‑149 (MQ=255)
                                                        aTCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGc                                                                                         >  5:179852/1‑149 (MQ=255)
                                                          cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGc                                                                                         >  7:186714/1‑147 (MQ=255)
                                                             gTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAggc                                                                                    <  3:282349/149‑1 (MQ=255)
                                                                  gcgcTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGccc                                                                               <  8:208074/149‑1 (MQ=255)
                                                                        cccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCAttt                                                                          >  2:226997/1‑148 (MQ=255)
                                                                         ccGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCa                                                                        <  2:258817/149‑1 (MQ=255)
                                                                           gcgAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                     <  3:106717/102‑1 (MQ=255)
                                                                           gcgAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCtgt                                                                                                                     >  4:106717/1‑102 (MQ=255)
                                                                                   cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCt                                                              >  5:94894/1‑149 (MQ=255)
                                                                                   cATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCGGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTAAATTTGCGGcg                                                              >  4:344337/1‑148 (MQ=255)
                                                                                       aaGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGt                                                                                                         >  8:336354/1‑102 (MQ=255)
                                                                                       aaGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGt                                                                                                         <  7:336354/102‑1 (MQ=255)
                                                                                        aGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGa                                                          <  1:226997/148‑1 (MQ=255)
                                                                                         ggCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACAGGGAAAACCAGCCCGGGAGGTGATGTGATAACGGTGGGGGCGGGGGGGGCCGCGAAGGGCCGCGAGGCGGGCCCATTTGAATTTACGGCGTTgtg                                                         >  4:243063/1‑146 (MQ=21)
                                                                                             tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCAtgt                                                                                                                              >  6:379436/1‑75 (MQ=255)
                                                                                                    ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGccc                                                                                                     >  6:194241/1‑93 (MQ=255)
                                                                                                    ggggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGccc                                                                                                     <  5:194241/93‑1 (MQ=255)
                                                                                                        ttgctTTTCAGCCCTTCGAGGCCCCGGCCCATGACACTTCATTCACGACACATGGTGGTCTTTTGATACGTGTGGGTGCGGGTGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCGGGCTTTGATACTGGTCCAGATTctc                                         <  6:19381/148‑3 (MQ=255)
                                                                                                        ctgctTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTcgc                                         <  3:262637/149‑1 (MQ=255)
                                                                                                                 aGCCCTTCGCGGGCACGCCCGCTCCCACCGGAATCACCAGAACGGGAGGTGAAGTGATACCGGGGGGAGCGGGCGGGCCCGCGGAGGGCCGCTAGGGGGGCCCCTTTCAATTTACGGCTTTGATAATGGGGCAGATTAGCGGACGAAt                                 >  6:200377/1‑148 (MQ=21)
                                                                                                                   cccTTCGCGGCCACCCCCGCTCCCACAGCAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACg                              <  3:208640/149‑1 (MQ=255)
                                                                                                                    ccTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACgg                             <  8:186714/149‑1 (MQ=255)
                                                                                                                     cTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGt                            >  7:252787/1‑149 (MQ=255)
                                                                                                                           gggCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCt                      >  6:98082/1‑149 (MQ=255)
                                                                                                                               aCGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGc                  >  2:61112/1‑149 (MQ=255)
                                                                                                                                       tCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCAGCGAAGGGCCGCAAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGAACGAATACGCTCGATCTTCGCCGGCATcc          >  5:298854/1‑148 (MQ=18)
                                                                                                                                         ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGc         >  6:15537/1‑148 (MQ=21)
                                                                                                                                         ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAAttt                                                                    >  5:101851/1‑89 (MQ=255)
                                                                                                                                         ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAAttt                                                                    <  6:101851/89‑1 (MQ=255)
                                                                                                                                         ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc        >  5:116299/1‑149 (MQ=255)
                                                                                                                                         ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc        <  6:116299/149‑1 (MQ=255)
                                                                                                                                         ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc        >  8:100357/1‑149 (MQ=255)
                                                                                                                                          cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt       <  6:94894/149‑1 (MQ=255)
                                                                                                                                          cacaGGAAGCACCAGACCTGGCGGTCATGTGATACCTGTGGTAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCt       <  3:243063/149‑1 (MQ=255)
                                                                                                                                           acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc      >  2:52204/1‑149 (MQ=25)
                                                                                                                                             aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCa    <  6:179852/149‑1 (MQ=255)
                                                                                                                                               gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGa                         >  3:142702/1‑126 (MQ=255)
                                                                                                                                               gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGa                         <  4:142702/126‑1 (MQ=255)
                                                                                                                                               gAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgc  <  8:252787/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTACCTGGCGCCGACCCGGGTCAACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGC  >  NC_002947/1013085‑1013375

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