Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,317 A→C 5.2% N21H (AAT→CAT)  yffL → CPZ‑55 prophage; uncharacterized protein YffL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3170AC 5.2% 219.9 / 3.2 77N21H (AAT→CAT) yffLCPZ‑55 prophage; uncharacterized protein YffL
Reads supporting (aligned to +/- strand):  ref base A (31/42);  new base C (2/2);  total (33/44)
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.71e-01

TCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGC  >  NC_000913/2560178‑2560456
                                                                                                                                            |                                                                                                                                           
tCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtg                                                                                                                                          <  2:1067110/143‑1 (MQ=255)
tCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtg                                                                                                                                          >  1:341048/1‑143 (MQ=255)
  cGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                                                                        <  1:165710/143‑1 (MQ=255)
    gCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgttgc                                                                                                                                      <  2:858655/143‑1 (MQ=255)
    gCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                                                      <  1:103503/143‑1 (MQ=255)
     cttttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                                                     >  1:899271/1‑143 (MQ=255)
      ttttcttttcCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                                                    <  1:917462/143‑1 (MQ=255)
      ttttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCACCGAAACATCCAGCGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGCAAcgccccaa                                                                                                                                    >  2:436908/1‑138 (MQ=255)
      ttttcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                                                    >  2:350595/1‑143 (MQ=255)
      tttgcttctTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                                                    <  2:1045911/143‑1 (MQ=255)
       ttccttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                                                                   <  1:956584/140‑1 (MQ=255)
         tcttttTCGCTAATTATTCCATTATTCAAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACAt                                                                                                                                 >  1:317157/1‑143 (MQ=255)
           tttttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                                                               >  2:730663/1‑143 (MQ=255)
             ctacGCTAATTATTCCATTATTCCAATTAAGGGGAATAAATTATCCATACTGGGAATCCTTTTTCGATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCg                                                                                                                             <  1:510174/140‑1 (MQ=255)
               tCGCTAATTATTCCATTATTCCAATTGAGTGGAATAAGTTGTTCAATAGCAATAACTTAAATCCATGTCTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                                                           <  1:249394/143‑1 (MQ=255)
               tCGCGAATTATTCCATTATTCCACTTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGAGCAACACAAATAAAAGCACCCGGAATGCTGCAACATGCCggg                                                                                                                           <  2:197241/143‑1 (MQ=255)
                cGCTACTAATTACAGTATTCCAATTAAGTGGAATAAATTATTCAATAGCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                                                          <  1:406198/143‑1 (MQ=255)
                cGCTAATTATTCCATTATTCTAATTAAGTGGAATAAATTAGTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                                                          <  1:1011840/143‑1 (MQ=255)
                  cTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                                                        >  1:201558/1‑143 (MQ=255)
                     gtgagtccagtgttccaAGTACGTGGTATAAAGTAGTCACGAGCAAGAACTGAAGTCCGGGTTTGCAACGAAAGATCCAGAGGTCGAAAAGAGGATCAACACAAATAAAAGCCGCCGGAATGGTGCAACATGCCGGGGGCGGt                                                                                                                     <  2:68884/138‑1 (MQ=255)
                     attattccattattccAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                                                     >  1:289443/1‑143 (MQ=255)
                       tattccattattccaATTAAGTGCAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCCACATGCCGGGGGCGGTAt                                                                                                                   <  1:319026/143‑1 (MQ=255)
                       tattccattatcgcaATTACGGAGAATAAATTATTCAATATCAAGAACTGAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                                                                   <  1:919151/143‑1 (MQ=255)
                         ttccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTGTga                                                                                                                 <  2:965471/143‑1 (MQ=255)
                         ttccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                                                 <  1:325831/143‑1 (MQ=255)
                          tccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAATCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                                                >  1:1010260/1‑143 (MQ=255)
                          tccattagtccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                                                <  1:970216/143‑1 (MQ=255)
                            cattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                                                              <  2:261705/143‑1 (MQ=255)
                                  tccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAACTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                                                        <  1:773256/143‑1 (MQ=255)
                                   ccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                                                       <  2:500964/143‑1 (MQ=255)
                                      atcacgTGGAATACACTACTCAACACCAAGAACCCAAAGCCCTGTTTCCAACGAAAGATCCAGAGTTCCAAAATAGGATCAACACACATAACAGCCCCCGGAATGCTGCCACATGCCGGGTGCGGTATGAGAAACAAGCTaaa                                                                                                    <  1:526740/138‑1 (MQ=255)
                                       ttAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGGCCGATCCACAACTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                                                   <  1:806859/143‑1 (MQ=255)
                                           gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                                                  >  1:172328/1‑124 (MQ=255)
                                           gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                                                                  <  2:172328/124‑1 (MQ=255)
                                           gTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCACCCGGAAAGCTGCAACATGCCGGGGGCGGCATGACAAACAAGCCAAAGgaga                                                                                               >  2:232997/1‑143 (MQ=255)
                                               aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGCTATGAGAAACAAGCTAAAGGAGAAttt                                                                                           >  1:457335/1‑143 (MQ=255)
                                                   acttatgCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                                                       <  2:109893/141‑1 (MQ=255)
                                                    tgtatgCAATAGCAAGAACTTAAATCCACGTTTGCAACGAACGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                                                      <  2:272397/141‑1 (MQ=255)
                                                       attCAATATCAAGAACTTAAATCCCTCATCGCAACGAAAGATCCAGAGTTCGATAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                                                   <  1:748086/143‑1 (MQ=255)
                                                       attCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                                                   <  1:72300/143‑1 (MQ=255)
                                                        ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                                                                  <  1:267010/143‑1 (MQ=255)
                                                         tcaagatcaagaACTTAAATCCACGTTTGCAACGAAAGATCCAGAGTTCGAAAAGAGGATCAACACAAATAAAAGACCCCCGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                                                 <  1:445528/143‑1 (MQ=255)
                                                         tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAGTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                                                 >  2:545262/1‑143 (MQ=255)
                                                            atatCAAGAACTTAAAGCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                                                              <  1:950128/143‑1 (MQ=255)
                                                               tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGATATGAGCAACAAGCTAACGGAGAATTTCTCATGCCGGATATGc                                                                           >  2:928741/1‑143 (MQ=255)
                                                               tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATAAACACAAATAAAACCCCCCGGAATGCTGCAACATGGCGGGGGCGCGAACAGAAACCCGCCAAACGCGAGTGTCTCATGGCTGctactc                                                                           >  1:1093147/1‑140 (MQ=255)
                                                                 aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                                                         >  2:812417/1‑143 (MQ=255)
                                                                 aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                                                         >  2:781511/1‑143 (MQ=255)
                                                                  aGATCTTACATCCATGTTCGCAACGAAAGCTCCAGAGTTCTAAAATAGGATCAACACCAATAAAAGCGCCCGGAATGCTGCGACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                                                        <  2:186008/143‑1 (MQ=255)
                                                                  aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                                                        >  1:766109/1‑143 (MQ=255)
                                                                  aGAACTGAAATCCATGGTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAAGAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                                                        <  1:445322/143‑1 (MQ=255)
                                                                        ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTCGCAGTa                                                                   >  1:640913/1‑143 (MQ=255)
                                                                        ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAACTTCTCATGTCTGAAATGCTAGCAGTa                                                                   >  1:205334/1‑143 (MQ=255)
                                                                         tAAATCCATGTTTGCAGCGAAAGATCCAGAGTTGGAAAATAGAATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                                                                  <  2:906606/143‑1 (MQ=255)
                                                                           aaTCGATGTTGGCAAGAAAAGATCGAGAGTTAGAGAATAGGATAAATACGAATAAAAGACCCCGGAATGCTGCAACATGCCGGGGGCGTTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATAGGCTAGCAGTAGaa                                                                <  2:748039/143‑1 (MQ=255)
                                                                           aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                                                                >  1:763474/1‑143 (MQ=255)
                                                                                aTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                                                           <  1:737925/143‑1 (MQ=255)
                                                                                  gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                                                         <  2:341048/143‑1 (MQ=255)
                                                                                        aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                                                   >  1:902431/1‑143 (MQ=255)
                                                                                        aaCGAAAGAACCAGAGTCCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                                                   <  1:299308/143‑1 (MQ=255)
                                                                                           gAAAGATCCAGAGTTCGAAAATAGGATCAAAACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGAGCAGTATGAGAAACAAACTAAAAGAGAATTTCTCATGTCTGATATGATAACAACAGAAAAAGAAACAAACAatg                                                >  2:893308/1‑143 (MQ=255)
                                                                                                  ccAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                                                      <  1:535512/130‑1 (MQ=255)
                                                                                                  ccAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                                                      >  2:535512/1‑130 (MQ=255)
                                                                                                        ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtt                                   >  2:173629/1‑143 (MQ=255)
                                                                                                                    aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAAt                       >  1:518596/1‑143 (MQ=255)
                                                                                                                        acacAGATAAAAGCCCCCGGAAGGCTGCAAGATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAAc                   <  2:436871/143‑1 (MQ=255)
                                                                                                                        acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAAc                   >  1:893920/1‑143 (MQ=255)
                                                                                                                          acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAGTTCTCGCACAAAATTAACgc                 >  1:952891/1‑142 (MQ=255)
                                                                                                                            aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGt               >  1:876634/1‑143 (MQ=255)
                                                                                                                               tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGAGATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTAc            <  1:545262/143‑1 (MQ=255)
                                                                                                                                 accgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACgg          >  1:781044/4‑143 (MQ=255)
                                                                                                                                      cgccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGAAATGATAGGAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATaa     <  2:727314/141‑1 (MQ=255)
                                                                                                                                      ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATaa     >  2:223086/1‑143 (MQ=255)
                                                                                                                                       cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATcaa    >  2:344887/1‑143 (MQ=255)
                                                                                                                                       cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATaaa    <  2:895294/143‑1 (MQ=255)
                                                                                                                                        ccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCCCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTCTAAACAAAATTAACGAGTTACAGAATAAAg   >  1:976437/1‑143 (MQ=255)
                                                                                                                                         cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGc  <  1:355865/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTTTAAACAAAATTAACGAGTTACGGAATAAAGC  >  NC_000913/2560178‑2560456

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