Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 715,019 Δ1 bp 6.1% coding (478/483 nt) PP_0607 ← type IV pili biogenesis protein FimT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947715,0190A. 6.1% 628.9 / 31.8 180coding (478/483 nt)PP_0607type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base A (88/81);  new base . (4/7);  total (92/88)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 185 total aligned reads displayed.

TGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAGC  >  NC_002947/714870‑715167
                                                                                                                                                      |                                                                                                                                                    
tGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTagca                                                                                                                                                      >  1:4490976/1‑150 (MQ=255)
   aCGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCg                                                                                                                                                   <  1:1048317/150‑1 (MQ=255)
     gggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                                 <  1:1281666/150‑1 (MQ=255)
     gggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGt                                                                                                                                                 <  1:301807/150‑1 (MQ=255)
       gcgcGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCa                                                                                                                                               >  2:2405333/1‑150 (MQ=255)
        cgcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCTCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCAt                                                                                                                                              <  1:2826380/150‑1 (MQ=255)
         gcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCAtt                                                                                                                                             <  1:4360670/150‑1 (MQ=255)
           gAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGCGCGGTCATTGc                                                                                                                                           <  1:1525036/150‑1 (MQ=255)
               agaACATCACCTTTTCAATGCCCAAAGAGCTGTGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAg                                                                                                                                       <  1:2115204/150‑1 (MQ=255)
               agaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAg                                                                                                                                       <  2:578813/150‑1 (MQ=255)
                gaACATCACCTTTTCAATGCCCACAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAgg                                                                                                                                      <  2:960143/150‑1 (MQ=255)
                    aTCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                                                  >  2:550681/1‑150 (MQ=255)
                    aTCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                                                  >  1:4519525/1‑150 (MQ=255)
                    aTCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcct                                                                                                                                  >  1:4441513/1‑150 (MQ=255)
                      cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGtcctcc                                                                                                                                <  2:818632/150‑1 (MQ=255)
                       aCCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                                                                               <  2:1234062/150‑1 (MQ=255)
                       aCCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                                                                               >  1:2008445/1‑150 (MQ=255)
                         cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTa                                                                                                                             >  1:4269432/1‑150 (MQ=255)
                         cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTAGCTGATCTGAAGGGCAGTAACGCCGCCAATGCAGGGGCTGCTTTTCATCCCTTTAGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGTTCATTGCTCAGGTCATCCGTa                                                                                                                             >  1:686164/1‑150 (MQ=255)
                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            <  1:3380845/150‑1 (MQ=255)
                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            >  1:3338560/1‑150 (MQ=255)
                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            <  2:3742361/150‑1 (MQ=255)
                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            <  2:1673425/150‑1 (MQ=255)
                          ttttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                                                                            <  2:4386153/150‑1 (MQ=255)
                           tttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAg                                                                                                                           <  1:247382/150‑1 (MQ=255)
                               aaTGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAAc                                                                                                                       >  2:2355220/1‑150 (MQ=255)
                                 tGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCt                                                                                                                     <  1:2265968/150‑1 (MQ=255)
                                   cccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                                   <  2:1815387/150‑1 (MQ=255)
                                   cccAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGa                                                                                                                   <  1:3307500/150‑1 (MQ=255)
                                    ccAAAGAGCTGCGGGTTCGCTCTCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTATCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                                                                                  >  1:2192649/1‑150 (MQ=255)
                                      aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGActc                                                                                                                >  1:3951015/1‑150 (MQ=255)
                                        agagCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTg                                                                                                              <  2:1536886/150‑1 (MQ=255)
                                        agagCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTg                                                                                                              <  1:2093469/150‑1 (MQ=255)
                                            cTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                                          <  1:2532502/150‑1 (MQ=255)
                                            cTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGc                                                                                                          <  1:3696891/150‑1 (MQ=255)
                                             tGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTc                                                                                                                                                <  2:3329444/111‑1 (MQ=255)
                                                 ggTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGAt                                                                                                     <  2:1777802/150‑1 (MQ=255)
                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                                                                                >  2:1108037/1‑150 (MQ=255)
                                                      gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                                                                                >  1:1664328/1‑150 (MQ=255)
                                                        tcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCGTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCCGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCCGca                                                                                              >  2:1867433/1‑150 (MQ=255)
                                                          gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                                                                            >  1:467632/1‑150 (MQ=255)
                                                          gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                                                                            >  1:1873491/1‑150 (MQ=255)
                                                          gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                                                                            >  1:2828379/1‑150 (MQ=255)
                                                          gctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                                                                            >  1:1312710/1‑150 (MQ=255)
                                                           ctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacca                                                                                           <  2:1393257/150‑1 (MQ=255)
                                                                cTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGa                                                                                      <  2:3444503/150‑1 (MQ=255)
                                                                cTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGa                                                                                      <  2:3482378/150‑1 (MQ=255)
                                                                cTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGa                                                                                      <  2:3952826/150‑1 (MQ=255)
                                                                 tGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGat                                                                                     >  1:1966116/1‑150 (MQ=255)
                                                                 tGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGat                                                                                     >  1:3862022/1‑150 (MQ=255)
                                                                    cGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATc                                                                                  <  2:1712685/150‑1 (MQ=255)
                                                                    cGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATc                                                                                  <  2:4144699/150‑1 (MQ=255)
                                                                     gggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                                 >  2:249700/1‑150 (MQ=255)
                                                                     gggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCg                                                                                 <  2:3040249/150‑1 (MQ=255)
                                                                      ggCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                                                                                >  1:500372/1‑150 (MQ=255)
                                                                       gCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGc                                                                               >  2:904645/1‑150 (MQ=255)
                                                                        cAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCt                                                                              <  1:1520097/150‑1 (MQ=255)
                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            <  2:751248/150‑1 (MQ=255)
                                                                          gTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGc                                                                            >  1:2484365/1‑150 (MQ=255)
                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGag                                                                         >  2:3821761/1‑150 (MQ=255)
                                                                             ccgccgCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGag                                                                         <  1:3356807/150‑1 (MQ=255)
                                                                              cgccgcCAATGCTGGGGCTGCTTTGCAGCCCTTGCGCGACACAATACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACACTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagc                                                                        >  2:2746504/1‑150 (MQ=255)
                                                                               gccgccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagca                                                                       >  1:1233642/1‑150 (MQ=255)
                                                                               gccgccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTACCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagca                                                                       >  1:1172006/1‑150 (MQ=255)
                                                                                  gccAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcg                                                                    >  2:2730794/1‑150 (MQ=255)
                                                                                      aTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc                                                                <  1:2687928/150‑1 (MQ=255)
                                                                                      aTGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc                                                                >  1:4041719/1‑150 (MQ=255)
                                                                                        gCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGt                                                              >  2:1626719/1‑150 (MQ=255)
                                                                                        gCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGt                                                              >  1:2435388/1‑150 (MQ=255)
                                                                                         cAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                             >  1:4128580/1‑150 (MQ=255)
                                                                                         cAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGtt                                                             >  1:231545/1‑150 (MQ=255)
                                                                                          aGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTc                                                            <  1:2177229/150‑1 (MQ=255)
                                                                                                tgctTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACAt                                                      >  1:3567101/1‑150 (MQ=255)
                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCa                                                 <  2:2262240/150‑1 (MQ=255)
                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCa                                                 <  2:2208035/150‑1 (MQ=255)
                                                                                                     tGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCa                                                 >  1:3712069/1‑150 (MQ=255)
                                                                                                         gCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGttt                                             <  1:3085081/150‑1 (MQ=255)
                                                                                                            cTTGCGCGACACAAGACCGATCCTGCAAAGACCTCGCTAGCAGCGGTCATTTCTCAGGTCCTCCGTAAGCAAACTGAATCTTCCGCTTGATGCCAGAACCACCTGATATAGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGATTTgc                                          >  2:3336982/1‑150 (MQ=255)
                                                                                                                    acacAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAAcc                                  >  2:718045/1‑150 (MQ=255)
                                                                                                                     cacaAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCa                                 <  2:2638880/150‑1 (MQ=255)
                                                                                                                      acaAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAg                                >  2:4046363/1‑150 (MQ=255)
                                                                                                                       caAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTTCGAGCAGCGGCGCGTTCGCTCATTTCCAGTTTTGCTCCCAACCAGc                               <  1:3941473/150‑1 (MQ=255)
                                                                                                                       caAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGc                               <  1:3941474/150‑1 (MQ=255)
                                                                                                                        aaGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCt                              >  2:5913/1‑150 (MQ=255)
                                                                                                                               cTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGccc                       <  1:565461/150‑1 (MQ=255)
                                                                                                                                tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg                      <  2:4490976/150‑1 (MQ=255)
                                                                                                                                tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg                      <  1:3138771/150‑1 (MQ=255)
                                                                                                                                tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg                      >  2:3592646/1‑150 (MQ=255)
                                                                                                                                tCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCg                      >  2:1133213/1‑150 (MQ=255)
                                                                                                                                 gcTGCATTGACCGCGCTAGCAGCGGTCTTTGCTCAGTTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCTTTCGCACATTTCCAGTTTTGCGCCCAACCTGCTGTTGCCCGa                     <  1:2514618/149‑1 (MQ=255)
                                                                                                                                  cTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGAc                    <  2:3656530/150‑1 (MQ=255)
                                                                                                                                       aTGACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCAt               >  2:1046943/1‑150 (MQ=255)
                                                                                                                                        tGACCGCTGCTAGCGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATg              >  1:1272897/1‑150 (MQ=255)
                                                                                                                                          aCCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATggg            >  1:345850/1‑150 (MQ=255)
                                                                                                                                            cgcgCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCa          >  1:1530094/1‑150 (MQ=255)
                                                                                                                                             gcgcTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATgctcac         >  2:3604812/1‑145 (MQ=255)
                                                                                                                                                gcTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGc       <  1:653287/150‑1 (MQ=255)
                                                                                                                                                gcTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCCGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGGTGCCCGACGTCATGGGCACGc       <  1:550681/150‑1 (MQ=255)
                                                                                                                                                  tAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGccc     >  2:211236/1‑150 (MQ=255)
                                                                                                                                                    gcagcGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAg   <  1:2582336/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCGCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACATTTCCAGTTTTGCGCCCAACCAGCTGTTGCCCGACGTCATGGGCACGCCCAGC  >  NC_002947/714870‑715167

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: