Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,935,301 C→A 5.5% intergenic (‑54/+1092) USA300HOU_1798 ← / ← hsdS2 hypothetical protein/type I site‑specific deoxyribonuclease specificity subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,935,3010CA 5.5% 215.0 / 2.6 75intergenic (‑54/+1092)USA300HOU_1798/hsdS2hypothetical protein/type I site‑specific deoxyribonuclease specificity subunit
Reads supporting (aligned to +/- strand):  ref base C (29/39);  new base A (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

GGATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAATA  >  CP000730/1935168‑1935441
                                                                                                                                      |                                                                                                                                             
ggATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCgg                                                                                                                                        <  2:851688/141‑1 (MQ=255)
ggATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAAGTCTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCgg                                                                                                                                        <  2:836960/141‑1 (MQ=255)
   tttGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTaa                                                                                                                                     >  2:317315/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAATCTTTAATCGATAATCTTCTTCCTTTTTATTTGTTAAATACTC‑TTTTTTATAGTGTCCTCTATTTTTTTCTTTTTTCACTTTGCGGTAACTTGaata                                                                                                                             >  1:683364/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGaata                                                                                                                             >  1:147174/1‑141 (MQ=255)
           aTTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGaata                                                                                                                             >  1:664715/1‑141 (MQ=255)
               ttATTTTTAATATTATTATTTTGATTACCGGTAAATATAAACTTTAATCGATAACCTTCTTCCAGTTTATGTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTTCTTTTTACACGTTGCGGTAACTTGaatactaa                                                                                                                         >  1:158667/1‑141 (MQ=255)
                     ttAATATTATTATTTTGATTTCCTATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGaataatgatgtcga                                                                                                                   <  2:367855/141‑6 (MQ=255)
                     ttAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTAt                                                                                                                   <  2:855297/141‑1 (MQ=255)
                     ttAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTAt                                                                                                                   <  1:569636/141‑1 (MQ=255)
                         tattattatTTTGTTTTCCAATAAATATAAACTTTATTCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGt                                                                                                               <  1:597832/141‑1 (MQ=255)
                            tattattTTGATTTCCAATCAATATAAACTTTACTCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatt                                                                                                            <  2:443928/141‑1 (MQ=255)
                            tattattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGTTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtctt                                                                                                            >  2:790031/1‑141 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           >  2:203784/1‑141 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           <  1:844171/141‑1 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           <  2:426559/141‑1 (MQ=255)
                             attattTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGtatta                                                                                                           >  1:970350/1‑141 (MQ=255)
                                  tttGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCtat                                                                                                      <  1:953838/141‑1 (MQ=255)
                                                  tataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGt                                                                                     >  1:845279/1‑141 (MQ=255)
                                                  tataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGt                                                                                     <  1:1044804/141‑1 (MQ=255)
                                                   ataAACTTTACTCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTTCACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTg                                                                                    <  2:1010581/141‑1 (MQ=255)
                                                   ataAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTg                                                                                    <  2:1058316/141‑1 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTCTTTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCa                                                                                  >  1:597518/1‑141 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAt                                                                                 >  1:104641/1‑141 (MQ=255)
                                                      aaCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAt                                                                                 <  1:778650/141‑1 (MQ=255)
                                                       aCTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAtt                                                                                <  2:891400/141‑1 (MQ=255)
                                                        cTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAttt                                                                               <  2:941157/141‑1 (MQ=255)
                                                          ttAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAttttt                                                                             <  1:500001/141‑1 (MQ=255)
                                                           tAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCAtttttt                                                                            <  2:272946/141‑1 (MQ=255)
                                                             aTCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTATTTTTTTGTTGTTTTAATAATGATTTTATTCTTATTTTTTATTATATGTAAT‑GCTGGTGTCtttttttt                                                                          >  2:912886/1‑140 (MQ=255)
                                                              tCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTTTTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGt                                                                         >  2:11372/1‑141 (MQ=255)
                                                                gATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATtttttttttt                                                                       >  2:590348/1‑141 (MQ=255)
                                                                gATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGttt                                                                       <  2:406526/141‑1 (MQ=255)
                                                                  tAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGttttt                                                                     >  2:1059899/1‑141 (MQ=255)
                                                                      cttcttCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTttatt                                                                 <  1:985136/141‑1 (MQ=255)
                                                                         cttcTATTTTATTTGTTAAATACTC‑TTTTTTTTTTTGTTTTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTttttgttttttgttttt                                                              >  2:951927/1‑138 (MQ=255)
                                                                          ttcCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGa                                                             <  2:1079245/141‑1 (MQ=255)
                                                                            ccATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGAtt                                                           <  1:479823/141‑1 (MQ=255)
                                                                            ccATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGAtt                                                           <  1:251018/141‑1 (MQ=255)
                                                                              aTTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATtct                                                         >  1:122138/1‑141 (MQ=255)
                                                                                tttATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCa                                                       <  2:970350/141‑1 (MQ=255)
                                                                                   aTTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGATAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTTCATTTTTTTGTTTTTTATTATGATTCTCaatt                                                    >  2:486158/1‑139 (MQ=255)
                                                                                      tGTTAAATACTC‑TTTTTTTTTTTTTCTTCTAATTTTTGCGGGTTATACGTTGTGGTAACTTGAGTAATGATTTTATTGGGATTATCTATAAGATCAAAT‑TCTGTTGCATTTTTTTGTTTTTtattatttttttcttgtgat                                                 >  2:52957/1‑127 (MQ=255)
                                                                                      tGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTTTTTTTTGTAACTTGAATAATAATTTTATTCGTATTATCTATATGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATTATtcttattttat                                                 >  2:404856/1‑133 (MQ=255)
                                                                                      tGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTTTTTTTTTTTATGATTCTCAAGTGaa                                                 >  2:902013/1‑141 (MQ=255)
                                                                                       gTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTTTTTTTTTTTCTTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCaatttatt                                                >  2:400070/1‑135 (MQ=255)
                                                                                            aTACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCa                                           >  2:257787/1‑141 (MQ=255)
                                                                                              aCTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCata                                         <  2:95962/141‑1 (MQ=255)
                                                                                               ctc‑tTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCatat                                        <  2:868838/141‑1 (MQ=255)
                                                                                                 c‑tTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATAc                                      <  2:276394/141‑1 (MQ=255)
                                                                                                 c‑tTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATAc                                      <  1:163884/141‑1 (MQ=255)
                                                                                                    tttttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACaa                                    <  2:997564/141‑1 (MQ=255)
                                                                                                    tttttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACaa                                    <  1:138688/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 >  1:780349/1‑141 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:504282/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:147174/141‑1 (MQ=255)
                                                                                                       ttAGAGTGTCCTCTATTTTTTGCTTTTTACAAGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAatat                                 <  2:664715/141‑1 (MQ=255)
                                                                                                        tAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTTTTATTCTTAAGTGAATACCCATATACattatt                                >  2:379712/1‑141 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGTGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  2:646367/141‑1 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  1:19785/141‑1 (MQ=255)
                                                                                                         agagTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATAttt                               <  1:826141/141‑1 (MQ=255)
                                                                                                           agTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTaa                             <  2:2760/141‑1 (MQ=255)
                                                                                                           agTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTaa                             >  1:719606/1‑141 (MQ=255)
                                                                                                            gtgtCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAAc                            >  1:530269/1‑141 (MQ=255)
                                                                                                              gtCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAg                          <  2:1028286/141‑1 (MQ=255)
                                                                                                                      ttttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATTATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaa                  >  2:427676/1‑141 (MQ=255)
                                                                                                                      ttttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaa                  >  1:820663/1‑141 (MQ=255)
                                                                                                                       tttttGCTTTTTACTCGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATTAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATTATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaa                 <  2:154290/141‑1 (MQ=255)
                                                                                                                        ttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaa                >  1:255304/1‑141 (MQ=255)
                                                                                                                        ttttGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaa                >  1:262495/1‑141 (MQ=255)
                                                                                                                           tGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTaaaaaaa             <  1:971137/141‑1 (MQ=255)
                                                                                                                            gCTTTTTACACGTTGAGGTAGCTTGAATAATAATTTTATTCGTATTATCTAAAAGATCATAT‑GCTAGTGCATTTTATTATTTTTTATTAGGATTCGCAAGTAAAAACCCATCAAAAATAATTAACAGGTATCTAAAAAAat            <  2:923656/141‑1 (MQ=255)
                                                                                                                             cTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAata           >  1:453790/1‑141 (MQ=255)
                                                                                                                              tttttACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAatat          >  1:931675/1‑141 (MQ=255)
                                                                                                                                 ttACTCGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTATAAAAAAATATTCg       >  2:636757/1‑141 (MQ=255)
                                                                                                                                    cacGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGaa    <  1:714218/141‑1 (MQ=255)
                                                                                                                                     acGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAat   >  1:506507/1‑141 (MQ=255)
                                                                                                                                      cGTTGTGGTAACTTGAATAATAATTTTATTCGTATTATCTACATATTCATTTCTTTGGTG‑TTTTTGGTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAATAGGTCTCTAAAAAAATATTCGGAata  <  1:431954/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                             
GGATTTGAAAAATTCTTATTTTTAATATTATTATTTTGATTTCCAATAAATATAAACTTTAATCGATAACCTTCTTCCATTTTATTTGTTAAATACTC‑TTTTTTAGAGTGTCCTCTATTTTTTGCTTTTTACACGTTGCGGTAACTTGAATAATAATTTTATTCGTATTATCTATAAGATCAAAT‑GCTGGTGCATTTTTTTGTTTTTTATTATGATTCTCAAGTGAATACCCATATACAATATTTAACAGGTCTCTAAAAAAATATTCGGAATA  >  CP000730/1935168‑1935441

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