Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 600,631 A→G M640V (ATG→GTG)  cusA → copper/silver efflux system, membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913600,6310AG100.0% 250.9 / NA 69M640V (ATG→GTG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/32);  total (37/32)

AACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGG  >  NC_000913/600486‑600773
                                                                                                                                                 |                                                                                                                                              
aaCCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACg                                                                                                                                           >  1:571463/1‑151 (MQ=255)
aaCCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACg                                                                                                                                           <  1:532244/151‑1 (MQ=255)
  ccGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGAt                                                                                                                                         <  1:184063/151‑1 (MQ=255)
   cGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATg                                                                                                                                        >  1:754298/1‑151 (MQ=255)
     aCAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATgg                                                                                                                                       >  2:298781/1‑150 (MQ=255)
     aCAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATgg                                                                                                                                       >  1:607477/1‑150 (MQ=255)
       aaGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACa                                                                                                                                    <  2:918476/151‑1 (MQ=255)
           tcaTTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCCAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGATGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGTCAGGCGTGACGATGGACaaaa                                                                                                                                 >  2:279643/3‑150 (MQ=255)
            aaTTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATc                                                                                                                               <  1:339673/151‑1 (MQ=255)
              ttATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCAt                                                                                                                             >  2:9650/1‑151 (MQ=255)
               tATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCAt                                                                                                                             <  2:20635/150‑1 (MQ=255)
                aTGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTg                                                                                                                           >  2:454389/1‑151 (MQ=255)
                   aGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAgg                                                                                                                        >  2:286180/1‑151 (MQ=255)
                     cGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGaa                                                                                                                      <  2:573857/151‑1 (MQ=255)
                       tACCTGAAGTGTCGCGGGTATTTGGCAAAACCGGTAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGACGCAGGAGCAGTGGCGGCCAGGCGTGACGCTGGACAAAATCATTGAGGAAc                                                                                                                     <  2:519024/150‑1 (MQ=255)
                           tGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGAt                                                                                                                >  2:340051/1‑151 (MQ=255)
                            gAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATAGAGGAACTGGAt                                                                                                                <  1:692483/150‑1 (MQ=255)
                              aGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATaa                                                                                                              >  2:705706/1‑150 (MQ=255)
                                   cgcgGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGt                                                                                                        >  1:867900/1‑151 (MQ=255)
                                    gcgGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTg                                                                                                       <  1:723350/151‑1 (MQ=255)
                                      gggTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCg                                                                                                     >  1:300537/1‑151 (MQ=255)
                                        gTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGc                                                                                                   <  1:322226/151‑1 (MQ=255)
                                          aTTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCt                                                                                                  <  2:229537/150‑1 (MQ=255)
                                           tttGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGc                                                                                                >  1:752592/1‑151 (MQ=255)
                                              ggCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCg                                                                                              <  2:720233/150‑1 (MQ=255)
                                                 aaaaCCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCgggg                                                                                           >  2:377027/1‑150 (MQ=255)
                                                       gggAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGa                                                                                    <  2:803326/151‑1 (MQ=255)
                                                          aaaGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATc                                                                                 >  1:599276/1‑151 (MQ=255)
                                                           aaGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCt                                                                                >  2:731765/1‑151 (MQ=255)
                                                              cGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGTTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCtgtg                                                                             >  1:462477/1‑151 (MQ=255)
                                                                gAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTggg                                                                           <  1:801369/151‑1 (MQ=255)
                                                                          ccGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAAt                                                                 <  2:231997/151‑1 (MQ=255)
                                                                           cGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAAtt                                                                <  2:412245/151‑1 (MQ=255)
                                                                            gATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTc                                                               <  2:191509/151‑1 (MQ=255)
                                                                            gATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTc                                                               >  2:897417/1‑151 (MQ=255)
                                                                            gATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTc                                                               >  1:798147/1‑151 (MQ=255)
                                                                             aTTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCg                                                              <  1:1016632/151‑1 (MQ=255)
                                                                                 tGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTaa                                                           >  2:254154/1‑150 (MQ=255)
                                                                                 tGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAAc                                                          >  2:1046866/1‑151 (MQ=255)
                                                                                  gCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTa                                                            <  2:1006668/148‑1 (MQ=255)
                                                                                   cTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCg                                                        >  2:655050/1‑151 (MQ=255)
                                                                                    tCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGt                                                       >  1:869370/1‑151 (MQ=255)
                                                                                    tCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGt                                                       >  1:244321/1‑151 (MQ=255)
                                                                                     ccGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGGCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGACGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGt                                                       >  2:663170/1‑150 (MQ=255)
                                                                                     ccGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTa                                                      <  2:198355/151‑1 (MQ=255)
                                                                                        cTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCg                                                   >  1:864937/1‑151 (MQ=255)
                                                                                          ggAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGa                                                  <  1:821146/150‑1 (MQ=255)
                                                                                           gagaTGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACTGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGat                                                 >  2:976485/1‑150 (MQ=255)
                                                                                             gaTGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATg                                              <  2:708175/151‑1 (MQ=255)
                                                                                                      aaCGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAAcc                                     >  1:527222/1‑151 (MQ=255)
                                                                                                              tCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTa                              >  1:339187/1‑150 (MQ=255)
                                                                                                               ccAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTaa                             >  1:956630/1‑150 (MQ=255)
                                                                                                               ccAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTaa                             <  1:244466/150‑1 (MQ=255)
                                                                                                                 aGCTTAAGCCGCAGGAGCAGTGGCGGCCCGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTaaa                            >  2:966064/1‑149 (MQ=255)
                                                                                                                    ttAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGAGCCGCCAATTCGTAACCGTGTCGATATGCTCTCAACCGGCATTAAAAGccc                       >  2:34066/1‑151 (MQ=255)
                                                                                                                     tAAGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGcc                        <  2:245157/149‑1 (MQ=255)
                                                                                                                       aGCCGCAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCAt                    <  2:412841/151‑1 (MQ=255)
                                                                                                                            cAGGAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCa               <  1:542479/151‑1 (MQ=255)
                                                                                                                              ggAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCAtt             >  1:943426/1‑151 (MQ=255)
                                                                                                                              ggAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCAtt             >  1:22868/1‑151 (MQ=255)
                                                                                                                               gAGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTa            <  1:168718/151‑1 (MQ=255)
                                                                                                                                aGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTaa           <  1:34067/151‑1 (MQ=255)
                                                                                                                                aGCAGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTa            >  2:348962/1‑150 (MQ=255)
                                                                                                                                   aGTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGt        >  1:153355/1‑151 (MQ=255)
                                                                                                                                    gTGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGtt       <  1:464672/151‑1 (MQ=255)
                                                                                                                                     tGGCGGCCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGttt      >  2:613685/1‑151 (MQ=255)
                                                                                                                                       gcggcCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGtt       <  2:134585/148‑1 (MQ=255)
                                                                                                                                        cggcCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCg   <  2:556826/151‑1 (MQ=255)
                                                                                                                                         ggcCAGGCGTGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCgg  <  1:562616/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAAACCGGGAAAGCGGAAACCGCCACCGATTCTGCTCCGCTGGAGATGGTAGAAACGACCATCCAGCTTAAGCCGCAGGAGCAGTGGCGGCCAGGCATGACGATGGACAAAATCATTGAGGAACTGGATAACACCGTGCGGCTGCCGGGGCTGGCGAATCTGTGGGTGCCGCCAATTCGTAACCGTATCGATATGCTCTCAACCGGCATTAAAAGCCCCATCGGCATTAAAGTTTCCGG  >  NC_000913/600486‑600773

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