Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,542 A→T 100% 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,5420AT100.0% 218.5 / ‑6.7 66Y67N (TAT→AAT) ydaCDUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (44/20);  minor base C (0/2);  total (44/22)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTC  >  NC_000913/1413404‑1413687
                                                                                                                                          |                                                                                                                                                 
tCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCCATGCAAtt                                                                                                                                         >  2:652346/1‑149 (MQ=255)
     ccGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCgg                                                                                                                                    >  2:859750/1‑148 (MQ=255)
       gTATTCAACCATCCACCCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAgg                                                                                                                                  <  1:449838/149‑1 (MQ=255)
          ttCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATCCAACTCTGGAGGACt                                                                                                                               >  2:1113264/1‑149 (MQ=255)
            cAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGa                                                                                                                             >  2:639685/1‑149 (MQ=255)
             aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTGTTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGCCTGaa                                                                                                                            >  1:233077/1‑149 (MQ=255)
             aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGATGCCTGaa                                                                                                                            >  1:139485/1‑149 (MQ=255)
             aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                            <  2:363585/149‑1 (MQ=255)
             aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                            >  1:864809/1‑149 (MQ=255)
             aaCCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGaa                                                                                                                            >  1:45158/1‑149 (MQ=255)
                    cACTCTTCGTTAAAGACTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                     <  1:89363/149‑1 (MQ=255)
                    cACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTtctc                                                                                                                     <  1:696298/149‑1 (MQ=255)
                         ttCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                                >  1:992657/1‑149 (MQ=255)
                         ttCATTAAAGATTACTTGTGCCATCGCATCAACTCTGGCATTCTCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAAt                                                                                                                >  1:1139300/1‑149 (MQ=255)
                             ttAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTaaa                                                                                                            >  2:852905/1‑149 (MQ=255)
                                acgATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAAtt                                                                                                         <  2:75190/147‑1 (MQ=255)
                                           tgCCATCGCATCCCCTCTGGCAGGCGCCAATGTTAGCCGGAAATTGATGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTa                                                                                              <  1:145790/149‑1 (MQ=255)
                                             ccATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACt                                                                                            <  2:602673/149‑1 (MQ=255)
                                             ccATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGATATGAAGTTTTAATTTGAATTCAATGCAATTCTCGGGGACGGACGTTTCTCGCAATTAAAATTTGTCAAGTTTACt                                                                                            >  2:966661/1‑149 (MQ=255)
                                              cATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACtt                                                                                           >  1:295380/1‑149 (MQ=255)
                                               aTCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTAACAGTTTTACttt                                                                                          >  2:831236/1‑149 (MQ=255)
                                                 cGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCt                                                                                        <  1:517951/149‑1 (MQ=255)
                                                  gCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTg                                                                                       >  1:1063995/1‑149 (MQ=255)
                                                    aTCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGct                                                                                     >  2:976277/1‑149 (MQ=255)
                                                      cACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGATATGTTGGCATGATTTGAATTCAATGCCATCCTCGGGGCCTGAAGTTTCTCGCAATTAAAAGTTCTCAGTTTTACTTTCTGctct                                                                                   >  2:291745/1‑149 (MQ=255)
                                                        ccTCTGGCAGGCGCCAATGTTAAACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctctct                                                                                 <  2:148031/149‑1 (MQ=255)
                                                         ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                                >  1:697985/1‑149 (MQ=255)
                                                         ctctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                                >  2:369910/1‑149 (MQ=255)
                                                           ctGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGa                                                                              >  2:912639/1‑149 (MQ=255)
                                                                 ggCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGcc                                                                        >  1:1074417/1‑149 (MQ=255)
                                                                    gcCAATGTTAGACTGAAATTGACGCCCGATTCTGATTATTAATAATCAGCTATGAAGTTTTAATTTTAATTCAATGTAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                     <  2:328953/149‑1 (MQ=255)
                                                                    gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCACTTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                     >  2:352393/1‑149 (MQ=255)
                                                                    gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                     <  1:463516/149‑1 (MQ=255)
                                                                    gcCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGATTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                     >  2:352396/1‑149 (MQ=255)
                                                                      cAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGctt                                                                   >  1:456188/1‑149 (MQ=255)
                                                                           ttAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTAACTTTCTGCTTTCTTGCAACGCCTGCTTCtttt                                                              >  2:139533/1‑149 (MQ=255)
                                                                                 tGAAATTGACGCCCGATGTTGATTATTAATAATGAGCTTTGAGGTTTTAATTTGAATTCAATGCAATTCTCGAGGATTGAATTTTCTCGCGATTCAAAATTGTCTGGTTTACTTTCTGCTCTCTGGACACGCCTGCTTCTTTTTTACCt                                                        >  2:862406/1‑149 (MQ=255)
                                                                                 tGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCt                                                        >  2:232241/1‑149 (MQ=255)
                                                                                   aaaTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTga                                                      <  2:1139300/149‑1 (MQ=255)
                                                                                     agtggctcccGATGTTTATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTgaga                                                    <  2:401281/142‑1 (MQ=255)
                                                                                     aTTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTgaga                                                    >  1:1019270/1‑149 (MQ=255)
                                                                                               gATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                          <  1:1003838/149‑1 (MQ=255)
                                                                                                atttttattattAATCATCAGCTCTGACTGTTTAACTTTAATCCAATGCATTTCTCGAGTACTGATGTTCCTCGCAAGTAAAATTTATCAGTTTTACTTTCTTCTCTCCGGTAAGTCCTTCTTCTTTTTTGCCTGAGAGCCTTTTTTCg                                         <  2:983532/143‑1 (MQ=255)
                                                                                                   ttGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTATTTCGCAt                                      >  2:801859/1‑149 (MQ=255)
                                                                                                    tGATTATTAATATTCCGCGTTGAAGTTTTAATTTGAATCCAATGCAATTCGCGAGGACTGAAGTCTCTCGCAATGAAAATTTATCAGTTTTATTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCAtt                                     <  2:338963/149‑1 (MQ=255)
                                                                                                    tGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCAtt                                     >  1:909922/1‑149 (MQ=255)
                                                                                                       ttattaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTg                                  >  1:209550/1‑149 (MQ=255)
                                                                                                         attaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGGAGTTTCTCGCAATTAAAATTTATCAGTTCTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTGTTACCTGAGACCATTTTTTCGCTTGCTgct                                >  2:680290/1‑147 (MQ=255)
                                                                                                          ttaATAATCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTTAtt                               >  1:664370/1‑149 (MQ=255)
                                                                                                           taataaTCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACTTGAGAGCCTTTTTTCGCATTCTGAttt                              >  2:1020264/1‑149 (MQ=255)
                                                                                                            aataatCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCGCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTc                             <  1:461878/149‑1 (MQ=255)
                                                                                                              taatCAGCTATGAAGTTTCAATTTGAATTCAATGCACTTCTCGAGGACTGAAGTTTCTCGCAATTGAAATTTATCAGTCATACTTTCTGCTCTCTGGAAACGCCTGCTTCGTTTTTACCTGAGAGCATTTTTTCGCCTTCCGATTTCGt                           >  2:209337/1‑149 (MQ=255)
                                                                                                                atCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTTGCTTTCTGTTTTCGTTa                         >  2:307075/1‑149 (MQ=255)
                                                                                                                 tCAGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCACTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAg                        <  2:849967/149‑1 (MQ=255)
                                                                                                                   aGCTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGATTTGCTTTCTGCTCTCTGGAAACGCCTGATTCTTTTTTGGCGGAGCGCGTTTTTTCGCATTCTGATTTCGTTAGtt                      >  1:1099713/1‑149 (MQ=255)
                                                                                                                     cTATGAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTa                    >  2:238735/1‑149 (MQ=255)
                                                                                                                         gAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGTATTCTGATTTCGTTAGTTTAGAtt                >  2:242648/1‑149 (MQ=255)
                                                                                                                         gAAGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAtt                >  1:267373/1‑149 (MQ=255)
                                                                                                                          aaGTTTTAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAttt               >  2:588504/1‑149 (MQ=255)
                                                                                                                           aGTTTTAATTTGCATTCAATGCAATTCTCGAGGACTTAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAtttt              <  2:139485/149‑1 (MQ=255)
                                                                                                                           aGTTTTAAGTTGAATTCAATGCAATTCTCGAGGACTGAATTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCTTGCTTCTTTGTTACCTGGAAGCATTTTTTCGCATTCTGATTTCGTTATTTTAGAtttt              <  2:233077/149‑1 (MQ=255)
                                                                                                                             ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGa            <  1:801859/149‑1 (MQ=255)
                                                                                                                             ttttAATTTGAATTCAATGCAATTCTCGAGGACTGAAGTTTATCGCAATTAAAATTTATCAGTTTTACTTTCTGCTATCTGGAAAAGCATGCTTCCTTTTTACTTGAGAGCATTTTTTCGCATTCCGATTTCGTTAGTTTAGATTTCGa            >  2:32550/1‑149 (MQ=255)
                                                                                                                              tttAAGTTGAATTCAATGCGATTCTCGAGGGCTGAAGTTTCTTGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGTAACGCCTGCTTCTTTTTTTCCTGAGAGCATTTTTTCGTATTCTGATTTCGTTAGTTTAGATTTTGaa           <  2:377135/149‑1 (MQ=255)
                                                                                                                                    ttGAATTCAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCtt     >  1:479002/1‑149 (MQ=255)
                                                                                                                                       aaTTCAATGCAATTCTCGAGGACTGAAGTTTCTCGTAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCATGCTTCTTTTTTACCTGAGAGCCTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGCATATCTTGTc  >  1:281209/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
TCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTC  >  NC_000913/1413404‑1413687

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: ‑