Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,130,198 A→T 8.9% intergenic (‑247/+57) gcvH‑I ← / ← PP_0990 glycine cleavage system protein H/DoxX family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,130,1980AT 8.9% 366.5 / 22.0 132intergenic (‑247/+57)gcvH‑I/PP_0990glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base A (59/60);  new base T (5/7);  total (64/68)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 135 total aligned reads displayed.

TATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTG  >  NC_002947/1130050‑1130345
                                                                                                                                                    |                                                                                                                                                   
tatttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTacgac                                                                                                                                                    <  1:2259216/150‑1 (MQ=255)
   ttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                                                                                 <  1:474504/150‑1 (MQ=255)
        aagaGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaa                                                                                                                                            <  2:677262/147‑1 (MQ=255)
        aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTATCGaa                                                                                                                                            >  1:1864246/1‑150 (MQ=255)
          tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTGCAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                                                                          <  1:2115322/150‑1 (MQ=255)
          tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCTTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                          <  1:4551186/150‑1 (MQ=255)
          tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                                                                          <  2:361933/150‑1 (MQ=255)
          tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                                                                          <  1:2115331/150‑1 (MQ=255)
          tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                          <  1:3017565/150‑1 (MQ=255)
          tAGATTACGTAAATGTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCGTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                          <  2:1628462/150‑1 (MQ=255)
              ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATTGTTTTCGCTACGTCGTAGCGAAAACGAt                                                                                                                                      <  2:243233/150‑1 (MQ=255)
              ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGAt                                                                                                                                      <  1:1652035/150‑1 (MQ=255)
              ttACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCCCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGTGGCTGTGTATCGTTTTCGCTACGACGGAGCGAAAACGAt                                                                                                                                      <  1:563461/150‑1 (MQ=255)
               tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGAAAACGATa                                                                                                                                     <  2:921071/150‑1 (MQ=255)
               tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGAAAACGATa                                                                                                                                     <  2:72330/150‑1 (MQ=255)
               tACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATa                                                                                                                                     <  1:4967796/150‑1 (MQ=255)
                     aaTTTACGGAGCAGCCGTATTTATTGTGTGGAATGATATCGTTACTGGTGCCTCTGAGAGGAAATGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtgtg                                                                                                                               <  2:1048084/150‑1 (MQ=255)
                     aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtgtg                                                                                                                               <  2:4092449/150‑1 (MQ=255)
                           cGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAACCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCt                                                                                                                         <  2:2323681/150‑1 (MQ=255)
                           cGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGCTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCt                                                                                                                         <  1:328330/150‑1 (MQ=255)
                              agaAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCACAAACACCTAGCGAAAACGAAACAAGTTGTGcaagga                                                                                                                      >  1:4301326/1‑145 (MQ=255)
                                 aGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg                                                                                                                   <  1:2983175/150‑1 (MQ=255)
                                 aGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTg                                                                                                                   <  1:878053/150‑1 (MQ=255)
                                   ccGTATTTTTTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                                 <  1:3460236/150‑1 (MQ=255)
                                   ccGTATTTAGGGTGTGGAATGATATCGTCACAGGTGCCTCTGAGAGGAAATGGGTGTGAGCATCAAACCCTGGGGCACCGACGGCGAGGTGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                                 <  1:4840352/150‑1 (MQ=255)
                                     gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                                                                               >  2:1130679/1‑150 (MQ=255)
                                     gTATTTATGGTGTGGAATGATATCGTTCCAGGTGCCTCTTACCGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACTACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                                                                               <  2:4271062/150‑1 (MQ=255)
                                     gTATTGAGTGTGTGGAATGATATCGTGACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAACCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGcac                                                                                                               <  2:5001305/150‑3 (MQ=255)
                                        tttattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCAACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGt                                                                                                            <  2:1684512/150‑1 (MQ=255)
                                          tattGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGc                                                                                                          <  1:2937131/150‑1 (MQ=255)
                                              gtgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTCGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                                                      <  1:1390964/150‑1 (MQ=255)
                                               tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGTAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                                     <  2:511008/150‑1 (MQ=255)
                                               tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                                     <  1:3615231/150‑1 (MQ=255)
                                                 tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGc                                                                                                   >  2:2022173/1‑150 (MQ=255)
                                                 tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGc                                                                                                   <  2:2323532/150‑1 (MQ=255)
                                                 tgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGc                                                                                                   >  2:120163/1‑150 (MQ=255)
                                                   gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCt                                                                                                 >  1:4520504/1‑150 (MQ=255)
                                                   gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCt                                                                                                 >  1:2756903/1‑150 (MQ=255)
                                                    aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCtt                                                                                                <  2:3309565/150‑1 (MQ=255)
                                                      tGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGTGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                                                                              <  1:3630177/150‑1 (MQ=255)
                                                      tGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                                                                              <  2:4866274/150‑1 (MQ=255)
                                                          atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCa                                                                                          >  1:2842221/1‑150 (MQ=255)
                                                          atCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCAACGACGTAGCGAAAACGATACAAGTGGTGCCTGGACAGGCGAGGTGCCAAAGGCCTTCGGtaa                                                                                          >  1:622750/1‑148 (MQ=255)
                                                           tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                                                                         <  2:2110765/150‑1 (MQ=255)
                                                           tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAt                                                                                         <  2:323536/150‑1 (MQ=255)
                                                             gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtct                                                                                       <  2:2167426/150‑1 (MQ=255)
                                                             gTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtct                                                                                       >  1:4720863/1‑150 (MQ=255)
                                                              ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGCCTGGCGAGGTGCCAATGGCCGTCGGTCAtctc                                                                                      >  1:2698823/1‑150 (MQ=255)
                                                              ttACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCAtctc                                                                                      >  1:1006483/1‑150 (MQ=255)
                                                                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACt                                                                                   <  1:1590169/150‑1 (MQ=255)
                                                                 cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACt                                                                                   <  1:1626971/150‑1 (MQ=255)
                                                                    gTGCCTCGGAGAGGAAAAGGGTTTGAGCCTCCAACCCGTGTGCACCGCCGGCGCGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGCAAACGCTACATGGGGTGCCTGGACTGGGGAGGTGCCCAGGGCCGTCGGTCCTCTCAcgtgc                                                                                >  2:2207678/1‑146 (MQ=255)
                                                                     tGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTACCGAAAACGATACATGTGGTGCCTGGACAGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTca                                                                               <  2:2954982/150‑1 (MQ=255)
                                                                      gCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCCCTACGACGTAGCGAAAAACAAACAAGTGGTGCCTGGACTTGCGAGGTGCCAATGGCCATCGGTCATCTCACATAcag                                                                              >  1:1500080/1‑150 (MQ=255)
                                                                        ctctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagca                                                                            <  1:1355813/150‑1 (MQ=255)
                                                                         tctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagcag                                                                           >  1:116190/1‑150 (MQ=255)
                                                                         tctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagcag                                                                           >  2:4879124/1‑150 (MQ=255)
                                                                         tctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagcag                                                                           >  2:4451555/1‑150 (MQ=255)
                                                                         tctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagcag                                                                           >  2:4451409/1‑150 (MQ=255)
                                                                          ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagcagg                                                                          >  1:3231594/1‑149 (MQ=255)
                                                                           tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                          >  2:4337660/1‑85 (MQ=255)
                                                                           tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                          <  1:4337779/85‑1 (MQ=255)
                                                                           tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgat                                                                         >  2:4407589/1‑150 (MQ=255)
                                                                            gagaGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAgatg                                                                        >  1:3392845/1‑150 (MQ=255)
                                                                               agGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc                                                                     >  1:823410/1‑150 (MQ=255)
                                                                                  acaaGGGTTTGATGCTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATACAg                                                                  >  1:4292577/3‑150 (MQ=255)
                                                                                  aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAg                                                                  >  2:1974836/1‑150 (MQ=255)
                                                                                  aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAg                                                                  >  2:5030137/1‑150 (MQ=255)
                                                                                    aaGGGTTTGATGCTCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGAGACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCg                                                                >  2:3737519/1‑150 (MQ=255)
                                                                                    aaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCg                                                                >  1:4400228/1‑150 (MQ=255)
                                                                                         tttGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAc                                                           <  2:3304510/150‑1 (MQ=255)
                                                                                              gCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCCATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGcc                                                      >  1:1591350/1‑150 (MQ=255)
                                                                                              gCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGcc                                                      >  1:3488103/1‑150 (MQ=255)
                                                                                               cATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCg                                                     >  2:2699968/1‑150 (MQ=255)
                                                                                                 tCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGGGGTGCCTGGACTGGCGAGGTGCCCATGGCCGTCGGTCATCTCACTTTCAGCAGCTGATCCAGCGACAACGCGCCGgc                                                   >  1:4447727/1‑150 (MQ=255)
                                                                                                   aaaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGgcgc                                                 >  2:694814/1‑150 (MQ=255)
                                                                                                    aaCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGACAACGCGCCGGCGcc                                                >  2:635678/1‑150 (MQ=255)
                                                                                                     aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCt                                               >  2:4333437/1‑150 (MQ=255)
                                                                                                      cccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTc                                              <  2:4885045/150‑1 (MQ=255)
                                                                                                      cccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTc                                              >  2:4139069/1‑150 (MQ=255)
                                                                                                       ccTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCccgcgccgcg                                             >  1:3332948/1‑150 (MQ=255)
                                                                                                               aCCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCa                                     <  2:4953536/150‑1 (MQ=255)
                                                                                                               aCCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACCACGTAGCCAAAACGATACATGTTGTGCCTGGACTGGCGAGGTGCCCATGGCCGTCGGTCATAGCACATTCAGCAGATGATCCACCGAAAAAGCGCCGGCGCCACGGGCAAccc                                     >  1:3899577/1‑149 (MQ=255)
                                                                                                                           gaggCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAaga                         <  2:823410/150‑1 (MQ=255)
                                                                                                                                    tATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGccc                >  1:2848888/1‑150 (MQ=255)
                                                                                                                                    tATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGccc                >  1:247820/1‑150 (MQ=255)
                                                                                                                                     aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCca               <  1:2964423/150‑1 (MQ=255)
                                                                                                                                     aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCca               >  2:1785750/1‑150 (MQ=255)
                                                                                                                                      tCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                                                      <  2:3514606/62‑1 (MQ=255)
                                                                                                                                      tCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCgg                                                    >  2:624898/1‑112 (MQ=255)
                                                                                                                                      tCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcag              <  2:2783770/150‑1 (MQ=255)
                                                                                                                                        gTTTTCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggt            >  1:887174/1‑150 (MQ=255)
                                                                                                                                        gTTTTCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggt            >  1:3383701/1‑150 (MQ=255)
                                                                                                                                         ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtc           >  1:1331016/1‑150 (MQ=255)
                                                                                                                                           ttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcag         >  1:2218905/1‑150 (MQ=255)
                                                                                                                                           ttCGCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcag         >  1:3158574/1‑150 (MQ=255)
                                                                                                                                            tCGCTACGCCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcagg        <  1:3390394/150‑1 (MQ=255)
                                                                                                                                                tACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtggg    >  1:2000838/1‑150 (MQ=255)
                                                                                                                                                 acgacgTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtgggt   <  1:56319/150‑1 (MQ=255)
                                                                                                                                                  cgacgTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtgggtg  <  1:1130679/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGGTG  >  NC_002947/1130050‑1130345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: