Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,474,324 A→G 8.3% N404N (AAT→AAC)  PP_3964 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,474,3240AG 8.3% 201.0 / 6.8 72N404N (AAT→AAC) PP_3964transposase
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (3/3);  total (38/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 = 9.39e-01

ATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGAGTG  >  NC_002947/4474176‑4474470
                                                                                                                                                    |                                                                                                                                                  
aTAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCa                                                                                                                                                    <  6:263253/149‑1 (MQ=255)
  aaccaggtcgagaaccagatacggccctgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCGTCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTAGCCGTGGCTGATCCGtt                                                                                                                                                  >  4:336503‑M2/54‑149 (MQ=255)
   ACCAGTTCGAGCGGCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCATGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTg                                                                                                                                                 >  4:67075‑M1/1‑149 (MQ=255)
   ACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTg                                                                                                                                                 >  6:124513‑M1/1‑149 (MQ=255)
        ttCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCAc                                                                                                                                            >  4:196486/1‑149 (MQ=255)
           gAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATGCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCa                                                                                                                                          <  1:109165/148‑1 (MQ=255)
               accagatacggccctgggcactcgggcgctcgaactggttGTTGTCGATGGGCACAGCCCCATCGTCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTAGCCGTGGCTGATCCGTTGGGCACCAAgtcg                                                                                                                                     <  3:336503‑M2/109‑2 (MQ=255)
                 ccGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCCCAGCCCCATCATCCCGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGGAATCGAGGGCTTTGGCCGTGGCTGCTCCATTTGGCACCAGGACTCg                                                                                                                                   >  6:194452/1‑149 (MQ=14)
                              gCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATAGGGCACCAGGTCTCGCTGAGCCAACATc                                                                                                                      <  1:20651/149‑1 (MQ=255)
                                    tCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCa                                                                                                                >  5:345036/1‑149 (MQ=255)
                                      tGGTTTTCGACCTGGTTGTTGTCGATGGACACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATg                                                                                                              >  5:329321/1‑149 (MQ=255)
                                               ccacagccaaagccctcgactacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAGTCGAGGGCTTTGGCTGTGGCTGATCCGTTGGGCACCAGGTCACGCTGGGCCAACATCCACTCATGCAGtgtttt                                                                                                     <  8:150200‑M1/100‑1 (MQ=255)
                                               aCCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttttt                                                                                                     <  3:218427/149‑1 (MQ=255)
                                                    gccaaagccctcgactacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAGTCGAGGGCTTTGGCTGTGGCTGATCCGTTGGGCACCAGGTCACGCTGGGCCAACATCCACTCATGCAGTGttttgctga                                                                                                >  7:150200‑M1/45‑145 (MQ=255)
                                                     gtGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGAt                                                                                               <  5:194452‑M2/147‑1 (MQ=255)
                                                         agccctcgactacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAGTCGAGGGCTTTGGCTGTGGCTGATCCGTTGGGCACCAGGTCACGCTGGGCCAACATCCACTCATGCAGtgttt                                                                                                      >  7:259206‑M1/40‑138 (MQ=255)
                                                         agccctcgactacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAGTCGAGGGCTTTGGCTGTGGCTGATCCGTTGGGCACCAGGTCACGCTGGGCCAACATCCACTCATGCAGtgttt                                                                                                      <  8:259206‑M1/99‑1 (MQ=255)
                                                          gtCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACACAGCGTTTGAGGCTGTAATCTAGGGCTTTGGCCGTGGGTGATACATTGGGCCCCAGGTGTCTCTGAGCCCACATCCAGTCATTCAGTTTTTTGATGATCGGCa                                                                                          >  2:166741/1‑149 (MQ=17)
                                                             ctcgattacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCg                                                                                       <  7:261151‑M1/114‑1 (MQ=255)
                                                             ctcgattacagcctgaaacgctgggtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCg                                                                                       <  5:333214‑M1/114‑1 (MQ=255)
                                                              aTGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGc                                                                                      >  1:153536/1‑149 (MQ=255)
                                                                  gCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAAt                                                                                  <  6:299143/149‑1 (MQ=255)
                                                                    acaGCCCCATCATCCAGGTAGCGCGTCACCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                                                                                <  1:166741/149‑1 (MQ=255)
                                                                         cccATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTc                                                                           <  2:153536/149‑1 (MQ=255)
                                                                         cccATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTc                                                                           <  3:196486/149‑1 (MQ=255)
                                                                             tcatcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttt                                                                                                       >  3:85950/1‑117 (MQ=255)
                                                                             tcatcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGtttt                                                                                                       <  4:85950/117‑1 (MQ=255)
                                                                             tcatcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTAt                                                                       <  7:255805/149‑1 (MQ=255)
                                                                             tcatcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTAt                                                                       >  7:15246/1‑149 (MQ=255)
                                                                               atcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTc                                                                     <  8:185938/149‑1 (MQ=255)
                                                                               atcCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTc                                                                     >  5:274968/1‑149 (MQ=255)
                                                                                   aGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGa                                                                                                   <  8:170267/115‑1 (MQ=255)
                                                                                   aGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGc                                                                   >  5:9387/1‑147 (MQ=255)
                                                                                   aGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGc                                                                   >  7:66713/1‑147 (MQ=255)
                                                                                    ggTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAg                                                                >  6:13639/1‑149 (MQ=255)
                                                                                     gtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAgc                                                               >  8:261151‑M1/12‑149 (MQ=255)
                                                                                     gtagcgctgacGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAgc                                                               >  6:333214‑M1/12‑149 (MQ=255)
                                                                                           cgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTg                                                                             >  2:208037/1‑129 (MQ=255)
                                                                                           cgTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTg                                                                             <  1:208037/129‑1 (MQ=255)
                                                                                            gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCtt                                                         >  1:232817/1‑148 (MQ=255)
                                                                                            gTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCtt                                                         >  5:10325/1‑148 (MQ=255)
                                                                                                 CGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                                                                                >  7:59333‑M1/1‑120 (MQ=255)
                                                                                                 CGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTc                                                                                <  8:59333‑M1/120‑1 (MQ=255)
                                                                                                   CTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTc                                                 >  6:345036‑M1/1‑149 (MQ=255)
                                                                                                    TACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCa                                                <  6:123823‑M1/149‑1 (MQ=255)
                                                                                                      CCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGctct                                                          <  6:81261‑M1/137‑1 (MQ=255)
                                                                                                      CCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGctct                                                          >  5:81261‑M1/1‑137 (MQ=255)
                                                                                                      CCCAGCGTTTCAGGCTGTAATCGAAGGCTTTGGCCATGGCTGATCCATTGGGAACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCAGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATg                                              <  8:51167‑M1/149‑1 (MQ=255)
                                                                                                        cAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTc                                            >  7:277808/1‑149 (MQ=255)
                                                                                                         aGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTcc                                           >  6:193411/1‑149 (MQ=255)
                                                                                                          gCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTccc                                          <  8:66713/149‑1 (MQ=255)
                                                                                                                 aGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCt                                    >  7:282379/1‑148 (MQ=255)
                                                                                                                  ggCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGa                                  <  6:10325/149‑1 (MQ=255)
                                                                                                                     tGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGt                               >  2:9999/1‑149 (MQ=255)
                                                                                                                       tAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTc                                                 <  4:36195/129‑1 (MQ=255)
                                                                                                                       tAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTc                                                 >  3:36195/1‑129 (MQ=255)
                                                                                                                       tAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTc                             <  1:9999/149‑1 (MQ=255)
                                                                                                                       tAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTc                             >  1:189772/1‑149 (MQ=255)
                                                                                                                       tAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTc                             <  6:9387/149‑1 (MQ=255)
                                                                                                                          tCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGAc                          <  5:193411/149‑1 (MQ=255)
                                                                                                                            gAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGAcc                         <  2:232817/148‑1 (MQ=255)
                                                                                                                            gAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCt                        >  5:159764/1‑149 (MQ=255)
                                                                                                                                 cTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCGCGCGCCTGACGTTCGACCTCGTAc                   <  4:94864/149‑1 (MQ=255)
                                                                                                                                 cTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTa                    >  3:94864/1‑148 (MQ=255)
                                                                                                                                      gCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGcc              >  7:101147/1‑149 (MQ=255)
                                                                                                                                       ccGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACaa                 <  8:257924/145‑1 (MQ=255)
                                                                                                                                       ccGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACaa                 >  7:257924/1‑145 (MQ=255)
                                                                                                                                           ggCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGa         >  5:252831/1‑149 (MQ=255)
                                                                                                                                              tGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATc       >  4:52409/1‑148 (MQ=255)
                                                                                                                                                aTCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGa     <  8:101147/148‑1 (MQ=255)
                                                                                                                                                aTCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGa     <  3:52409/148‑1 (MQ=255)
                                                                                                                                                aTCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGAg    <  6:159764/149‑1 (MQ=255)
                                                                                                                                                  ccATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGAGTg  >  3:165208/1‑149 (MQ=255)
                                                                                                                                                    aTTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGTGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGa     >  2:232415/1‑144 (MQ=255)
                                                                                                                                                    aTTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGTGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGa     <  1:232415/144‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATAACCAGTTCGAGCGCCCGAGTGCCCATGGCCGTATCTGGTTTTCGACCTGGTTGTTGTCGATGGGCACAGCCCCATCATCCAGGTAGCGCGTCAGCGCTACCCAGCGTTTCAGGCTGTAATCGAGGGCTTTGGCCGTGGCTGATCCATTGGGCACCAGGTCTCGCTGAGCCAACATCCAGTCATGCAGTTTTTTGAGGATCGGCACCGCCAATTCTTGTCGTATTCGCCAGCGCTCTTCATCACTCATGTCCCGCGCCTGACGTTCGACCTCGTACAAGCCGCTGATCGAGTG  >  NC_002947/4474176‑4474470

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