Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,200,197 C→A 100% D204Y (GAT→TAT)  intE ← e14 prophage; putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,1970CA100.0% 414.8 / NA 111D204Y (GAT→TAT) intEe14 prophage; putative integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (57/54);  total (57/54)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 112 total aligned reads displayed.

CTCCCGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATCGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACATCAATCAGGACGGAGC  >  NC_000913/1200052‑1200340
                                                                                                                                                 |                                                                                                                                               
cTCCCGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGaa                                                                                                                                              <  3:838105/149‑1 (MQ=255)
 tcacGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAg                                                                                                                                             >  3:614667/4‑149 (MQ=255)
 tCCCGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAg                                                                                                                                             >  2:30451/1‑149 (MQ=255)
  cccGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATagaaga                                                                                                                                            <  1:843023/149‑1 (MQ=255)
   ccGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAt                                                                                                                                           >  4:370737/1‑149 (MQ=255)
   ccGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAt                                                                                                                                           >  4:338132/1‑149 (MQ=255)
   ccGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAt                                                                                                                                           >  1:672509/1‑149 (MQ=255)
    cGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAtt                                                                                                                                          <  4:901313/149‑1 (MQ=255)
       tCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAttttt                                                                                                                                       <  4:435562/149‑1 (MQ=255)
         ttttCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGAttttttt                                                                                                                                     <  3:738894/149‑1 (MQ=255)
              cTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAtt                                                                                                                                <  1:72996/149‑1 (MQ=255)
               tGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               <  2:693031/149‑1 (MQ=255)
                  aTGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttcttc                                                                                                                            <  1:18645/149‑1 (MQ=255)
                  aTGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttcttc                                                                                                                            >  2:540226/1‑149 (MQ=255)
                     aCATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAg                                                                                                                         <  1:487470/149‑1 (MQ=255)
                      cATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAgg                                                                                                                        <  2:826985/149‑1 (MQ=255)
                        tGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCt                                                                                                                      >  4:745555/1‑149 (MQ=255)
                         gAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTc                                                                                                                     <  4:232598/149‑1 (MQ=255)
                           aGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAg                                                                                                                   >  2:511379/1‑149 (MQ=255)
                                     cccAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCgg                                                                                                         <  3:499962/149‑1 (MQ=255)
                                      ccAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCggg                                                                                                        >  3:693453/1‑149 (MQ=255)
                                      ccAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCggg                                                                                                        >  1:437180/1‑149 (MQ=255)
                                       cAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               >  3:921384/1‑125 (MQ=255)
                                        aaaTATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTa                                                                                                      >  4:831964/1‑149 (MQ=255)
                                              cGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAAttttt                                                                                                <  3:475579/149‑1 (MQ=255)
                                               gCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTttc                                                                                               >  3:655957/1‑149 (MQ=255)
                                               gCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTttc                                                                                               >  2:125708/1‑149 (MQ=255)
                                               gCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTttc                                                                                               <  2:437180/149‑1 (MQ=255)
                                                cTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTttct                                                                                              <  3:847436/149‑1 (MQ=255)
                                                cTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTttct                                                                                              >  4:382436/1‑149 (MQ=255)
                                                    aTTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTcttctt                                                                                          <  1:211019/149‑1 (MQ=255)
                                                     tttCATACGCGAAATATCACTTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTg                                                                                         <  3:760382/149‑1 (MQ=255)
                                                     tttCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTg                                                                                         >  3:40005/1‑149 (MQ=255)
                                                             gcgAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               >  3:579845/1‑103 (MQ=255)
                                                             gcgAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               <  4:579845/103‑1 (MQ=255)
                                                              cgAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAg                                                                                <  1:522957/149‑1 (MQ=255)
                                                               gAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGgttgt                                                                                    <  3:444205/144‑1 (MQ=255)
                                                               gAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGgttgt                                                                                    >  4:444205/1‑144 (MQ=255)
                                                               gAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTa                                                                                 >  3:694594/1‑147 (MQ=255)
                                                               gAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTa                                                                                 <  4:694594/147‑1 (MQ=255)
                                                               gAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGt                                                                               <  1:209217/149‑1 (MQ=255)
                                                                   tatCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCt                                                                           >  3:183630/1‑149 (MQ=255)
                                                                    atCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTa                                                                          >  3:894179/1‑149 (MQ=255)
                                                                     tCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTGCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTaa                                                                         <  3:587631/149‑1 (MQ=255)
                                                                       aCCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATg                                                                       >  3:548947/1‑149 (MQ=255)
                                                                          tAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTg                                                                              <  3:757243/139‑1 (MQ=255)
                                                                          tAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTg                                                                              >  4:757243/1‑139 (MQ=255)
                                                                          tAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTg                                                                    >  4:853858/1‑149 (MQ=255)
                                                                           aaCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTgg                                                                   <  3:294734/149‑1 (MQ=255)
                                                                           aaCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTgg                                                                   <  1:283127/149‑1 (MQ=255)
                                                                           aaCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTgg                                                                   <  1:474052/149‑1 (MQ=255)
                                                                           aaCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTgg                                                                   <  1:540226/149‑1 (MQ=255)
                                                                            aCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTggg                                                                  >  4:34998/1‑149 (MQ=255)
                                                                            aCCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTggg                                                                  <  1:117895/149‑1 (MQ=255)
                                                                               gCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGTAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTta                                                               <  4:548947/149‑1 (MQ=255)
                                                                               gCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTACGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTta                                                               >  1:389573/1‑149 (MQ=255)
                                                                                    ccAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGcc                                                          <  3:622329/149‑1 (MQ=255)
                                                                                     cAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCa                                                         >  4:112292/1‑149 (MQ=255)
                                                                                      aGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTc                                                                                                                                    <  4:213559/73‑1 (MQ=255)
                                                                                      aGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCag                                                        <  3:40040/149‑1 (MQ=255)
                                                                                       gTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               >  2:872438/1‑77 (MQ=255)
                                                                                       gTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               <  1:872438/77‑1 (MQ=255)
                                                                                        tAACCAGCGCTAAAAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCagga                                                      <  4:45207/149‑1 (MQ=255)
                                                                                         aaCCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCaggag                                                     <  2:672509/149‑1 (MQ=255)
                                                                                           ccAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTg                                                                              <  2:884142/122‑1 (MQ=255)
                                                                                           ccAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTg                                                                              >  1:884142/1‑122 (MQ=255)
                                                                                           ccAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATCGCCAGCAGGt                                                   >  4:248583/1‑149 (MQ=255)
                                                                                           ccAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGt                                                   >  3:138160/1‑149 (MQ=255)
                                                                                            cAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTa                                                  <  3:71490/149‑1 (MQ=255)
                                                                                             aGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTAc                                                 <  3:745555/149‑1 (MQ=255)
                                                                                             aGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTACGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTAc                                                 <  2:389573/149‑1 (MQ=255)
                                                                                                  tAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTtctc                                            >  1:581219/1‑149 (MQ=255)
                                                                                                   aaCAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTcc                                           <  2:138425/149‑1 (MQ=255)
                                                                                                    aCAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCg                                          >  4:764520/1‑149 (MQ=255)
                                                                                                      aGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTa                                        <  1:223571/149‑1 (MQ=255)
                                                                                                        cATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGt                                      <  4:614667/149‑1 (MQ=255)
                                                                                                         aTGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCt                                                                     >  3:288686/1‑117 (MQ=255)
                                                                                                         aTGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCt                                                                     <  4:288686/117‑1 (MQ=255)
                                                                                                         aTGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGtg                                     >  3:894165/1‑149 (MQ=255)
                                                                                                         aTGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGtg                                     >  3:452584/1‑149 (MQ=255)
                                                                                                            gCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTa                                                  >  4:406049/1‑133 (MQ=255)
                                                                                                            gCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTa                                                  <  3:406049/133‑1 (MQ=255)
                                                                                                            gCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGtgctg                                  >  1:703454/1‑149 (MQ=255)
                                                                                                               tttCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTgcgc                               >  2:920301/1‑149 (MQ=255)
                                                                                                                 tCCCATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCct                             <  3:445276/149‑1 (MQ=255)
                                                                                                                  cccATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               >  3:855629/1‑50 (MQ=255)
                                                                                                                  cccATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCAttc                                                                                                                               <  4:855629/50‑1 (MQ=255)
                                                                                                                    cATGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCtt                          >  4:236320/1‑149 (MQ=255)
                                                                                                                     aTGTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCttt                         >  1:428514/1‑149 (MQ=255)
                                                                                                                       gTAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGa                       <  1:511379/149‑1 (MQ=255)
                                                                                                                        tAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGaa                      >  2:543536/1‑149 (MQ=255)
                                                                                                                        tAACGATGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGaa                      >  1:887596/1‑149 (MQ=255)
                                                                                                                             aTGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTAc                                                 >  3:843757/1‑117 (MQ=255)
                                                                                                                             aTGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTAc                                                 <  4:843757/117‑1 (MQ=255)
                                                                                                                             aTGAGTAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACat                 >  2:288738/1‑149 (MQ=255)
                                                                                                                                  tAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCGAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTg                        <  2:863099/137‑1 (MQ=255)
                                                                                                                                  tAGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCGAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTg                        >  1:863099/1‑137 (MQ=255)
                                                                                                                                   aGCGTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACatcaatca           <  3:740344/149‑1 (MQ=255)
                                                                                                                                      gTCTGCGATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACATCAATCAGGa        <  3:338132/149‑1 (MQ=255)
                                                                                                                                            gATATAGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACATCAATCAGGACGGAGc  <  3:142488/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTCCCGGTCTTTTCCTGAATGACATGAAGATGATCATCCCAAATATCGCTAAATTTCATACGCGAAATATCACCTAACCGCTGACCAGTAACCAGCGCTAACAGCATGGCATTTCCCATGTAACGATGAGTAGCGTCTGCGATATCGAAGATTTTTTTCCATTCTTCAAGGCTCAGCCGTTGTCGGGTAATTTTTCTTCTTGGTTGTTTAGTGGCTAATGCTGGGTTATAGCCAGGAGGTACTTCTCCGTAGTGCTGCGCCTCTTTGAAAACATCAATCAGGACGGAGC  >  NC_000913/1200052‑1200340

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