Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,230,788 G→A 8.8% S2N (AGT→AAT)  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,7880GA 8.8% 192.8 / 5.9 68S2N (AGT→AAT) PP_1969molybdenum cofactor biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (3/3);  total (33/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.80e-01

GCAAGCGCGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACGGAAAGCGGCTGGTG  >  NC_002947/2230644‑2230936
                                                                                                                                                |                                                                                                                                                    
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                                                                    >  8:334502/1‑134 (MQ=255)
                                                                                             gTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTgcgc                                                                    <  7:334502/134‑1 (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)
                                                                                                                                                gTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACGGAAAGCGGCtggtg  >  4:52848/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCAAGCGCGACCCGGAACCCTCGAACTGAGCATTTAGCGAAACGGCGTCGTGCGTTCGCTGCGCGGCAGTGCAGCCCGAAGCGTTTGGCTATCGTTTGGAAATTCGGCATTGCGCATTGGGGCTTGCAAAAACCTGATTTATGAGTCAGGTTTTGCCGGCCCCTTCATCCATGCCGGAGTCATTCATGATTGTCGATCGTCAAGGCAGGCGCTTTCGCAACCTGCGCGTCAGCCTGACGGCTGCCTGCAACTACGCCTGTACCTACTGCGTACCTGACGGAAAGCGGCTGGTG  >  NC_002947/2230644‑2230936

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