Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,413,542 A→T Y67N (TAT→AAT)  ydaC ← DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,5420AT98.4% 211.1 / NA 64Y67N (TAT→AAT) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (31/32);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATC  >  NC_000913/1413402‑1413682
                                                                                                                                            |                                                                                                                                             
cATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCaa                                                                                                                                       <  1:946310/149‑1 (MQ=255)
cATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCaa                                                                                                                                       >  2:902079/1‑149 (MQ=255)
    aGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATCATTAATAATCAGCTATGCAGTTTTAATTTGAATTCAATGCAATtct                                                                                                                                   >  2:483332/1‑149 (MQ=255)
       ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGa                                                                                                                                <  2:483705/149‑1 (MQ=255)
           aTTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGAc                                                                                                                            >  2:784372/1‑149 (MQ=255)
               aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCGGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGGTTTAATTTGAATTCAATGCAATTCTCGAGGACTGag                                                                                                                        >  2:693287/1‑148 (MQ=255)
               aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                        <  2:737322/149‑1 (MQ=255)
                aCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAg                                                                                                                       >  1:823137/1‑149 (MQ=255)
                   aTCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCCGGCGCCCAAGTTACAGACAAATTGACAACCCGTGGTGATTATTAATAAGCAGCAATGAAGTGGTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGttt                                                                                                                    >  1:257479/1‑149 (MQ=255)
                    tcccccctTAATTAGAGATTATCTTTGCCATCGCATCACCTCTGGCAGGGGCCAATTGTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtc                                                                                                                   <  2:752470/143‑1 (MQ=255)
                    tCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtc                                                                                                                   <  2:126929/149‑1 (MQ=255)
                     ccACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtct                                                                                                                  >  1:814369/1‑149 (MQ=255)
                      cACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                 <  1:952403/149‑1 (MQ=255)
                      cACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                 <  2:504418/149‑1 (MQ=255)
                          ctTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCaa                                                                                                             >  2:705127/1‑149 (MQ=255)
                           ttCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                            >  1:327813/1‑149 (MQ=255)
                                ttaaGATTGCTTTTGCCATCGCATCACCTCTGGCATTCTCCAATGTTATACTGAAATTGACGCCCTATGTTTATTATTAAGAAGCAGCTAGGAAGTTTGAATTTGAATTCAATGCAAGTCTCGAGGACTGAAGTTTCTCGCAATTaaaa                                                                                                       <  1:495511/147‑1 (MQ=255)
                                tAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGCTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaaa                                                                                                       <  2:996384/149‑1 (MQ=255)
                                 aaaGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAt                                                                                                      >  1:581838/1‑149 (MQ=255)
                                   aGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAttt                                                                                                    <  1:428919/149‑1 (MQ=255)
                                       tACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCa                                                                                               <  1:797395/149‑1 (MQ=255)
                                          ttGTGCCATCGCATCCCATCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGtt                                                                                            <  1:788496/149‑1 (MQ=255)
                                          ttGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGtt                                                                                            >  2:947452/1‑149 (MQ=255)
                                            gtgCCATCGCATCGCCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATGATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGGGGACTGAAGTTTGTCCCAATTAAAATTTACCGGtttg                                                                                          >  2:1011676/1‑148 (MQ=255)
                                                  tCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTc                                                                                    <  1:527003/149‑1 (MQ=255)
                                                      aTCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTATTAATCGGCTTTGAGGTTTTATTTTGAATTCATTGCATTTTCCGAGGACGGAAGTTTCTGGCAATTAAATTTATCCATTTTTACTTTTTGct                                                                                >  1:1051889/1‑149 (MQ=255)
                                                       tCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctc                                                                               >  1:481757/1‑149 (MQ=255)
                                                           ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                           <  1:721807/149‑1 (MQ=255)
                                                           ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                           >  1:777802/1‑149 (MQ=255)
                                                           ctcgGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGCGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                           >  2:287912/1‑149 (MQ=255)
                                                              tGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTggca                                                                        >  1:716504/1‑147 (MQ=255)
                                                               ggCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGaaa                                                                       <  2:1003704/149‑1 (MQ=255)
                                                                gCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAAc                                                                      >  1:817393/1‑149 (MQ=255)
                                                                gCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAACTTAACAGTTTTACTTTCTGCTCTCTGGAAAc                                                                      >  2:498505/1‑149 (MQ=255)
                                                                 cAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCCGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAATTATCAGTTGTTCTTTCTGCTCTCTGGAAACg                                                                     >  2:810930/1‑149 (MQ=255)
                                                                 cAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACg                                                                     >  2:760569/1‑149 (MQ=255)
                                                                  aGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGc                                                                    <  1:1150/149‑1 (MQ=255)
                                                                  aGGCGCCAATGTTAGACTGAAATGGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGc                                                                    <  1:123026/149‑1 (MQ=255)
                                                                  aGGCGCAAATGTTAGACTGAAATTGCCGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGc                                                                    <  2:614353/149‑1 (MQ=255)
                                                                   ggCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGcc                                                                   <  1:963084/149‑1 (MQ=255)
                                                                    gcgcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCt                                                                  >  1:361772/1‑149 (MQ=255)
                                                                         aaagTTAAACTGAAATTGACCGCCGATGTTGATTATTAATAATTAGCTATTAAGCTTTAATTTGAAGACAATGCAATTCTCCATGACTGAAGGATCGCGCAATTAAAATTTATCAGTCTTACTGTCTGCTCTCTGGAACCGCCTGcttc                                                             <  2:124861/146‑1 (MQ=255)
                                                                          aTGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAGTTCATCAGTTATACTATCTGCTCTCTGGAAACGCCTGcttct                                                            >  1:417165/1‑149 (MQ=255)
                                                                            gTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGATTCAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAGTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAGACGTCTgacccgat                                                          >  2:654992/1‑142 (MQ=255)
                                                                                  cTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCTCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAcc                                                    <  1:678308/149‑1 (MQ=255)
                                                                                     aaaTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTga                                                 <  1:784372/149‑1 (MQ=255)
                                                                                           actcccGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAt                                           <  2:83725/146‑1 (MQ=255)
                                                                                             gCCCGATGTTGATTCTTAACAAGCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAttt                                         <  2:544434/149‑1 (MQ=255)
                                                                                             gCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGACGTTTCTCGCAATTAAGATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAgtt                                         >  2:359846/1‑149 (MQ=255)
                                                                                             gCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAttt                                         >  2:1075260/1‑149 (MQ=255)
                                                                                                 gATGTTGATTATTAATAATCAGCTATGCAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                     <  1:761905/149‑1 (MQ=255)
                                                                                                       gATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTgg                                                                          <  2:512559/106‑1 (MQ=255)
                                                                                                       gATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTgg                                                                          >  1:512559/1‑106 (MQ=255)
                                                                                                         ttattaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTg                             <  1:464233/149‑1 (MQ=255)
                                                                                                                     aGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGtt                 <  1:381599/149‑1 (MQ=255)
                                                                                                                       cTATGCAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATGAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTa               <  2:74562/149‑1 (MQ=255)
                                                                                                                         aTGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGGATTCTGCTTTCGTTAGTTTAGa             >  2:823560/1‑149 (MQ=255)
                                                                                                                           gAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCCTTCTGATTTCGTTAGTTTAGCtt           >  2:571251/1‑149 (MQ=255)
                                                                                                                             aGTTTTAATGTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAtttt         <  1:1064317/149‑1 (MQ=255)
                                                                                                                                tttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGaa      <  1:756152/149‑1 (MQ=255)
                                                                                                                                tttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGaa      <  2:1019933/149‑1 (MQ=255)
                                                                                                                                tttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGaa      >  2:219356/1‑149 (MQ=255)
                                                                                                                                    aTTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATc  >  1:573878/1‑149 (MQ=255)
                                                                                                                                    aTTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATc  <  1:405377/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
CATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATC  >  NC_000913/1413402‑1413682

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