Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,574 C→T 7.6% T41M (ACG→ATG) ‡ ypjF → CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,5740CT 7.6% 456.9 / 18.4 145T41M (ACG→ATG) ‡ypjFCP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (55/79);  new base T (7/4);  total (62/83)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 149 total aligned reads displayed.

ACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTT  >  NC_000913/2777426‑2777720
                                                                                                                                                    |                                                                                                                                                  
aCGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGacac                                                                                                                                                    <  3:857629/149‑1 (MQ=255)
  gCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGc                                                                                                                                                  <  1:542252/149‑1 (MQ=255)
    taactaaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCg                                                                                                                                                >  2:664962/1‑149 (MQ=255)
       ctaaTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTc                                                                                                                                             <  3:445295/149‑1 (MQ=255)
          aTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTTGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGt                                                                                                                                          <  1:59006/149‑1 (MQ=255)
          aTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGt                                                                                                                                          <  1:251105/149‑1 (MQ=255)
          aTTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAACCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGt                                                                                                                                          <  4:655449/149‑1 (MQ=255)
           tttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtg                                                                                                                                         >  4:69999/1‑149 (MQ=255)
            ttCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtga                                                                                                                                        <  3:723037/149‑1 (MQ=255)
             tCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGtgat                                                                                                                                       <  2:264650/149‑1 (MQ=255)
                  agagCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAAc                                                                                                                                  <  4:33453/149‑1 (MQ=255)
                   gagCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCCGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGTAACt                                                                                                                                 >  4:714207/1‑149 (MQ=255)
                     gcaagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGt                                                                                                                               <  4:446276/149‑1 (MQ=255)
                       aagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                             >  2:350341/1‑149 (MQ=255)
                       aagcaTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtta                                                                                                                             >  2:939643/1‑149 (MQ=255)
                         gcaTGAACGCTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGttatt                                                                                                                           >  3:966595/1‑149 (MQ=255)
                         gcaTGAACACTCTACCTGCTACAATTTCGCAGGCTGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACGGttatt                                                                                                                           >  3:700668/1‑149 (MQ=255)
                           aTGATCACACTACCTGCTACAATTTCGCAGGCGGCGAAGTCCTGCCTATCGCCAGTGGCTGTCTGGCAAATGCTACTCACACGCCTGCTCGAACAACTCTATGGCCTGACACTGAACGACACGTCGCTCAGTGATGAAACTGTTATTaa                                                                                                                         <  4:655036/149‑1 (MQ=255)
                           aTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCTGTGATGAAACTGTTATTaa                                                                                                                         <  4:889851/149‑1 (MQ=255)
                           aTGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTaa                                                                                                                         <  3:677861/149‑1 (MQ=255)
                            tGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAg                                                                                                                        <  3:883003/149‑1 (MQ=255)
                              aaCACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTTACACGCCTGCTGGAACAACACTATGGCCTGTCACTGACCGTCACGCCTTTTAGCGATGAAACTGTTATTCAggt                                                                                                                      <  4:854112/149‑2 (MQ=255)
                               acacTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGaa                                                                                                                     >  3:32443/1‑149 (MQ=255)
                                 acTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACa                                                                                                                   >  1:627919/1‑149 (MQ=255)
                                  ctctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACat                                                                                                                  <  4:735952/149‑1 (MQ=255)
                                   tctACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACata                                                                                                                 >  2:225750/1‑149 (MQ=255)
                                            tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                        <  3:11582/149‑1 (MQ=255)
                                            tACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                        >  1:949695/1‑149 (MQ=255)
                                               aaTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTGCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTa                                                                                                     <  3:298009/149‑1 (MQ=255)
                                                aTTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTAt                                                                                                    <  2:975079/149‑1 (MQ=255)
                                                       aGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg                                                                                             <  3:949991/149‑1 (MQ=255)
                                                       aGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTg                                                                                             >  4:949991/1‑149 (MQ=255)
                                                           ggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                              <  2:155630/112‑1 (MQ=255)
                                                           ggcgAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGtt                                                                                                                              >  1:155630/1‑112 (MQ=255)
                                                              gAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGtgtca                                                                                                                                      >  3:1010479/1‑98 (MQ=38)
                                                              gAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGGCGTGTCGTTCAGtgtca                                                                                                                                      <  4:1010479/101‑4 (MQ=38)
                                                              gAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATg                                                                                      <  1:225750/149‑1 (MQ=255)
                                                               aaGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATgt                                                                                     <  1:389819/149‑2 (MQ=255)
                                                                ccgcctgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCa                                                                                    <  2:401462/146‑1 (MQ=255)
                                                                 gCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAg                                                                                   <  3:638969/149‑1 (MQ=255)
                                                                 gCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAg                                                                                   <  1:1006632/149‑1 (MQ=255)
                                                                     tgcctgTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTTGTATCACTCTGGCCGATGCAGTCaa                                                                               <  4:1036252/149‑1 (MQ=255)
                                                                          gTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTc                                                                          <  2:134695/149‑1 (MQ=255)
                                                                             gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtgg                                                                       <  2:115305/149‑1 (MQ=255)
                                                                             gCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtgg                                                                       >  2:476723/1‑149 (MQ=255)
                                                                              ccAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATTGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      <  2:713307/149‑1 (MQ=255)
                                                                              ccAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      <  4:966595/149‑1 (MQ=255)
                                                                              ccAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      <  2:493913/149‑1 (MQ=255)
                                                                              ccAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCtggt                                                                      <  1:276731/149‑1 (MQ=255)
                                                                                   ggCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGaaa                                                                 <  2:980099/149‑1 (MQ=255)
                                                                                        tCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTAc                                                            <  1:804048/149‑1 (MQ=255)
                                                                                          tGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                          >  2:810348/1‑149 (MQ=255)
                                                                                          tGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                          >  1:975354/1‑149 (MQ=255)
                                                                                          tGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGa                                                          >  1:582429/1‑149 (MQ=255)
                                                                                            gCAAATGCTACTGACACGCCTGCTGGAACAACACTATTGCCGGACACTGAACGACACGCCGTTCAGTGATGAAGCTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                        <  2:776210/149‑1 (MQ=255)
                                                                                            gCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                        >  4:640099/1‑149 (MQ=255)
                                                                                              aaaTGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTg                                                      <  1:350341/149‑1 (MQ=255)
                                                                                                  gatACTGACACGCCTGATGCAACAACACCATGGCATGACACTGAACGACACGCCGTTCAGTGATGAAACTGCTATTAATGAAAATATCGATGCTGGTATCACCCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTAc                                                  <  4:1032756/147‑1 (MQ=255)
                                                                                                   cTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACg                                                 >  2:746599/1‑149 (MQ=255)
                                                                                                    tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAg                                                                                                                        <  3:791047/77‑1 (MQ=255)
                                                                                                    tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAg                                                                                                                        >  4:791047/1‑77 (MQ=255)
                                                                                                    tACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGt                                                >  1:831001/1‑149 (MQ=255)
                                                                                                      cTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACa                                                                                                                   <  1:519439/80‑1 (MQ=255)
                                                                                                       tGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATc                                             >  4:360506/1‑149 (MQ=255)
                                                                                                         acacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                        <  2:167929/88‑1 (MQ=255)
                                                                                                         acacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTg                                                                                                        >  1:167929/1‑88 (MQ=255)
                                                                                                          cacGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGAt                                          <  1:939643/149‑1 (MQ=255)
                                                                                                           acGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATc                                         >  1:38334/1‑149 (MQ=255)
                                                                                                           acGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATc                                         <  2:26616/149‑1 (MQ=255)
                                                                                                             gCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATcac                                       >  3:477674/1‑149 (MQ=255)
                                                                                                                   tgGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGa                                 <  2:831001/149‑1 (MQ=255)
                                                                                                                    ggAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGAt                                                                                       <  2:553579/94‑1 (MQ=255)
                                                                                                                     gAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGAAAAAGTACGAACTGGTACGTATCGAGCACAGAGGAtt                               >  4:778539/1‑149 (MQ=255)
                                                                                                                       acaacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                         <  4:651238/105‑1 (MQ=255)
                                                                                                                       acaacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCt                                                                         >  3:651238/1‑105 (MQ=255)
                                                                                                                         aacaCTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                        >  1:264728/1‑120 (MQ=255)
                                                                                                                           cacTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTg                         <  1:224062/149‑1 (MQ=255)
                                                                                                                            acTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTgg                        <  3:338231/149‑1 (MQ=255)
                                                                                                                             cTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGc                       <  1:810348/149‑1 (MQ=255)
                                                                                                                              tATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                        <  3:1034044/115‑1 (MQ=255)
                                                                                                                              tATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAAc                                                        >  4:1034044/1‑115 (MQ=255)
                                                                                                                                  gCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAAGATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGcaacaa                  <  1:56252/149‑1 (MQ=255)
                                                                                                                                   ccTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCaacaac                 <  2:272410/149‑1 (MQ=255)
                                                                                                                                    cTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAacaaca                <  3:532174/149‑1 (MQ=255)
                                                                                                                                     tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAg               >  1:933133/1‑149 (MQ=255)
                                                                                                                                     tGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAg               >  1:705809/1‑149 (MQ=255)
                                                                                                                                       tcacTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtc             >  2:456779/2‑149 (MQ=255)
                                                                                                                                       acacTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtc             >  2:448393/1‑149 (MQ=255)
                                                                                                                                         acTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGTAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGAt                                          <  4:417673/118‑1 (MQ=255)
                                                                                                                                         acTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGAt                                          >  3:417673/1‑118 (MQ=255)
                                                                                                                                         acTGAACGGCGTGTCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc           >  3:181970/1‑149 (MQ=255)
                                                                                                                                         acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTc                            >  1:872367/1‑132 (MQ=255)
                                                                                                                                         acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc           <  2:933133/149‑1 (MQ=255)
                                                                                                                                         acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGtctc           <  1:1057144/149‑1 (MQ=255)
                                                                                                                                         acTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGAAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTc                            <  2:872367/132‑1 (MQ=255)
                                                                                                                                          cTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTcc          <  3:221470/149‑1 (MQ=255)
                                                                                                                                           tGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCa         <  2:582429/149‑1 (MQ=255)
                                                                                                                                             aaCGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCata       >  1:675577/1‑149 (MQ=255)
                                                                                                                                                  cacGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATAttt  >  4:809326/1‑149 (MQ=255)
                                                                                                                                                   acGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACaga                                    >  1:1017480/1‑114 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTACTGACACGCCTGCTGGAACAACACTATGGCCTGACACTGAACGACACGCCGTTCAGTGATGAAACTGTTATTAAGGAACATATCGATGCTGGTATCACTCTGGCCGATGCAGTCAATTTTCTGGTGGAAAAGTACGAACTGGTACGTATCGATCACAGAGGATTTTCGTGGCAACAACAGTCTCCATATATTT  >  NC_000913/2777426‑2777720

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