Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,536,736 A→G Y287H (TAC→CAC)  nanS ← N‑acetyl‑9‑O‑acetylneuraminate esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,7360AG100.0% 315.0 / NA 90Y287H (TAC→CAC) nanSN‑acetyl‑9‑O‑acetylneuraminate esterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/53);  total (37/53)

GAATAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTACGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAAT  >  NC_000913/4536599‑4536871
                                                                                                                                         |                                                                                                                                       
gAATAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCg                                                                                                                                       <  1:1343181/140‑1 (MQ=255)
 aaTAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTgcgc                                                                                                                                      <  2:1329507/140‑1 (MQ=255)
   tAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTg                                                                                                                                    >  2:507268/1‑140 (MQ=255)
   tAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTg                                                                                                                                    <  1:720288/140‑1 (MQ=255)
   tAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTg                                                                                                                                    <  1:15860/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  1:1014544/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  2:396578/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  2:911670/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  1:945798/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  1:888898/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  1:512384/140‑1 (MQ=255)
     aaCTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  1:315520/140‑1 (MQ=255)
          cccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTa                                                                                                                             >  2:907600/1‑140 (MQ=255)
          cccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTa                                                                                                                             >  2:503037/1‑140 (MQ=255)
          cccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTa                                                                                                                             >  1:445606/1‑140 (MQ=255)
           ccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTaa                                                                                                                            >  2:746643/1‑140 (MQ=255)
           ccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTaa                                                                                                                            >  2:1184475/1‑140 (MQ=255)
           ccGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTaa                                                                                                                            >  2:894993/1‑140 (MQ=255)
            cGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  <  2:919308/133‑1 (MQ=255)
            cGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGaa                                                                                                                                  >  1:919308/1‑133 (MQ=255)
            cGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAat                                                                                                                           >  2:825049/1‑140 (MQ=255)
            cGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAat                                                                                                                           >  1:633731/1‑140 (MQ=255)
              cTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTg                                                                                                                                    >  2:1109331/1‑129 (MQ=255)
              cTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTg                                                                                                                                    <  1:1109331/129‑1 (MQ=255)
              cTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAatat                                                                                                                         >  1:126235/1‑140 (MQ=255)
                 aCCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATccc                                                                                                                      <  2:295885/140‑1 (MQ=255)
                 aCCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATccc                                                                                                                      >  2:164803/1‑140 (MQ=255)
                      tCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCt                                                                                                                 >  1:190713/1‑140 (MQ=255)
                             aaaaCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCg                                                                                                          >  2:356380/1‑140 (MQ=255)
                             aaaaCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCg                                                                                                          >  2:13147/1‑140 (MQ=255)
                             aaaaCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCg                                                                                                          >  2:509888/1‑140 (MQ=255)
                             aaaaCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCg                                                                                                          >  2:662384/1‑140 (MQ=255)
                             aaaaCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCg                                                                                                          >  1:119513/1‑140 (MQ=255)
                                aCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTcc                                                                                                       <  1:1378843/140‑1 (MQ=255)
                                aCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTcc                                                                                                       <  1:164803/140‑1 (MQ=255)
                                aCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTcc                                                                                                       <  2:1319130/140‑1 (MQ=255)
                                aCTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTcc                                                                                                       <  2:1208464/140‑1 (MQ=255)
                                   gCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCggg                                                                                                    <  1:828363/140‑1 (MQ=255)
                                             ttCTACTAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAg                                                                                          <  1:1079489/140‑1 (MQ=255)
                                               cTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGgc                                                                                        <  1:798266/140‑1 (MQ=255)
                                               cTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGgc                                                                                        <  2:428609/140‑1 (MQ=255)
                                               cTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGgc                                                                                        <  2:257735/140‑1 (MQ=255)
                                                tACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGgcg                                                                                       <  1:950769/140‑1 (MQ=255)
                                                  cAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGcgcg                                                                                     <  1:503037/140‑1 (MQ=255)
                                                  cAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGcgcg                                                                                     >  2:898722/1‑140 (MQ=255)
                                                   aaaTCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGt                                                                                    <  1:825049/140‑1 (MQ=255)
                                                      tCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCg                                                                                 >  2:535300/1‑140 (MQ=255)
                                                      tCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCg                                                                                 >  2:461014/1‑140 (MQ=255)
                                                       cTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGt                                                                                <  1:49749/140‑1 (MQ=255)
                                                       cTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGt                                                                                <  2:983989/140‑1 (MQ=255)
                                                           cAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGt                                                                            <  2:190713/140‑1 (MQ=255)
                                                            aGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTc                                                                           <  2:348264/140‑1 (MQ=255)
                                                            aGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTc                                                                           <  2:1176999/140‑1 (MQ=255)
                                                                 taatCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCtt                                                                      <  1:535300/140‑1 (MQ=255)
                                                                 taatCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCtt                                                                      <  1:461014/140‑1 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGt                                                                                >  1:1223575/1‑127 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGt                                                                                <  2:1223575/127‑1 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCa                                                                   >  2:790243/1‑140 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCa                                                                   >  2:1040236/1‑140 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCa                                                                   >  2:319322/1‑140 (MQ=255)
                                                                    tCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCa                                                                   >  2:239031/1‑140 (MQ=255)
                                                                         cGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCttg                                                              <  1:13147/140‑1 (MQ=255)
                                                                              ggCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCt                                                         >  2:578698/1‑140 (MQ=255)
                                                                                        cTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTcgac                                               <  2:279197/140‑1 (MQ=255)
                                                                                        cTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTcgac                                               <  1:356380/140‑1 (MQ=255)
                                                                                        cTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTcgac                                               <  1:578698/140‑1 (MQ=255)
                                                                                        cTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTcgac                                               >  2:377672/1‑140 (MQ=255)
                                                                                        cTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTcgac                                               <  2:986708/140‑1 (MQ=255)
                                                                                                actacTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAaataa                                       >  2:500648/5‑140 (MQ=255)
                                                                                                aCTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTTAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAaataa                                       >  1:495116/1‑140 (MQ=255)
                                                                                                aCTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAaataat                                      >  2:103227/1‑140 (MQ=255)
                                                                                                  tGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAataata                                     >  2:960597/1‑140 (MQ=255)
                                                                                                  tGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAataata                                     <  1:907600/140‑1 (MQ=255)
                                                                                                         cAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTa                              <  1:1398534/140‑1 (MQ=255)
                                                                                                         cAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTa                              <  2:633731/140‑1 (MQ=255)
                                                                                                              ccGTCGTCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACa                         <  1:618687/140‑1 (MQ=255)
                                                                                                               cgtcgtCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACat                        <  1:1184475/140‑1 (MQ=255)
                                                                                                               cgtcgtCCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACat                        <  1:1290663/140‑1 (MQ=255)
                                                                                                                 tcgtcCAGTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGCAGTCGACGAAAATAATATTGGCTAAAACatta                      >  1:165361/1‑140 (MQ=255)
                                                                                                                        gTTCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGAt               >  2:242158/1‑140 (MQ=255)
                                                                                                                         ttCTCCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  1:103227/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  2:628515/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  2:445606/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  2:43433/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  2:1407757/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  1:380164/140‑1 (MQ=255)
                                                                                                                          tctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATaa             <  1:1099464/140‑1 (MQ=255)
                                                                                                                           ctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAAt            >  2:114129/1‑140 (MQ=255)
                                                                                                                           ctcCGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAAt            >  2:1097228/1‑140 (MQ=255)
                                                                                                                              cGGTGACCGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGc         >  1:540435/1‑140 (MQ=255)
                                                                                                                                    ccGGTGCGCTGAACCGTAATACCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATaaa   <  2:1284460/140‑1 (MQ=255)
                                                                                                                                     cGGTGCGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAAt  <  1:377672/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
GAATAAACTCCCCGCTCACCTTTCGCGCCAAAACTGCAAAATTGCTTCTACAAATCTGTCAGAAATAATCCCCCGACGGGCTGCAGTGCTGAAATGACTGCTTCGCAGTGCCGTCGTCCAGTTCTCCGGTGACCGGTACGCTGAACCGTAATATCCCGTGCTTAAATCGTCCGGGTCTTCATCAGGCGCGTTCGTCAGTCCTCTTTCACCTTGTTGCTGGAAGTCGACGAAAATAATATTGGCTAAAACATTATTTTGATAATTGCCATAAAT  >  NC_000913/4536599‑4536871

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