Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,767,668 C→A 5.4% intergenic (+315/‑40) rnlB → / → yfjP CP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,6680CA 5.4% 218.5 / 3.3 74intergenic (+315/‑40)rnlB/yfjPCP4‑57 prophage; uncharacterized protein/CP4‑57 prophage; 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGA  >  NC_000913/2767520‑2767814
                                                                                                                                                    |                                                                                                                                                  
tcatGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTcgc                                                                                                                                                    <  1:219336/147‑1 (MQ=255)
    gTCACCTCTCTTCTGTAAATCCTATCATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGtt                                                                                                                                                <  2:1219421/149‑1 (MQ=255)
      cACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCt                                                                                                                                              <  2:31047/149‑1 (MQ=255)
       acttaTTTGCGGGGAATCCTATGAGGCTGTTCCCTGTAAAGCGTCCTGGGAGTTCAGGCAGTTTCTGCAATCTATGTGCATTCACAAAAGCTTCAGATAGATATCCTTCTCCTGAACAATCCGATCACCGGCAGAACCTCGCGGTTCTg                                                                                                                                             <  2:522467/145‑1 (MQ=255)
              ttCTGGTAATCCTATGGTGCCTTTCCCTGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAACATCTTCCGATACCTATCTTTCTCCTGACCACTTCGATTAACAGCAGAACCTCTAGGTTCTGCCGttca                                                                                                                                      <  2:123613/149‑3 (MQ=255)
              tgatGGAGAGCCTGTGATGCCTTGCCCGGTAAATTGTGCTTTGAGGTCAGAAAGATTCTGCAATCGATGTGCAATGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCGCGAGGTTCTGCCGttca                                                                                                                                      <  2:852488/146‑3 (MQ=255)
                  ggAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtt                                                                                                                                  >  1:128253/1‑149 (MQ=255)
                    aaaTCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGtttt                                                                                                                                >  1:1035229/1‑149 (MQ=255)
                     actcttATGACTCCTTTCCCGGTAACTGGTGTCGTGAGTTCAGAAAGATTCTCCATTCTATGTGCATTGAAAAAATCTCCCGATAGATACCATTCTCCTGAACTATCCGATTAATGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTa                                                                                                                               <  2:1346535/144‑1 (MQ=255)
                      aTCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTAc                                                                                                                              <  1:93351/149‑1 (MQ=255)
                       tCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAACATCTTCAGATAGATATCATTCTCCTGAACAAGCCGATGAACGGCAGAACCTCGCGGTCATGCCGTTTTTGTTTTGCt                                                                                                                             >  2:1127380/1‑149 (MQ=255)
                         cTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTCTGTGCATTGAAAAAATCTCCAGCTATTTGTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCCGTTGTGCCCTTTTTGTTTTACttt                                                                                                                           >  1:320842/1‑149 (MQ=255)
                         cTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACttt                                                                                                                           >  1:79609/1‑149 (MQ=255)
                           atgatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTc                                                                                                                         <  1:801303/149‑1 (MQ=255)
                             gatgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGACCGGCAGAACCTCGCGGTTCTGCCGTTGTTGTTTTACTTTTAag                                                                                                                       >  1:445378/1‑149 (MQ=255)
                              atgCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCaga                                                                                                                      >  2:1037366/1‑149 (MQ=255)
                                 cctcacccGGTAAATGGTGCTGTGAGGTCAGAAAGATTCTGCAATCTAGGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGTACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagga                                                                                                                   <  1:391541/144‑1 (MQ=255)
                                 ccTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgagga                                                                                                                   <  2:97361/149‑1 (MQ=255)
                                  cTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAgaggag                                                                                                                  >  1:1352216/1‑149 (MQ=255)
                                     tCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtgg                                                                                                               <  1:169010/149‑1 (MQ=255)
                                      cccGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAgtggt                                                                                                              >  1:1079702/1‑149 (MQ=255)
                                        cGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtga                                                                                                            >  1:35671/1‑149 (MQ=255)
                                         ggTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTGGCGGTGCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                                                                           >  1:33625/1‑149 (MQ=255)
                                         ggTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGtgat                                                                                                           <  1:1047894/149‑1 (MQ=255)
                                                ggTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCGGAGGAGTGGTGATGAAAAAt                                                                                                    >  1:401091/1‑149 (MQ=255)
                                                  tGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAACCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTCTCGGTGTTACTTTTCAGAGGAGTGGTGATGAAAAAttt                                                                                                  <  1:659162/149‑1 (MQ=255)
                                                   gCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAAtttt                                                                                                 <  2:1042695/149‑1 (MQ=255)
                                                    cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                                                                <  1:337681/149‑1 (MQ=255)
                                                    cTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTg                                                                                                >  2:74567/1‑149 (MQ=255)
                                                      gtgAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                                                                              <  1:231250/149‑1 (MQ=255)
                                                      gtgAGTTCAGAAAGATTCTGCAATCTATGGGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                                                                              <  1:359199/149‑1 (MQ=255)
                                                      gtgAGTGCAGAAAGATTCTGCAATCTATGTGCATGGAAAAAATGTTCAGATAGATATCATTCTCCTGAACAGTTAGAGGAACGGATAAACCGTGCGGGCCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAACAATTTTGaa                                                                                              >  2:6288/1‑149 (MQ=255)
                                                      gtgAGGTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGaa                                                                                              <  1:1015658/149‑1 (MQ=255)
                                                        gAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGt                                                                                            <  1:286695/149‑1 (MQ=255)
                                                             cAGAAAGATTCTGCAATCTATGGCCCTTTAAAAAACCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGTACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTTAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGc                                                                                       <  2:782174/149‑1 (MQ=255)
                                                              aGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGCTGAAAAATTTTGAAGTTCTGCa                                                                                      >  1:455958/1‑149 (MQ=255)
                                                                      tCTGCAATCTATGTGCATTGAAAAAATCTTCAGTTAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGGGGTGATGAAAACTTTTGAAGTTCTGCAGCCACGAc                                                                              >  2:1306442/1‑149 (MQ=255)
                                                                          cgaTCGATGTGCAGTGAAAAAATCTTCAGATAGATATCATTCTCCTGCACAATCGGATGAACGGCAGAACCTCTCGGTTCTGCGGTTTTTGTTTTAGTTTTCAGAGGAGTGGTGATGAAAAACGTTGACGTTCTGCAGCCACTACatat                                                                          <  2:529129/147‑4 (MQ=255)
                                                                             tCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAAACTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTGCAGAACTc                                                                       >  1:664100/1‑149 (MQ=255)
                                                                              cTATGTGCATTGAAAAAATCTTCATATAGATATCATTATCCTGAACAACCCGATGAACGGCAGAACCTCGCGGTTCCGCCGTTTTGGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTcc                                                                      <  1:1135638/149‑1 (MQ=255)
                                                                                aTGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCt                                                                    <  1:261931/149‑1 (MQ=255)
                                                                                   tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACCATCCGCTGAACGGCAAGACCTCGCGGTTCGTCCGTTTTTGTTTTACTTTTCAGAGGAGCGGTGCTGGCACAGTTTGTAGTTCTTCAGCCACTACACACCTCCCTTTc                                                                 >  2:704884/1‑149 (MQ=255)
                                                                                   tgCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGAGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTc                                                                 >  1:1057674/1‑149 (MQ=255)
                                                                                    gCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTcc                                                                <  1:146850/149‑1 (MQ=255)
                                                                                       ttGAAACAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTATGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCCCTACAGAACTCCCTTTCCGGt                                                             >  2:834406/1‑149 (MQ=255)
                                                                                       ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGt                                                             <  2:1002712/149‑1 (MQ=255)
                                                                                       ttGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTGCTGCAGCCACTACAGAACTCCCTTTCCGGt                                                             >  1:1240107/1‑149 (MQ=255)
                                                                                        tGAAAAAATCTTCAGATAGATATCAGTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTc                                                            <  1:823891/149‑1 (MQ=255)
                                                                                         gAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTCCAGAACTCCCTTTCCGGTCt                                                           >  1:664032/1‑149 (MQ=255)
                                                                                          aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCtc                                                          >  2:374935/1‑148 (MQ=255)
                                                                                          aaaaaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTa                                                          >  2:1182609/1‑149 (MQ=255)
                                                                                             aaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCg                                                       <  2:1207306/149‑1 (MQ=255)
                                                                                             aaaTCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAAACTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTCCAGAACTCCCTTTCCGGTCttacg                                                       >  1:340169/1‑145 (MQ=255)
                                                                                                     agatagatATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGt                                               <  2:128253/149‑1 (MQ=255)
                                                                                                     agatagatATCATTCTCCTGAACAATCAGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTGTTCCGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTACGGGCTCCCCGTATGGGt                                               >  1:43750/1‑149 (MQ=255)
                                                                                                     agatagatATCATTCTCCGGAACAATCCGATGAACGGCAGAACCGCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGt                                               <  1:408140/149‑1 (MQ=255)
                                                                                                       atagataTCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTAt                                             <  1:1037366/149‑1 (MQ=255)
                                                                                                                    tcCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGGTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATAtt                                <  1:783385/149‑1 (MQ=255)
                                                                                                                      cTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAAGGTATATTgc                              >  1:524616/1‑149 (MQ=255)
                                                                                                                       tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCCGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCACTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgca                             >  1:602341/1‑149 (MQ=255)
                                                                                                                       tGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgca                             <  2:1035229/149‑1 (MQ=255)
                                                                                                                          aCAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTgcagca                          <  1:1042451/149‑1 (MQ=255)
                                                                                                                           cAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGcagcag                         >  2:1179440/1‑149 (MQ=255)
                                                                                                                              tCCGCTGAACGGCAGAACCTCGCGGTTCGGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCTCTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATa                      <  1:840672/149‑1 (MQ=255)
                                                                                                                               ccGATGAACGGCAGAGCCTCGCGGTTCTGCCGTTTTTGTTTTATTGTTCAGGGGAGTGGTGATGAAATATGTGGGATTTCTGCAGCCCCTACAGAACCCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATaa                     <  2:401091/149‑1 (MQ=255)
                                                                                                                                cGATGAACGGCAGAACCTCGCGGTTCTGCCGGTGTTGTGTTACTTTTGAGAGGAGAGGGGGTGAAAAAGTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTAGCTGAACGTATATTGCAGCAGATaaa                    >  2:658921/1‑149 (MQ=255)
                                                                                                                                cGATGAACGGCAGAACCTCGCGGTTCAGAGCTGGCGCTGGTACTTTTCAGAGGAGTGGTGATGAAAAATTGTGAAGTTCTGCAGCAACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGGATATTGCCGCAGAGaaa                    >  1:661003/1‑149 (MQ=255)
                                                                                                                                   tGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCa                 <  1:1105612/149‑1 (MQ=255)
                                                                                                                                    taaCGGCAGCACGTCGCGGTTCTGCCGTTTTTGTTTTACGTTTCAGAGGAGTGGTGATGGAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGATCGTATATTGCAGCAGATAAATCAg                <  2:602665/148‑1 (MQ=255)
                                                                                                                                      aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGtt              >  1:331628/1‑149 (MQ=255)
                                                                                                                                      aCGGCAGAACCGCGAGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGAAGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGCTAAATCAGtt              >  2:862992/1‑149 (MQ=255)
                                                                                                                                        ggCAGAACCTCGCGGTTCTGCCGTTTTTTTTTTTTTTTTCAGAGGAAGGGTGATGAAAAATTTTGAACTTCTGCAGCCACTACAGACCTCCCTTTCCCGTCTACCGCTAGGGGTATCTGACCGTATATTCCAGCCGGTCAATCAGTTaa            >  2:245224/1‑149 (MQ=255)
                                                                                                                                        ggCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCCATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTaa            >  2:935665/1‑149 (MQ=255)
                                                                                                                                                cTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGCGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAAcccccgac    >  2:187633/1‑144 (MQ=255)
                                                                                                                                                 tCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACg   <  2:664032/149‑1 (MQ=255)
                                                                                                                                                  cgcgGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGa  <  1:843117/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCATTCTCCTGAACAATCCGATGAACGGCAGAACCTCGCGGTTCTGCCGTTTTTGTTTTACTTTTCAGAGGAGTGGTGATGAAAAATTTTGAAGTTCTGCAGCCACTACAGAACTCCCTTTCCGGTCTACCGCTATGGGTATCTGAACGTATATTGCAGCAGATAAATCAGTTAACCCACTACGA  >  NC_000913/2767520‑2767814

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