Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,653,161 A→G 6.3% L1203P (CTG→CCG)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,653,1610AG 6.3% 182.8 / 2.7 63L1203P (CTG→CCG) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base G (2/2);  total (29/34)
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.37e-01

CGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCATGCCGCCACCCAGG  >  NC_002947/1653016‑1653304
                                                                                                                                                 |                                                                                                                                               
cgAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGGTGCGggtcaggg                                                                                                                                              >  4:6580/1‑149 (MQ=255)
 gagaTGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTggtcagggt                                                                                                                                             <  4:166745/149‑1 (MQ=255)
      gTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAgtg                                                                                                                                        >  4:12775/1‑149 (MQ=255)
         gAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGc                                                                                                                                     >  3:251430/1‑149 (MQ=255)
         gAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGc                                                                                                                                     >  7:105864/1‑149 (MQ=255)
          aGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGcc                                                                                                                                    <  1:372422/149‑1 (MQ=255)
               gtcgcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGc                                                                                                                               <  6:234543/149‑1 (MQ=255)
                tcgcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCt                                                                                                                              <  7:158931/149‑1 (MQ=255)
                tcgcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCt                                                                                                                              <  3:430218/149‑1 (MQ=255)
                  gcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTgg                                                                                                                            <  1:310925/149‑1 (MQ=255)
                  gcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTgg                                                                                                                            >  2:310925/1‑149 (MQ=255)
                  gcgtcgcTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTg                                                                                                                             <  8:160691/148‑1 (MQ=255)
                  gcgtcgcTCCTGGCTGGAAGAGGATTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCTCTGCTGGTCAGGGTAAGTGTGCCGTTGCTgg                                                                                                                            >  3:351497/1‑149 (MQ=255)
                        cTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACt                                                                                                                      >  8:131003/1‑149 (MQ=255)
                           ctggctggAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTgg                                                                                                                    <  5:99936/148‑1 (MQ=255)
                           ctggctggAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTgg                                                                                                                    <  6:208071/148‑1 (MQ=255)
                           ctggctggAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTgg                                                                                                                    <  7:131003/148‑1 (MQ=255)
                           ctggctggAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTg                                                                                                                     <  4:275791/147‑1 (MQ=255)
                              gctggAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCACTGTGCCGTTGCTGGCAAACTGGGCAt                                                                                                                <  1:232803/149‑1 (MQ=255)
                                     gagGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGc                                                                                                         >  6:217135/1‑149 (MQ=255)
                                         cTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGa                                                                                                     >  7:372063/1‑149 (MQ=255)
                                            cAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGAc                                                                                                    <  1:69841/147‑1 (MQ=255)
                                             aaCCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGt                                                                                                  <  5:331380/148‑1 (MQ=255)
                                              aCCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGc                                                                                                <  8:189582/149‑1 (MQ=255)
                                                          ccATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGAAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTg                                                                                    >  3:352878/1‑149 (MQ=255)
                                                                 aGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGt                                                                                                  <  2:187462/128‑1 (MQ=255)
                                                                 aGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGt                                                                                                  >  1:187462/1‑128 (MQ=255)
                                                                    cGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTgg                                                                          >  2:14010/1‑149 (MQ=255)
                                                                      tCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCt                                                                        >  3:379037/1‑149 (MQ=255)
                                                                      tCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCt                                                                        >  7:234809/1‑149 (MQ=255)
                                                                            gTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTg                                                                  <  1:33806/149‑1 (MQ=255)
                                                                             tGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCCTTGCTGGCAAACTGGGCATCGACCGCCTGACGTCCATGTTCGTCCTGCTTGCCGTGGCTGAAGTgg                                                                 >  2:284326/1‑149 (MQ=255)
                                                                                  tGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGgccgcc                                                            >  6:52843/1‑149 (MQ=255)
                                                                                    ttCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCGGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTa                                                          >  2:118470/1‑149 (MQ=255)
                                                                                    ttCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCGGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTa                                                          >  2:112673/1‑149 (MQ=255)
                                                                                    ttCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTa                                                          >  1:257253/1‑149 (MQ=255)
                                                                                       ggCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATgg                                                       <  3:12775/149‑1 (MQ=255)
                                                                                       ctgcaaGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATgg                                                       >  2:407776/6‑149 (MQ=255)
                                                                                         cctgcaccTGCAAGGCCTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCt                                                     >  4:394897/7‑149 (MQ=255)
                                                                                         ccGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCt                                                     <  4:357509/149‑1 (MQ=255)
                                                                                         ccGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCt                                                     <  8:234809/149‑1 (MQ=255)
                                                                                         ccGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGACCGCCTATGGCt                                                     >  8:263170/1‑149 (MQ=255)
                                                                                                cTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCt                                              <  1:407776/149‑1 (MQ=255)
                                                                                                  gCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTgg                                            <  1:14010/149‑1 (MQ=255)
                                                                                                   cAGGCCTTGCAGGTGCAGTGCGATGTCCTCGCGGGCGCTGCTGGGCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTgg                                            <  8:376922/148‑1 (MQ=255)
                                                                                                   cAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCGGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTgg                                            <  1:112673/148‑1 (MQ=255)
                                                                                                   cAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCGGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTgg                                            <  1:118470/148‑1 (MQ=255)
                                                                                                      gCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTGCTGCTTGCCGTGGCTGAAGTGgcc                                                               <  8:152540/126‑1 (MQ=255)
                                                                                                      gCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGgcc                                                               >  7:152540/1‑126 (MQ=255)
                                                                                                      gCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCTTGGGTCt                                        >  1:380050/1‑149 (MQ=255)
                                                                                                                tgcaggGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCtt                              <  8:372063/149‑1 (MQ=255)
                                                                                                                    gggCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGgcc                          <  4:379037/149‑1 (MQ=255)
                                                                                                                      gCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGACCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGgccgc                        <  7:263170/149‑1 (MQ=255)
                                                                                                                       cGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAATCTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGgccgcc                       <  6:132792/149‑1 (MQ=255)
                                                                                                                             ccTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTc                  >  6:24821/1‑148 (MQ=255)
                                                                                                                              cTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTgg                                                                          <  7:374879/91‑1 (MQ=255)
                                                                                                                              cTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTgg                                                                          >  8:374879/1‑91 (MQ=255)
                                                                                                                                 gcgGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCAt               <  2:257253/147‑1 (MQ=255)
                                                                                                                                 gcgGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCATGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCAt               <  2:365740/147‑1 (MQ=255)
                                                                                                                                   gggCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCATgccg           <  3:394897/149‑1 (MQ=255)
                                                                                                                                    ggCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCATgccgc          >  4:380098/1‑149 (MQ=255)
                                                                                                                                         tgctgGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCATGCCGCCAcc      >  6:315291/1‑148 (MQ=255)
                                                                                                                                           ctgGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTgg                            >  7:49073/1‑124 (MQ=255)
                                                                                                                                           ctgGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTgg                            <  8:49073/124‑1 (MQ=255)
                                                                                                                                            tgGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCCTGCCGCCACCCAgg  >  7:289727/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGAGATGTCGAGGTTGTCGCGTCGCTCCTGGCTGGAAGAGGCTTCAACCAGCATGCCGCCATTCGAGGCGTCGAGGGTGATCTGTTCGGCCGAGGCCTGCAGGCCTTGCAGGTGCAGGGCGATGTCCTCGCGGGCGCTGCTGGTCAGGGTCAGTGTGCCGTTGCTGGCAAACTGGGCATCGACCGCCTGACGTGCATGTTCGTCCTGCTTGCCGTGGCTGAAGTGGCCGCCTATGGCTCCCCCCTGGGTCTGGCCAGCCTTGGCCGCCAGCTCCATGCCGCCACCCAGG  >  NC_002947/1653016‑1653304

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