Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 965,473 T→A 5.0% L385Q (CTG→CAG)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913965,4730TA 5.0% 318.5 / 3.0 97L385Q (CTG→CAG) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/50);  new base A (2/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATATTTTGGTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTA  >  NC_000913/965328‑965621
                                                                                                                                                 |                                                                                                                                                    
ttGATATTTTGGTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCAGAc                                                                                                                                                   >  5:268861/1‑149 (MQ=255)
      ttttGGTCGCAGTGGGTTCCCTGTCCTGAGTGGCACCTGCCGCCGGTCTTGCGTGCAGTTGTTTCCCTCAGCCATCTGCACCTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGtt                                                                                                                                             <  2:287673/149‑1 (MQ=255)
      ttttGGTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGt                                                                                                                                              <  3:201658/148‑1 (MQ=255)
         tGGTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTa                                                                                                                                           <  6:32086/148‑1 (MQ=255)
           gTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTAtt                                                                                                                                         <  3:247257/148‑1 (MQ=255)
           gTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTAtt                                                                                                                                         <  1:40754/148‑1 (MQ=255)
           gTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTg                                                                                                                                        >  5:48215/1‑149 (MQ=255)
                aGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAAt                                                                                                                                   >  8:334536/1‑149 (MQ=255)
                aGTGGTTTCCCTGTCATGAGTGGCAACTGACGCCGGTATTGCGTGCAGTTGTTTCCCTCATCAATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCATATCGTCATATTTCATGGCATTAGTCGGACCTCGTTTATTGCAAAt                                                                                                                                   >  2:170165/1‑149 (MQ=255)
                 gTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTa                                                                                                                                           >  1:21921/1‑140 (MQ=255)
                 gTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTa                                                                                                                                           <  2:21921/140‑1 (MQ=255)
                 gTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATc                                                                                                                                  >  1:180897/1‑149 (MQ=255)
                  tGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCGCGTTTATTGCAAATCt                                                                                                                                 <  3:7966/149‑1 (MQ=255)
                    gTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTa                                                                                                                                <  1:283784/148‑1 (MQ=255)
                    gTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTAt                                                                                                                               <  2:180897/149‑1 (MQ=255)
                     tttCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCAGACCTCGTTTATTGCAAATCTAtt                                                                                                                              <  4:146824/149‑1 (MQ=255)
                        cccTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAg                                                                                                                            <  6:48215/148‑1 (MQ=255)
                            gTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTa                                                                                                                             >  3:82366/1‑143 (MQ=255)
                            gTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTa                                                                                                                             <  4:82366/143‑1 (MQ=255)
                             tCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCt                                                                                                                                                      >  2:186386/1‑117 (MQ=255)
                             tCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCt                                                                                                                                                      <  1:186386/117‑1 (MQ=255)
                             tCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCAGACCTCGTTTATTGCAAATCTATTAGCAAtt                                                                                                                       <  6:268861/148‑1 (MQ=255)
                               cTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTa                                                                                                                             >  6:370787/1‑140 (MQ=255)
                               cTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTa                                                                                                                             <  5:370787/140‑1 (MQ=255)
                                tGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTAt                                                                                                                               >  8:86581/1‑137 (MQ=255)
                                tGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTAt                                                                                                                               <  7:86581/137‑1 (MQ=255)
                                    tGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTg                                                                                                                 >  5:207926/1‑147 (MQ=255)
                                    tGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTg                                                                                                                 <  6:207926/147‑1 (MQ=255)
                                    tGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGt                                                                                                               <  8:195557/149‑1 (MQ=255)
                                      gCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTg                                                                                                              >  8:298633/1‑148 (MQ=255)
                                      gCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTg                                                                                                              <  7:298633/148‑1 (MQ=255)
                                          cTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACAtt                                                                                                         <  7:132346/149‑1 (MQ=255)
                                               gccgGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATcg                                                                                                    >  7:104354/1‑148 (MQ=255)
                                               gccgGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCa                                                                                                    >  1:3111/1‑149 (MQ=255)
                                                   gTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCGCGGt                                                                                                >  6:185982/1‑149 (MQ=255)
                                                            gCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGAt                                                                                       <  2:3111/149‑1 (MQ=255)
                                                              aGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATcc                                                                                     >  7:253173/1‑149 (MQ=255)
                                                                ttgttTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCAGACCTCGTTTATTGCAAATCTATTGGCAATTCCCTTGGt                                                                                                               <  6:285342/121‑1 (MQ=255)
                                                                ttgttTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCAGACCTCGTTTATTGCAAATCTATTGGCAATTCCCTTGGt                                                                                                               >  5:285342/1‑121 (MQ=255)
                                                                   tttCCCTCATCCATCTTCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATAGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTcgc                                                                                 <  7:179422/148‑1 (MQ=255)
                                                                     tCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACa                                                                                                           <  8:149277/120‑1 (MQ=255)
                                                                     tCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACa                                                                                                           >  7:149277/1‑120 (MQ=255)
                                                                     tCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTcgccgc                                                                              >  1:174071/1‑149 (MQ=255)
                                                                       ccTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTc                                                                                   <  6:148819/142‑1 (MQ=255)
                                                                       ccTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTc                                                                                   >  5:148819/1‑142 (MQ=255)
                                                                       ccTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGAAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGa                                                                            >  2:152200/1‑149 (MQ=255)
                                                                       ccTCATCCATCTGCAACTGGGAATCACAATTATGCTTATGCCAGTGCAAATCGTCATATTTCATGGAATTAGTCTGACCTCGTTTATTGCAATTCTATTATCAATTCCCTTGGTGACATTTATCACGGTTACGTTGATCCTCGCcccg                                                                             >  4:91305/1‑148 (MQ=255)
                                                                        cTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGAt                                                                           >  2:76367/1‑149 (MQ=255)
                                                                           atccatcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGt                                                                        >  5:274771/1‑149 (MQ=255)
                                                                           aaccatcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATATCATGGCATTACTCTGACCTCCTATATTGCAAATCTATTAGCAATTCCCTTTGTGACATTTATCACGGTTCCGTTGATCCTCCCCGCGATGGt                                                                        >  4:238846/3‑149 (MQ=255)
                                                                            tccatcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCGCGGTTCCGTTGATCCTCGCCGCGATGGtt                                                                       <  5:185982/149‑1 (MQ=255)
                                                                              catcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCGCGGTTCCGTTGATCCTCGCCGCGATGGTtgt                                                                     <  8:104354/149‑1 (MQ=255)
                                                                              catcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTtgt                                                                     <  8:253173/149‑1 (MQ=255)
                                                                                  tGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGCCCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGCTGACATTTATCACGGTTCCGTTGATCCTAGCCGGGAGGGTTGTGCAt                                                                 >  3:221265/1‑149 (MQ=255)
                                                                                  tGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCtt                                                                                                                  <  8:324489/100‑1 (MQ=255)
                                                                                  tGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCAt                                                                 >  5:342292/1‑149 (MQ=255)
                                                                                  tGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCAt                                                                 >  8:11590/1‑149 (MQ=255)
                                                                                   gCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATgg                                                                         <  8:207885/140‑1 (MQ=255)
                                                                                   gCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATgg                                                                         >  7:207885/1‑140 (MQ=255)
                                                                                      aCTGGTAATCAAACTTCTGGTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCGTTGGTTACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTaa                                                             <  4:29095/149‑1 (MQ=255)
                                                                                        tGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGtt                                                                       >  6:297463/1‑137 (MQ=255)
                                                                                        tGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGtt                                                                       <  5:297463/137‑1 (MQ=255)
                                                                                        tGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGc                                                           >  2:202705/1‑149 (MQ=255)
                                                                                        tGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGc                                                           >  2:153558/1‑149 (MQ=255)
                                                                                          ggAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGc                                                           <  1:355095/147‑1 (MQ=255)
                                                                                          ggAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCgg                                                         <  2:174071/149‑1 (MQ=255)
                                                                                            aaTCACACTTCTGCTTTTGCACGTGCAAATCGTCATATTTCATGGCAATAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATAACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCTGGc                                                       <  3:238846/149‑1 (MQ=255)
                                                                                             aTCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGcc                                                      >  8:38224/1‑149 (MQ=255)
                                                                                                    ttCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGACGCGATGGTTGTGCATTTAAGCGGGCCGTTAATc                                               >  3:52932/1‑149 (MQ=255)
                                                                                                      cTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATc                                               <  7:334536/147‑1 (MQ=255)
                                                                                                         cTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                          <  1:202705/149‑1 (MQ=255)
                                                                                                             tGCCCGTGCAAATCGTCTTATTTCATGGCTTTTTTCTGACCTCGTTTATTGCAAATCTTTTAGCATTTCGCTTTGTTACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCaa                                      <  3:91305/149‑1 (MQ=255)
                                                                                                               cccGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAgg                                    <  7:11590/149‑1 (MQ=255)
                                                                                                                ccGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAggg                                   >  7:194682/1‑149 (MQ=255)
                                                                                                                ccGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAggg                                   <  5:144177/149‑1 (MQ=255)
                                                                                                                 cGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGt                                  <  1:153558/149‑1 (MQ=255)
                                                                                                                   tGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTa                                <  4:52932/149‑1 (MQ=255)
                                                                                                                    gaaaaTCGGCATATTTAATGGCATGAGTCTGGCCTCGTTTATTGCAAATCTATTAGCAATTACCTTTGTTAAAGTTATCTCTGGTCCTTTGATACTAGAAGCGATTTTTGTGCATTTAAGCGGGCCGTTAAGCCTGGATCAAGGGTGAt                               >  5:95550/3‑149 (MQ=255)
                                                                                                                    gCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCAtt                                                                >  2:209081/1‑116 (MQ=255)
                                                                                                                    gCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCAtt                                                                <  1:209081/116‑1 (MQ=255)
                                                                                                                     cAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGTAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATg                              <  1:152200/149‑1 (MQ=255)
                                                                                                                     cAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTAt                               <  1:76367/148‑1 (MQ=255)
                                                                                                                      aaaTCGTCATATTTCATGGCATTAGACTCACCTCGTTTATTGCAAATCTATTATCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATgg                             >  1:254669/1‑149 (MQ=255)
                                                                                                                       cttCGTCATATCTCATGGAATTCGTCTGTCCTCGTGTCTTTCAAATCGATTCGCAATTGCGTTGGTGACTTTTCTCACGGTTCCGTTGATCCTCGCCGCGATGGATGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTGAGGGt                            <  4:42768/147‑1 (MQ=255)
                                                                                                                        aTCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                          >  2:101802/1‑134 (MQ=255)
                                                                                                                        aTCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                          <  1:101802/134‑1 (MQ=255)
                                                                                                                            tCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCtt                       <  1:122891/149‑1 (MQ=255)
                                                                                                                            tCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGAGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCtt                       >  8:379641/1‑149 (MQ=255)
                                                                                                                            tCATATTTCATGGCAGTAGTCTGCCCTCGTTTATTGCAAATCTATTAGCAAGTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAGTCCTGGAGCAAGGGTTATGGTTTCtt                       >  8:20118/1‑149 (MQ=255)
                                                                                                                                atTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGtt                                                                       <  7:269800/97‑1 (MQ=255)
                                                                                                                                atTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGtt                                                                       >  8:269800/1‑97 (MQ=255)
                                                                                                                                atTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTtgt                                                                     >  6:7776/1‑99 (MQ=255)
                                                                                                                                atTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTtgt                                                                     <  5:7776/99‑1 (MQ=255)
                                                                                                                                 tttCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGc                                                                   <  3:151862/100‑1 (MQ=255)
                                                                                                                                 tttCATGGCATTAGTCTGACATCGTTTATTGAAAATCTATTACCAATTACCTTGGTGACATTTATCACGGTTACGTTGATCCTCGCCGCCATGGTAGTGc                                                                   >  4:151862/1‑100 (MQ=255)
                                                                                                                                   tCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTATTCCTGGAGCAAGGGTTATTGTTTCTTGCCGAc                 <  2:254669/148‑1 (MQ=255)
                                                                                                                                          aTTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCttt         <  1:170165/149‑1 (MQ=255)
                                                                                                                                             aGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTgg       <  7:38224/148‑1 (MQ=255)
                                                                                                                                              gTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCt     >  6:285292/1‑149 (MQ=255)
                                                                                                                                                 tGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTa  >  5:132649/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTGATATTTTGGTATCAGTGGTTTCCCTGTCCTGAGTGGCAACTGCCGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTA  >  NC_000913/965328‑965621

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