Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,189,213 C→A 5.6% P93P (CCC→CCA PP_5490 → site‑specific DNA recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,189,2130CA 5.6% 896.4 / 28.7 265P93P (CCC→CCAPP_5490site‑specific DNA recombinase
Reads supporting (aligned to +/- strand):  ref base C (123/127);  new base A (3/12);  total (126/139)
Fisher's exact test for biased strand distribution p-value = 3.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 273 total aligned reads displayed.

GAGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCACG  >  NC_002947/2189064‑2189362
                                                                                                                                                     |                                                                                                                                                     
 agTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCt                                                                                                                                                      <  2:5958525/150‑1 (MQ=255)
 agTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCt                                                                                                                                                      <  1:4872200/150‑1 (MQ=255)
  gTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTg                                                                                                                                                     <  2:6202203/150‑1 (MQ=255)
       tgaGGCTGCCTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCa                                                                                                                                                <  2:530909/150‑1 (MQ=255)
         aGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCaga                                                                                                                                              <  2:3787258/150‑3 (MQ=255)
         aGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCaga                                                                                                                                              <  1:4569002/150‑3 (MQ=255)
           gctgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAacac                                                                                                                                            >  1:1509539/1‑150 (MQ=255)
             tgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAAcacca                                                                                                                                          >  2:6844702/1‑150 (MQ=255)
                ttGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTa                                                                                                                                       <  1:5769686/150‑1 (MQ=255)
                     catcTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgct                                                                                                                                  <  2:213797/150‑1 (MQ=255)
                       tctcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgc                                                                                                                                <  1:5483240/150‑1 (MQ=255)
                       tctcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgc                                                                                                                                <  2:569897/150‑1 (MQ=255)
                            ggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTc                                                                                                                           <  1:4564879/150‑1 (MQ=255)
                             gccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCa                                                                                                                          >  2:4912810/1‑150 (MQ=255)
                              ccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAg                                                                                                                         >  2:17070/1‑150 (MQ=255)
                                aggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgat                                                                                                                       <  1:3727175/150‑1 (MQ=255)
                                 ggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgatg                                                                                                                      >  2:2470612/1‑150 (MQ=255)
                                      tgtggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAAt                                                                                                                 <  1:5784087/150‑1 (MQ=255)
                                      tgtggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAAt                                                                                                                 <  2:1091044/150‑1 (MQ=255)
                                        tggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGa                                                                                                               <  1:2774523/150‑1 (MQ=255)
                                          gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACg                                                                                                             >  2:5601290/1‑150 (MQ=255)
                                          gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACg                                                                                                             >  1:2586832/1‑150 (MQ=255)
                                              tggtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGatat                                                                                                         >  1:195504/1‑150 (MQ=255)
                                                gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                       >  2:4072538/1‑150 (MQ=255)
                                                gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                       <  1:3818069/150‑1 (MQ=255)
                                                gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                       <  2:1264704/150‑1 (MQ=255)
                                                gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                       >  2:2555721/1‑150 (MQ=255)
                                                 tgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTg                                                                                                      <  2:6545256/150‑1 (MQ=255)
                                                   tgtcgtcTTGATGGGCTCGCTCGGTCTACTATTGATTTGCTGTCTATTGCAGAGCAGCTGTATCGCTTGGGCGCTTGGTTGGAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACTGCTGCTGGTCAGATGATAATGACGATATTTGCa                                                                                                    <  1:5924229/149‑1 (MQ=255)
                                                   aGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCa                                                                                                    >  2:727409/1‑150 (MQ=255)
                                                   aGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCa                                                                                                    >  2:5233622/1‑150 (MQ=255)
                                                    gtcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAg                                                                                                   <  1:681910/150‑1 (MQ=255)
                                                    gtcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAg                                                                                                   <  1:2107297/150‑1 (MQ=255)
                                                    gtcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAg                                                                                                   <  2:2169630/150‑1 (MQ=255)
                                                     tcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATCTTTGCAgg                                                                                                  >  1:5557583/1‑150 (MQ=255)
                                                     tcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAgg                                                                                                  >  1:1006113/1‑150 (MQ=255)
                                                     tcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAgg                                                                                                  >  1:3375859/1‑150 (MQ=255)
                                                        tcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCAt                                                                                               <  1:6373510/150‑1 (MQ=255)
                                                         cTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATc                                                                                              <  1:80004/150‑1 (MQ=255)
                                                          ttGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCg                                                                                             >  1:4857518/1‑150 (MQ=255)
                                                          ttGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCg                                                                                             >  1:74399/1‑150 (MQ=255)
                                                             aTCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTa                                                                                          >  1:4304493/1‑150 (MQ=255)
                                                                    cGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGa                                                                                   <  2:4304493/150‑1 (MQ=255)
                                                                     gCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAg                                                                                  <  2:4553905/150‑1 (MQ=255)
                                                                      cACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAgc                                                                                 <  2:5158878/150‑1 (MQ=255)
                                                                         ggTCGACAATTGCTTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCa                                                                              <  1:5439861/150‑1 (MQ=255)
                                                                            cGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGcc                                                                           <  2:1344460/150‑1 (MQ=255)
                                                                             gACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCt                                                                          >  1:6014002/1‑150 (MQ=255)
                                                                             gACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCACAGCCt                                                                          >  1:127977/1‑150 (MQ=255)
                                                                                aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                       >  1:6271733/1‑150 (MQ=255)
                                                                                aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                       <  1:3745145/150‑1 (MQ=255)
                                                                                aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                       >  1:2716733/1‑150 (MQ=255)
                                                                                 aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                      <  1:2570375/150‑1 (MQ=255)
                                                                                 aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                      <  2:2674457/150‑1 (MQ=255)
                                                                                 aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                      <  1:2570446/150‑1 (MQ=255)
                                                                                 aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                      <  1:6701010/150‑1 (MQ=255)
                                                                                    gATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTa                                                                   <  2:3164368/150‑1 (MQ=255)
                                                                                         gCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCg                                                              <  1:6108573/150‑1 (MQ=255)
                                                                                          cTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCgg                                                             >  1:3406248/1‑150 (MQ=255)
                                                                                           tGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGa                                                            >  1:1630924/1‑150 (MQ=255)
                                                                                             tCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACg                                                          >  2:1542947/1‑150 (MQ=255)
                                                                                               tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                        >  2:5504283/1‑150 (MQ=255)
                                                                                               tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                        >  2:4495690/1‑150 (MQ=255)
                                                                                               tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                        >  1:2174940/1‑150 (MQ=255)
                                                                                                  tGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGc                                                     <  1:1125424/150‑1 (MQ=255)
                                                                                                    cAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCgg                                                   <  1:1414228/150‑1 (MQ=255)
                                                                                                     agagCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCggg                                                  >  1:1110948/1‑150 (MQ=255)
                                                                                                      gagCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGt                                                 <  1:3223905/150‑1 (MQ=255)
                                                                                                         cagcTGATTCGCTTGGGCGCTGGTTTGCAGTCTCTGTATGAGCCCTGTCTCAACACCACTACGGCTGCTGGTCAGATGATATTGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCga                                              <  1:1776691/150‑1 (MQ=255)
                                                                                                         cagcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCga                                              <  1:6600400/150‑1 (MQ=255)
                                                                                                          agcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTAAGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACAGGACGAGAGCGGGTCgag                                             >  1:1889157/1‑150 (MQ=255)
                                                                                                           gcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCgaga                                            <  2:3314265/150‑1 (MQ=255)
                                                                                                           gcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCgaga                                            <  1:5160041/150‑1 (MQ=255)
                                                                                                           gcTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCgaga                                            <  2:4993195/150‑1 (MQ=255)
                                                                                                            cTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGagag                                           >  1:3104967/1‑150 (MQ=255)
                                                                                                             tGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAgg                                          >  2:64412/1‑150 (MQ=255)
                                                                                                             tGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCCTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAgg                                          >  2:64315/1‑150 (MQ=255)
                                                                                                               caTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGAGAATGACGATATTGGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCa                                        <  2:4426581/149‑1 (MQ=255)
                                                                                                               aaTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCa                                        >  1:5272950/1‑150 (MQ=255)
                                                                                                               aaTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCa                                        >  1:701442/1‑150 (MQ=255)
                                                                                                                aTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAg                                       <  1:5886185/150‑1 (MQ=255)
                                                                                                                 tCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgc                                      <  2:3023672/150‑1 (MQ=255)
                                                                                                                   gCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgcgc                                    >  2:3032213/1‑150 (MQ=255)
                                                                                                                     ttGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCag                                  >  2:1972690/1‑150 (MQ=255)
                                                                                                                       gggCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGa                                                                                   <  2:4686082/99‑1 (MQ=255)
                                                                                                                       gggCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCagag                                >  2:2783536/1‑150 (MQ=255)
                                                                                                                       gggCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCagag                                >  2:3102907/1‑150 (MQ=255)
                                                                                                                           gcTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAg                            >  2:5404270/1‑150 (MQ=255)
                                                                                                                            cTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAgg                           <  2:6029868/150‑1 (MQ=255)
                                                                                                                                    cAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATaaaa                   >  1:2985339/1‑150 (MQ=255)
                                                                                                                                     aGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAAt                  >  2:557417/1‑150 (MQ=255)
                                                                                                                                        ctctGTCTGAGCCCTGGCTCAACACCACTACGGCTACTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCg               <  1:4952528/150‑1 (MQ=255)
                                                                                                                                         tctgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCgg              >  1:2525987/1‑150 (MQ=255)
                                                                                                                                           tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGACGGGCATAAAATTCGGTc            >  2:6094738/1‑150 (MQ=255)
                                                                                                                                           tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTc            >  1:923434/1‑150 (MQ=255)
                                                                                                                                           tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTc            >  2:121308/1‑150 (MQ=255)
                                                                                                                                            tgttgAGCCAGGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTcc           >  1:1458733/4‑150 (MQ=255)
                                                                                                                                             tctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCa          >  2:1714307/1‑150 (MQ=255)
                                                                                                                                             gttgAGCCAGGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCa          >  2:3471877/3‑150 (MQ=255)
                                                                                                                                                   cccTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCa    >  1:4202729/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCACG  >  NC_002947/2189064‑2189362

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