Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 219,624 A→G C124C (TGT→TGC rcsF ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913219,6240AG100.0% 245.0 / NA 66C124C (TGT→TGCrcsFputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/39);  total (27/39)

ATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGC  >  NC_000913/219479‑219770
                                                                                                                                                 |                                                                                                                                                  
ataaagaaaaGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGatgcata                                                                                                                                                >  1:335235/10‑150 (MQ=255)
aTGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGatgcatac                                                                                                                                               <  1:331907/151‑1 (MQ=255)
     aTGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGcc                                                                                                                                           <  1:513471/150‑1 (MQ=255)
           ttCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGAc                                                                                                                                       <  2:180839/148‑1 (MQ=255)
               ccGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATagg                                                                                                                                 <  1:22296/150‑2 (MQ=255)
               ccGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATagc                                                                                                                                 <  1:1017307/150‑1 (MQ=255)
               ccGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATagc                                                                                                                                 <  2:536734/150‑1 (MQ=255)
                  ttGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTGCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGcc                                                                                                                             <  1:1019375/151‑1 (MQ=255)
                      cgcTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGc                                                                                                                         <  2:955738/151‑1 (MQ=255)
                          ttttACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTa                                                                                                                      >  1:85788/1‑150 (MQ=255)
                          ttttACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTAc                                                                                                                     <  1:2661/151‑1 (MQ=255)
                            ttACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCg                                                                                                                   <  1:671562/151‑1 (MQ=255)
                             tACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGc                                                                                                                  <  1:451635/151‑1 (MQ=255)
                                  gACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTg                                                                                                             >  2:641393/1‑151 (MQ=255)
                                             cgcgGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGc                                                                                                   <  2:819767/150‑1 (MQ=255)
                                                   aCGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAg                                                                                            <  2:840944/151‑1 (MQ=255)
                                                       cGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTaa                                                                                         <  2:96297/150‑1 (MQ=255)
                                                       cGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTaa                                                                                         <  1:717155/150‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCNNNCGAATAACNNNNNNTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  1:1392/151‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  2:48288/151‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  2:1050074/151‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  2:174411/151‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  2:543731/151‑1 (MQ=255)
                                                         aaCTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATAc                                                                                      <  1:795391/151‑1 (MQ=255)
                                                           cTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAg                                                                                    <  1:423574/151‑1 (MQ=255)
                                                           cTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAg                                                                                    <  1:237902/151‑1 (MQ=255)
                                                           cTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAg                                                                                    <  2:566168/151‑1 (MQ=255)
                                                           cTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAg                                                                                    <  1:131698/151‑1 (MQ=255)
                                                             tttCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCa                                                                                  >  2:597090/1‑151 (MQ=255)
                                                               tCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCAtt                                                                                >  2:899034/1‑151 (MQ=255)
                                                               tCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCAtt                                                                                >  1:989860/1‑151 (MQ=255)
                                                                cTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACCGCATTg                                                                               >  2:494062/1‑151 (MQ=255)
                                                                cTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTg                                                                               >  1:576445/1‑151 (MQ=255)
                                                                cTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTg                                                                               >  1:359566/1‑151 (MQ=255)
                                                                 tttATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTgg                                                                              >  1:350160/1‑151 (MQ=255)
                                                                  ttATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGc                                                                             <  1:397571/151‑1 (MQ=255)
                                                                            agcgaATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTa                                                                   >  2:451848/1‑151 (MQ=255)
                                                                                aaTAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAgg                                                               <  2:808563/151‑1 (MQ=255)
                                                                                     cGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGTCGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTg                                                          >  1:900136/1‑151 (MQ=255)
                                                                                        cTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGAtt                                                       >  1:541030/1‑151 (MQ=255)
                                                                                                   aaCTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCtt                                            >  1:124975/1‑151 (MQ=255)
                                                                                                       ggAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGt                                        >  1:962608/1‑151 (MQ=255)
                                                                                                         aaaCTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTg                                       >  2:265985/1‑150 (MQ=255)
                                                                                                            cTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGt                                   <  1:382433/151‑1 (MQ=255)
                                                                                                              gCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGATTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTg                                 >  2:982364/1‑151 (MQ=255)
                                                                                                                   tttCGCCGTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATg                            <  2:845941/151‑1 (MQ=255)
                                                                                                                          gTAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTcggcg                     >  2:890263/1‑151 (MQ=255)
                                                                                                                           tAATGTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTcggcgg                    <  1:620086/151‑1 (MQ=255)
                                                                                                                               gTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGagag                >  2:38066/1‑151 (MQ=255)
                                                                                                                               gTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGagag                >  1:691258/1‑151 (MQ=255)
                                                                                                                               gTTAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGATTTCATTTTAGCGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGagag                >  2:664942/1‑151 (MQ=255)
                                                                                                                                ttAAGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCGGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGt               <  2:264913/151‑1 (MQ=255)
                                                                                                                                  aaGCGCAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTcc             <  1:685897/151‑1 (MQ=255)
                                                                                                                                    gcgcAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTc              >  1:876257/1‑148 (MQ=255)
                                                                                                                                    gcgcAGAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTc              >  1:688384/1‑148 (MQ=255)
                                                                                                                                         gAACCGATGCATACAGCCTGACGATAGCCGCCTGGGGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCCTCCGCTTACGTGCGGTTGGAATGCTCGGCGGGGAGTCCTGATTag      >  2:317285/1‑151 (MQ=255)
                                                                                                                                         gAACCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTag      >  1:533490/1‑151 (MQ=255)
                                                                                                                                          aaCCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt        <  2:688640/148‑1 (MQ=255)
                                                                                                                                          aaCCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt        <  2:85797/148‑1 (MQ=255)
                                                                                                                                          aaCCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt        <  2:977997/148‑1 (MQ=255)
                                                                                                                                          aaCCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt        <  1:983727/148‑1 (MQ=255)
                                                                                                                                          aaCCGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTag      <  1:790729/150‑1 (MQ=255)
                                                                                                                                            ccGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAgg   <  2:824501/151‑1 (MQ=255)
                                                                                                                                            ccGATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAgg   <  2:357959/151‑1 (MQ=255)
                                                                                                                                              gATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc                           >  1:63854/1‑125 (MQ=255)
                                                                                                                                              gATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGc                           <  2:63862/125‑1 (MQ=255)
                                                                                                                                              gATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGc  >  1:598321/1‑150 (MQ=255)
                                                                                                                                              gATGCATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGc  >  1:140851/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATGAGATGCAGTTCGCCGTTGGCGCTTTTTACCAGACCTGGCTGGCGCGGAACGGCGAACTTTTCTTTATAGGGCGAGCGAATAACGCCTATTTGCTCGAACTGGAAACTGCTCATTTCGCCGTAATGTTAAGCGCAGAACCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGC  >  NC_000913/219479‑219770

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