Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,189,214 T→G 5.5% W94G (TGG→GGG)  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,2140TG 5.5% 922.4 / 28.5 269W94G (TGG→GGG) PP_5490site‑specific DNA recombinase
Reads supporting (aligned to +/- strand):  ref base T (128/125);  new base G (3/12);  total (132/137)
Fisher's exact test for biased strand distribution p-value = 3.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

AGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCACGG  >  NC_002947/2189065‑2189363
                                                                                                                                                     |                                                                                                                                                     
agTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCt                                                                                                                                                       <  2:5958525/150‑1 (MQ=255)
agTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCt                                                                                                                                                       <  1:4872200/150‑1 (MQ=255)
 gTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTg                                                                                                                                                      <  2:6202203/150‑1 (MQ=255)
 gTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTg                                                                                                                                                      <  2:1202832/150‑1 (MQ=255)
      tgaGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCa                                                                                                                                                 <  2:530062/150‑1 (MQ=255)
       gaGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCag                                                                                                                                                <  1:5893751/150‑2 (MQ=255)
        aGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCAGGGCTCaga                                                                                                                                               <  2:4326274/150‑3 (MQ=255)
          gctgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAacac                                                                                                                                             >  1:1509539/1‑150 (MQ=255)
            tgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAAcacca                                                                                                                                           >  2:6844702/1‑150 (MQ=255)
            tgctTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAAcacca                                                                                                                                           <  1:3552215/150‑1 (MQ=255)
               ttGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTa                                                                                                                                        <  1:5769686/150‑1 (MQ=255)
                      tctcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgc                                                                                                                                 <  2:569897/150‑1 (MQ=255)
                      tctcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGgctgc                                                                                                                                 <  1:5483240/150‑1 (MQ=255)
                        tcAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGctgctg                                                                                                                               <  1:725036/150‑1 (MQ=255)
                         caggccaggcGATGTGGTTGTGGTGAGTCGACTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTgg                                                                                                                              >  2:3386208/1‑150 (MQ=255)
                          aggccaggcGATGTGGGTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGt                                                                                                                             >  1:6339319/1‑150 (MQ=255)
                           ggccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTc                                                                                                                            <  1:4564879/150‑1 (MQ=255)
                            gccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTcg                                                                                                                           >  2:4005822/1‑149 (MQ=255)
                            gccaggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCa                                                                                                                           >  2:4912815/1‑150 (MQ=255)
                               aggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgat                                                                                                                        <  1:3727175/150‑1 (MQ=255)
                                ggcGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAgatg                                                                                                                       >  2:2470612/1‑150 (MQ=255)
                                    aTGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATaa                                                                                                                   <  2:2997339/150‑1 (MQ=255)
                                     tgtggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAAt                                                                                                                  <  1:5784087/150‑1 (MQ=255)
                                       tggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGa                                                                                                                <  1:2774523/150‑1 (MQ=255)
                                       tggttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGGTGCCGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATgc                                                                                                                <  1:383044/150‑2 (MQ=255)
                                         gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACg                                                                                                              >  1:2586832/1‑150 (MQ=255)
                                         gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACg                                                                                                              >  2:5601290/1‑150 (MQ=255)
                                         gttgtggtGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACg                                                                                                              >  2:4530314/1‑150 (MQ=255)
                                             tggtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGatat                                                                                                          >  1:195504/1‑150 (MQ=255)
                                               gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                        <  2:1264704/150‑1 (MQ=255)
                                               gtgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATAttt                                                                                                        >  2:2555721/1‑150 (MQ=255)
                                                tgAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTg                                                                                                       <  2:6545256/150‑1 (MQ=255)
                                                 gAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGc                                                                                                      <  2:3131811/150‑1 (MQ=255)
                                                 gAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGc                                                                                                      <  2:2577263/150‑1 (MQ=255)
                                                  aGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCa                                                                                                     >  2:5233622/1‑150 (MQ=255)
                                                  aGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCa                                                                                                     >  2:727409/1‑150 (MQ=255)
                                                   gtcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAg                                                                                                    <  2:2169630/150‑1 (MQ=255)
                                                   gtcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAg                                                                                                    <  1:2107297/150‑1 (MQ=255)
                                                    tcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAgg                                                                                                   >  1:1006113/1‑150 (MQ=255)
                                                    tcgtcTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAgg                                                                                                   >  1:3375859/1‑150 (MQ=255)
                                                        cTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATc                                                                                               <  1:80004/150‑1 (MQ=255)
                                                         ttGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCg                                                                                              >  1:4857518/1‑150 (MQ=255)
                                                           tggcGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCt                                                                                            <  2:1611894/147‑1 (MQ=255)
                                                            aTCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCCCTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTa                                                                                           <  1:2572937/150‑1 (MQ=255)
                                                            aTCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTa                                                                                           >  1:4304493/1‑150 (MQ=255)
                                                             tCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAg                                                                                          >  2:4649305/1‑150 (MQ=255)
                                                             tCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAg                                                                                          >  2:4129308/1‑150 (MQ=255)
                                                                     cACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAgc                                                                                  <  2:5158878/150‑1 (MQ=255)
                                                                        ggTCGACAATTGCTTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCa                                                                               <  1:5439861/150‑1 (MQ=255)
                                                                        ggTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCa                                                                               <  2:306489/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:3745145/150‑1 (MQ=255)
                                                                               aaTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGAt                                                                        >  1:2716733/1‑150 (MQ=255)
                                                                                aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                       <  1:2570375/150‑1 (MQ=255)
                                                                                aTTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATc                                                                       <  1:6701010/150‑1 (MQ=255)
                                                                                 ttGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCg                                                                      >  1:5569083/1‑150 (MQ=255)
                                                                                       tGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGAcc                                                                <  1:6391995/150‑1 (MQ=255)
                                                                                        gCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCg                                                               <  1:6108573/150‑1 (MQ=255)
                                                                                          tGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGa                                                             >  1:1630924/1‑150 (MQ=255)
                                                                                              tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                         >  1:2174940/1‑150 (MQ=255)
                                                                                              tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                         >  2:5504283/1‑150 (MQ=255)
                                                                                              tATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACgag                                                         >  2:4454708/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                                                  >  2:1529914/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)
                                                                                                           cTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGagag                                            >  1:3104967/1‑150 (MQ=255)
                                                                                                            tGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAgg                                           >  2:64412/1‑150 (MQ=255)
                                                                                                            tGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAgg                                           >  2:3546964/1‑150 (MQ=255)
                                                                                                            tGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCCTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAgg                                           >  2:64315/1‑150 (MQ=255)
                                                                                                              aaTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCa                                         >  1:5272950/1‑150 (MQ=255)
                                                                                                              aaTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCa                                         >  1:701442/1‑150 (MQ=255)
                                                                                                               aTCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAg                                        <  1:3000392/150‑1 (MQ=255)
                                                                                                                  gCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAgcgc                                     >  2:3032069/1‑150 (MQ=255)
                                                                                                                   cTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCa                                    >  2:479462/1‑150 (MQ=255)
                                                                                                                    ttGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCag                                   >  2:1972690/1‑150 (MQ=255)
                                                                                                                    ttGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCag                                   >  1:5930987/1‑150 (MQ=255)
                                                                                                                      gggCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGa                                                                                    >  1:4686082/1‑99 (MQ=255)
                                                                                                                        gcgcTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGa                               >  2:2122093/1‑150 (MQ=255)
                                                                                                                          gcTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAg                             >  2:5404270/1‑150 (MQ=255)
                                                                                                                           cTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAgg                            <  2:6029868/150‑1 (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:6833490/150‑1 (MQ=255)
                                                                                                                                          tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTc             >  1:923434/1‑150 (MQ=255)
                                                                                                                                          tgtctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTc             >  2:2725457/1‑150 (MQ=255)
                                                                                                                                           tgttgAGCCAGGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTcc            >  1:1458733/4‑150 (MQ=255)
                                                                                                                                            tctgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCa           >  2:1714307/1‑150 (MQ=255)
                                                                                                                                              tgAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCaaa         <  2:3276224/150‑1 (MQ=255)
                                                                                                                                               gAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAg        >  2:1090472/1‑150 (MQ=255)
                                                                                                                                                aGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGc       >  1:5183419/1‑150 (MQ=255)
                                                                                                                                                  cccTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCa     >  1:4202729/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTTGATGAGGCTGCTTGATCATCTCAGGCCAGGCGATGTGGTTGTGGTGAGTCGTCTTGATCGGCTCGCACGGTCGACAATTGATTTGCTGTCTATTGCAGAGCAGCTGAATCGCTTGGGCGCTGGGTTGCAGTCTCTGTCTGAGCCCTGGCTCAACACCACTACGGCTGCTGGTCAGATGATAATGACGATATTTGCAGGCATCGCTAGCTTCGAGCGCAGCCTGATCGTAGACCGGACGAGAGCGGGTCGAGAGGCAGCGCAGAGGAAGGGCATAAAATTCGGTCCAAAGCCACGG  >  NC_002947/2189065‑2189363

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: