Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 923,487 A→G E75G (GAG→GGG)  clpA → ATP‑dependent Clp protease ATP‑binding subunit ClpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913923,4870AG98.5% 231.3 / ‑6.1 66E75G (GAG→GGG) clpAATP‑dependent Clp protease ATP‑binding subunit ClpA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/35);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGCGCTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGAGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAGGAGTCGCAGGC  >  NC_000913/923356‑923622
                                                                                                                                   |                                                                                                                                       
aCTGGCGCTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGaca                                                                                                                                 <  1:134040/140‑1 (MQ=255)
    gcgcTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGAcacaca                                                                                                                             >  2:299658/1‑140 (MQ=255)
     cgcTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGTGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAg                                                                                                                            >  2:504983/1‑140 (MQ=255)
     cgcTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAg                                                                                                                            >  1:259774/1‑140 (MQ=255)
     cgcTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAg                                                                                                                            >  2:513522/1‑140 (MQ=255)
       ctgctCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGcc                                                                                                                          >  2:837580/1‑140 (MQ=255)
         gctCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGa                                                                                                                        >  1:20972/1‑140 (MQ=255)
         gctCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGa                                                                                                                        >  1:506262/1‑140 (MQ=255)
          ctCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGAc                                                                                                                       <  1:299658/140‑1 (MQ=255)
                    cATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGtttt                                                                                                             >  2:262667/1‑140 (MQ=255)
                    cATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGtttt                                                                                                             >  2:473324/1‑140 (MQ=255)
                      tCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCa                                                                                                           <  1:473324/140‑1 (MQ=255)
                       cTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAg                                                                                                          >  2:601197/1‑140 (MQ=255)
                       cTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAg                                                                                                          >  1:75318/1‑140 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:218307/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:774346/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:513522/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  2:740484/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:748828/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:601197/140‑1 (MQ=255)
                        tGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  2:506262/140‑1 (MQ=255)
                         gCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         >  2:800599/1‑139 (MQ=255)
                         gCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGc                                                                                                         <  1:800599/139‑1 (MQ=255)
                         gCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGAGCGCGACACACAGCCGACGCTGAGTTTTCAGCg                                                                                                        >  2:946593/1‑140 (MQ=255)
                             gggAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTa                                                                                                    >  2:698423/1‑140 (MQ=255)
                             gggAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTa                                                                                                    <  2:624342/140‑1 (MQ=255)
                               gAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCa                                                                                                           >  1:720805/1‑131 (MQ=255)
                               gAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCa                                                                                                           <  2:720805/131‑1 (MQ=255)
                                 ggCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGc                                                                                                <  2:467186/140‑1 (MQ=255)
                                    gcTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAAc                                                                                             >  1:263575/1‑140 (MQ=255)
                                      tGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGt                                                                                           <  1:583377/140‑1 (MQ=255)
                                       ggAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTg                                                                                          <  1:278490/140‑1 (MQ=255)
                                        gAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGc                                                                                         >  2:490509/1‑140 (MQ=255)
                                          aGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCgg                                                                                       <  1:411633/140‑1 (MQ=255)
                                             gtgtTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCt                                                                                    >  1:882595/1‑140 (MQ=255)
                                             gtgtTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCt                                                                                    >  1:873435/1‑140 (MQ=255)
                                                   tgtgGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATg                                                                              >  2:847936/1‑140 (MQ=255)
                                                             gTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGtcctc                                                                    <  1:374065/140‑1 (MQ=255)
                                                                 cgcTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATg                                                                              <  2:191085/126‑1 (MQ=255)
                                                                 cgcTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATg                                                                              >  1:191085/1‑126 (MQ=255)
                                                                 cgcTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGt                                                                <  1:262667/140‑1 (MQ=255)
                                                                     ccGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCa                                                            <  2:873435/140‑1 (MQ=255)
                                                                     ccGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCa                                                            <  2:263575/140‑1 (MQ=255)
                                                                       gTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAAt                                                          >  1:734715/1‑140 (MQ=255)
                                                                                tGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTAcc                                                 >  2:120200/1‑140 (MQ=255)
                                                                                tGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTAcc                                                 <  2:259774/140‑1 (MQ=255)
                                                                                  gAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCgg                                               <  1:847936/140‑1 (MQ=255)
                                                                                              gAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCt                                   <  1:490509/140‑1 (MQ=255)
                                                                                              gAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCt                                   <  1:120200/140‑1 (MQ=255)
                                                                                                aCAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTgg                                 <  2:882595/140‑1 (MQ=255)
                                                                                                    aCCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGc                             >  1:749924/1‑140 (MQ=255)
                                                                                                    aCCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGc                             >  1:902194/1‑140 (MQ=255)
                                                                                                    aCCACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGc                             >  2:794196/1‑140 (MQ=255)
                                                                                                     ccACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCt                            <  2:902194/140‑1 (MQ=255)
                                                                                                     ccACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCt                            <  1:395918/140‑1 (MQ=255)
                                                                                                     ccACACCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCt                            <  2:734715/140‑1 (MQ=255)
                                                                                                      cacaCCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTa                           >  1:132881/1‑140 (MQ=255)
                                                                                                         acCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCt                        >  2:280539/1‑140 (MQ=255)
                                                                                                         acCCGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCt                        >  2:804587/1‑140 (MQ=255)
                                                                                                            cGTTCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTa                     <  2:749924/140‑1 (MQ=255)
                                                                                                              ttCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCt                                   >  1:167617/1‑124 (MQ=255)
                                                                                                              ttCTGCCTGCCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCt                                   <  2:167617/124‑1 (MQ=255)
                                                                                                                ctgcctgcCAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGa                 >  2:845890/1‑140 (MQ=255)
                                                                                                                   cctgccAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACa              <  1:845890/140‑1 (MQ=255)
                                                                                                                     tgccAGTGAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAgg            >  1:591934/1‑140 (MQ=255)
                                                                                                                            gAAGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAGGAGTCGCa     <  1:280539/140‑1 (MQ=255)
                                                                                                                             aaGAGGGGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAGGAGTCGCAg    <  1:794196/140‑1 (MQ=255)
                                                                                                                               gaggggCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAGGAGTCGCAggc  <  2:20972/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
ACTGGCGCTGCTCAGTAACCCATCTGCCCGGGAGGCGCTGGAAGCGTGTTCTGTGGATTTGGTTGCGCTCCGTCAGGAACTGGAAGCCTTTATTGAACAAACCACACCCGTTCTGCCTGCCAGTGAAGAGGAGCGCGACACACAGCCGACGCTGAGTTTTCAGCGTGTACTGCAACGTGCGGTCTTCCATGTCCAGTCCTCCGGTCGCAATGAGGTTACCGGTGCAAACGTTCTGGTCGCTATCTTTAGCGAACAGGAGTCGCAGGC  >  NC_000913/923356‑923622

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