Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,230,786 G→A 8.9% M1M (ATG→ATA) † PP_1969 → molybdenum cofactor biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,230,7860GA 8.9% 191.8 / 7.0 67M1M (ATG→ATA) †PP_1969molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base A (3/3);  total (32/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 = 9.54e-01

GCAAGCGCGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACGGAAAGCGGCTGG  >  NC_002947/2230644‑2230934
                                                                                                                                              |                                                                                                                                                    
gcaagcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAg                                                                                                                                                >  4:23047/1‑149 (MQ=255)
gcaagcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAg                                                                                                                                                >  6:122993/1‑149 (MQ=255)
gcaagcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAg                                                                                                                                                >  2:222955/1‑149 (MQ=255)
  aagcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGt                                                                                                                                              >  3:188171/1‑149 (MQ=255)
   agcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGtt                                                                                                                                             <  8:205705/149‑1 (MQ=255)
    gcgcGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGttt                                                                                                                                            >  2:88836/1‑149 (MQ=255)
      gcgACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAGGttttt                                                                                                                                          <  1:222955/149‑2 (MQ=255)
      gcgACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAGGttttt                                                                                                                                          <  3:23047/149‑2 (MQ=255)
      gcgACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGACTCATAAATCAGGttttt                                                                                                                                          <  5:122993/149‑2 (MQ=255)
         aCCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCTGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCg                                                                                                                                       <  3:84543/149‑1 (MQ=255)
                aCCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCtt                                                                                                                                <  5:206079/149‑1 (MQ=255)
                 cccTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTc                                                                                                                               <  8:234698/149‑1 (MQ=255)
                  ccTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCa                                                                                                                              >  5:187388/1‑149 (MQ=255)
                   cTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCAt                                                                                                                             <  4:188171/149‑1 (MQ=255)
                   cTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCAt                                                                                                                             <  8:155875/149‑1 (MQ=255)
                    tCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATc                                                                                                                            >  2:242620/1‑149 (MQ=255)
                      gAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATc                                                                                                                            <  2:306466/147‑1 (MQ=255)
                         cTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGc                                                                                                                       <  2:71860/149‑1 (MQ=255)
                            aGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCg                                                                                                                     >  2:238451/1‑148 (MQ=255)
                             gCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGa                                                                                                                   <  1:88836/149‑1 (MQ=255)
                             gCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGa                                                                                                                   <  4:348723/149‑1 (MQ=255)
                                 ttAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGtca                                                                                                               >  1:189236/1‑149 (MQ=255)
                                      gAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGtcattcat                                                                                                          <  1:157540/149‑1 (MQ=255)
                                          cGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGAtt                                                                                                      <  6:187388/149‑1 (MQ=255)
                                            gCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGt                                                                                                    <  1:242620/149‑1 (MQ=255)
                                               tcgtGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGa                                                                                                                   >  6:139720/1‑131 (MQ=255)
                                               tcgtGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGa                                                                                                                   <  5:139720/131‑1 (MQ=255)
                                                 gtGAGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCAATGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATc                                                                                               >  6:170219/1‑149 (MQ=255)
                                                  tGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCg                                                                                              >  8:54611/1‑149 (MQ=255)
                                                     gTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTc                                                                                            >  1:313891/1‑148 (MQ=255)
                                                      ttCGCTGCGCGGCATTGCATCCCGAATCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCaa                                                                                          <  7:233076/149‑1 (MQ=255)
                                                        cGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAag                                                                                         <  2:189236/148‑1 (MQ=255)
                                                         gCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGa                                                                                                 <  7:150344/139‑1 (MQ=255)
                                                         gCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGa                                                                                                 >  8:150344/1‑139 (MQ=255)
                                                           tGCGCGGCAGTGCAGCCCGAAGCGTTTGGCAATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAaggcag                                                                                     >  1:30132/1‑149 (MQ=255)
                                                                  cAGTGCAGCCCGAAGTGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCttt                                                                              >  8:258195/1‑149 (MQ=255)
                                                                   aGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTc                                                                             <  7:194155/149‑1 (MQ=255)
                                                                        aGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAAc                                                                        <  1:257954/149‑1 (MQ=255)
                                                                             gAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGC‑AAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTgcgc                                                                  >  7:285532/1‑149 (MQ=255)
                                                                                     tgtcTATCGTTTGGAAATTCGGGATTGCGCATTGGGGCTTGCAAAAACGTGATTTATGTGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCTTTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGc                                                            <  5:164727/145‑1 (MQ=255)
                                                                                     tGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGc                                                            <  1:238451/148‑1 (MQ=255)
                                                                                       gCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCt                                                          >  3:152259/1‑148 (MQ=255)
                                                                                           tCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGa                                                        <  2:30132/146‑1 (MQ=255)
                                                                                           tCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGAc                                                       <  2:313891/147‑1 (MQ=255)
                                                                                             gTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTgcgc                                                                  <  7:334502/134‑1 (MQ=255)
                                                                                             gTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTgcgc                                                                  >  8:334502/1‑134 (MQ=255)
                                                                                             gTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGct                                                   >  6:103142/1‑149 (MQ=255)
                                                                                                 ggAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGctgcct                                               >  5:73384/1‑149 (MQ=255)
                                                                                                     aTTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCaa                                           <  5:170219/149‑1 (MQ=255)
                                                                                                       tCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAAc                                          >  5:62835/1‑148 (MQ=255)
                                                                                                        cGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACt                                         <  7:258195/148‑1 (MQ=255)
                                                                                                              tGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGGGTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGCTCGTCAAGGAAGGCGATTTCGCAACCTGAGCGTCCGACGGTCGGCTGCCTGCAACTCCGACTg                                  >  8:240417/1‑149 (MQ=255)
                                                                                                                   aTTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACgg                                                     >  3:321972/1‑125 (MQ=255)
                                                                                                                   aTTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACgg                                                     <  4:321972/125‑1 (MQ=255)
                                                                                                                   aTTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGctgc                                                 <  4:10618/129‑1 (MQ=255)
                                                                                                                   aTTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGctgc                                                 >  3:10618/1‑129 (MQ=255)
                                                                                                                   aTTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTAcc                              >  8:236565/1‑148 (MQ=255)
                                                                                                                             gCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTAc                    <  7:236565/148‑1 (MQ=255)
                                                                                                                             gCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTAc                    >  2:277813/1‑148 (MQ=255)
                                                                                                                              cAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCt                  <  6:204888/149‑1 (MQ=255)
                                                                                                                              cAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCt                  >  5:204888/1‑149 (MQ=255)
                                                                                                                               aaaaaCCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTg                 >  7:20016/1‑149 (MQ=255)
                                                                                                                                aaaaCCTGATTTAAGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGa                >  7:246150/1‑149 (MQ=255)
                                                                                                                                  aaCCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACg              <  6:73384/149‑1 (MQ=255)
                                                                                                                                     cTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACgg             >  6:26487/1‑147 (MQ=255)
                                                                                                                                     cTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACgg             <  5:26487/147‑1 (MQ=255)
                                                                                                                                          ttATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGACGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCCCTTTAGCAACCTGCGCGTCAGCATGACTGCTGCCTGCAACTAAGCCTGTACCTACTGCGTACCTGACGGAAAGCgg      >  3:344571/1‑149 (MQ=255)
                                                                                                                                              gAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATAGTCAAGGCAGGCGCTTTCGCAACCTGCGAGTCAGCCTGACGGCTGCCTGCAACTAAGCCTGTCCCTACTGCGTACCTTACGGAAAGCGGCtgg  >  7:188665/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCAAGCGCGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACGGAAAGCGGCTGG  >  NC_002947/2230644‑2230934

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