Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 277,333 A→G T207A (ACG→GCG)  mmuM → CP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913277,3330AG100.0% 289.8 / NA 75T207A (ACG→GCG) mmuMCP4‑6 prophage; S‑methylmethionine:homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/40);  total (35/40)

GATGCCGGGGCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTACGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGGTGGTGTATCCGAAC  >  NC_000913/277189‑277479
                                                                                                                                                |                                                                                                                                                  
gATGCCGGGGCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGcgccgc                                                                                                                                              <  1:10811/151‑1 (MQ=255)
  tGCCGGGGCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTg                                                                                                                                            >  1:460379/1‑151 (MQ=255)
   gCCGGGGCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTg                                                                                                                                            <  2:267791/150‑1 (MQ=255)
      ggggCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTg                                                                                                                                        <  2:382494/151‑1 (MQ=255)
            gATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTgg                                                                                                                                  >  1:607243/1‑151 (MQ=255)
            gATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTgg                                                                                                                                  >  2:405039/1‑151 (MQ=255)
             aTCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCATGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgt                                                                                                                                 >  2:173537/1‑151 (MQ=255)
                 gctgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgttgc                                                                                                                              <  2:94222/150‑1 (MQ=255)
                 gctgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgttg                                                                                                                               <  1:181027/149‑1 (MQ=255)
                 gctgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgtt                                                                                                                                <  2:844704/148‑1 (MQ=255)
                 gctgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGAATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgtt                                                                                                                                <  2:500756/148‑1 (MQ=255)
                  ctgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgttgcgt                                                                                                                            >  1:658287/1‑151 (MQ=255)
                   tgGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgttgcgtt                                                                                                                           >  1:488545/1‑151 (MQ=255)
                    ggCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGgttgcgtt                                                                                                                           <  2:122624/150‑1 (MQ=255)
                          cGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTgg                                                                                                                      <  1:1078615/149‑1 (MQ=255)
                            aaaCCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCggg                                                                                                                  <  2:1157814/151‑1 (MQ=255)
                              aCCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGtt                                                                                                                >  2:459439/1‑151 (MQ=255)
                              aCCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGtt                                                                                                                <  1:737589/151‑1 (MQ=255)
                                   gCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATcc                                                                                                            >  1:938055/1‑150 (MQ=255)
                                    ccGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGc                                                                                                          <  1:10533/151‑1 (MQ=255)
                                    ccGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGGTATCCGc                                                                                                          >  2:367882/1‑151 (MQ=255)
                                     cGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCa                                                                                                         <  1:1052809/151‑1 (MQ=255)
                                     cGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCa                                                                                                         <  1:502285/151‑1 (MQ=255)
                                        aTTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAgg                                                                                                       >  2:1167638/1‑150 (MQ=255)
                                          ttttCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAggtgg                                                                                                    >  2:1020064/1‑151 (MQ=255)
                                             tCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAgggggggg                                                                                                 >  1:427497/1‑151 (MQ=255)
                                                gagaTTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGgcgc                                                                                              >  2:818839/1‑151 (MQ=255)
                                                  gaTTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGgcg                                                                                               >  1:761204/1‑148 (MQ=255)
                                                   aTTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAg                                                                                           <  2:269283/151‑1 (MQ=255)
                                                   aTTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAg                                                                                           <  1:741549/151‑1 (MQ=255)
                                                   aTTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAg                                                                                           <  1:985133/151‑1 (MQ=255)
                                                    ttGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAgg                                                                                          >  2:1027722/1‑151 (MQ=255)
                                                      gAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGc                                                                                         >  2:999095/1‑150 (MQ=255)
                                                        ggCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCAtt                                                                                      <  2:710129/151‑1 (MQ=255)
                                                        ggCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCAtt                                                                                      <  2:751355/151‑1 (MQ=255)
                                                          cGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTaa                                                                                    <  2:552789/151‑1 (MQ=255)
                                                          cGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTaa                                                                                    <  2:904015/151‑1 (MQ=255)
                                                          cGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGGCGCTAGGCATTa                                                                                     >  2:89348/1‑150 (MQ=255)
                                                              ggCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGGTATCCGCAGGTGGTGGCGCTAGGCCTTAACTGt                                                                                >  2:852631/1‑151 (MQ=255)
                                                               gCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTGCCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTa                                                                               <  1:347370/151‑1 (MQ=255)
                                                               gCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTa                                                                               <  2:1083658/151‑1 (MQ=255)
                                                                 cGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTAtt                                                                             >  2:224540/1‑151 (MQ=255)
                                                                 cGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTAtt                                                                             <  2:736958/151‑1 (MQ=255)
                                                                 cGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGGCGCTAGGCATTAACTGTAtt                                                                             >  1:1027293/1‑151 (MQ=255)
                                                                  gAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTg                                                                            <  2:836915/151‑1 (MQ=255)
                                                                    gCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTgc                                                                           <  1:416664/150‑1 (MQ=255)
                                                                      tGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGGGGCGCTAGGCATTAACTGTATTgcg                                                                          >  2:474224/1‑149 (MQ=255)
                                                                         tGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTgg                                                                      >  2:618801/1‑150 (MQ=255)
                                                                            ccGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGaaaa                                                                  <  2:97308/151‑1 (MQ=255)
                                                                             cGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTCGGTTTTAACTGTATTTCGCTGGaaaa                                                                  >  2:877173/1‑150 (MQ=255)
                                                                                  atCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAAcaccac                                                            <  2:1126736/151‑1 (MQ=255)
                                                                                      gcgTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCt                                                        <  1:172607/151‑1 (MQ=255)
                                                                                       cgTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTg                                                       <  1:679780/151‑1 (MQ=255)
                                                                                       cgTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTg                                                       <  2:1103390/151‑1 (MQ=255)
                                                                                         tGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTg                                                       >  1:1122631/1‑149 (MQ=255)
                                                                                          gcgcGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCg                                                     >  1:636107/1‑150 (MQ=255)
                                                                                            gcgGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCCCCGCTGCGtt                                                   >  1:171166/1‑150 (MQ=255)
                                                                                              gggCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGCGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGtt                                                   >  1:1087063/1‑148 (MQ=255)
                                                                                                gCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTgcagc                                              <  2:1183756/151‑1 (MQ=255)
                                                                                                   tGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCAtt                                           >  2:663263/1‑151 (MQ=255)
                                                                                                    ggTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCAtt                                           <  1:290221/150‑1 (MQ=255)
                                                                                                        ctcATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTacac                                      >  2:796310/1‑151 (MQ=255)
                                                                                                               aCCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTa                                <  1:476928/150‑1 (MQ=255)
                                                                                                               aCCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGGTTa                                >  2:794622/1‑150 (MQ=255)
                                                                                                               aCCCTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGGGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGttt                                 >  2:1054224/1‑149 (MQ=255)
                                                                                                                 ccTGCGCGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACg                             >  2:139811/1‑151 (MQ=255)
                                                                                                                    gcgcGACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACgg                            <  1:1084155/149‑1 (MQ=255)
                                                                                                                     cgcgACAGCGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGtgc                         >  1:9521/1‑151 (MQ=255)
                                                                                                                             cGAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCtg                 <  2:708702/151‑1 (MQ=255)
                                                                                                                              gAACACCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGGGCTGCCGCtgg                >  1:487135/1‑151 (MQ=255)
                                                                                                                                acacCTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCtggtg              <  2:339442/151‑1 (MQ=255)
                                                                                                                                   ccTGAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGgtggtg           <  2:354315/151‑1 (MQ=255)
                                                                                                                                      gAGCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGGTGGTGTa         <  1:474110/150‑1 (MQ=255)
                                                                                                                                        gCGACGGTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGGTGGTGTATcc      <  2:761377/151‑1 (MQ=255)
                                                                                                                                            aggTGCGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGGTGGTGTATCCGAAc  >  2:531495/2‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GATGCCGGGGCCGATCTGCTGGCCTGCGAAACCCTGCCGAATTTTTCCGAGATTGAGGCGTTGGCCGAGCTGTTGACCGCATATCCGCGTGCGCGGGCGTGGTTCTCATTTACCCTGCGCGACAGCGAACACCTGAGCGACGGTACGCCGCTGCGTGACGTGGTTGCGTTGCTGGCGGGTTATCCGCAGGTGGTGGCGCTAGGCATTAACTGTATTGCGCTGGAAAACACCACCGCTGCGTTGCAGCATTTACACGGTTTAACGGTGCTGCCGCTGGTGGTGTATCCGAAC  >  NC_000913/277189‑277479

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