Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 397,273 C→T 5.9% A212V (GCC→GTC)  sbmA → peptide antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,2730CT 5.9% 578.9 / 18.0 170A212V (GCC→GTC) sbmApeptide antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base C (81/79);  new base T (6/4);  total (87/83)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 176 total aligned reads displayed.

GTGTTCCGCTGGCGTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTA  >  NC_000913/397125‑397421
                                                                                                                                                    |                                                                                                                                                    
 tgtTCCGCTGGCGTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGcc                                                                                                                                                     >  1:815620/1‑149 (MQ=255)
   ttCCGCTGGCGTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcat                                                                                                                                                   >  1:753704/1‑149 (MQ=255)
   ttCCGCTGGCGTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcat                                                                                                                                                   >  3:273025/1‑149 (MQ=255)
                cAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATc                                                                                                                                      <  1:1057601/149‑1 (MQ=255)
                  gCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGc                                                                                                                                    <  4:273025/149‑1 (MQ=255)
                   cGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGcc                                                                                                                                   <  3:1297389/149‑1 (MQ=255)
                      tGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTc                                                                                                                                <  3:527001/149‑1 (MQ=255)
                         aCGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTg                                                                                                                             >  3:643260/1‑149 (MQ=255)
                           gAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGcc                                                                                                                           >  1:153310/1‑149 (MQ=255)
                           gAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGcc                                                                                                                           <  2:815620/149‑1 (MQ=255)
                           gAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGcc                                                                                                                           <  1:1463506/149‑1 (MQ=255)
                             atatTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCgg                                                                                                                         >  1:56507/1‑149 (MQ=255)
                              tatTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGt                                                                                                                        <  2:56507/149‑1 (MQ=255)
                                    aTGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGt                                                                                                                  <  2:153310/149‑1 (MQ=255)
                                    aTGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGt                                                                                                                  <  2:753704/149‑1 (MQ=255)
                                     tGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTa                                                                                                                 >  3:799237/1‑149 (MQ=255)
                                      ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTaa                                                                                                                <  1:382425/149‑1 (MQ=255)
                                      ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTaa                                                                                                                <  4:799237/149‑1 (MQ=255)
                                       gCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAAc                                                                                                               >  2:693916/1‑149 (MQ=255)
                                             tGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGa                                                                                                                                              <  3:291723/112‑1 (MQ=255)
                                             tGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGa                                                                                                                                              >  4:291723/1‑112 (MQ=255)
                                                 aacaacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATc                                                                                                                                      >  2:765252/1‑116 (MQ=255)
                                                    aacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTaa                                                                                                                >  1:683818/1‑135 (MQ=255)
                                                    aacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGc                                                                                                             >  4:1202582/1‑138 (MQ=255)
                                                    aacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGc                                                                                                             <  3:1202582/138‑1 (MQ=255)
                                                           tCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCgg                                                                                           >  3:1150368/1‑149 (MQ=255)
                                                               atCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCt                                                                                       >  2:447201/1‑149 (MQ=255)
                                                                tCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGc                                                                                                                            >  2:777348/1‑111 (MQ=255)
                                                                 cGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGc                                                                                     <  4:1150368/149‑1 (MQ=255)
                                                                  gAAGGGGCCGCTCAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGcc                                                                                    >  1:767734/1‑149 (MQ=255)
                                                                      gggCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGAtt                                                                                >  4:1392719/1‑149 (MQ=255)
                                                                        gCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcatc                                                                                                                                                  <  2:192388/81‑1 (MQ=255)
                                                                        gCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcatc                                                                                                                                                  <  2:402988/81‑1 (MQ=255)
                                                                        gCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcatc                                                                                                                                                  >  1:402988/1‑81 (MQ=255)
                                                                        gCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCcatc                                                                                                                                                  >  1:192388/1‑81 (MQ=255)
                                                                           gCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCgg                                                                           <  1:763708/149‑1 (MQ=255)
                                                                           gCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCgg                                                                           >  4:326730/1‑149 (MQ=255)
                                                                           gCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCgg                                                                           <  1:447201/149‑1 (MQ=255)
                                                                            cacaGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGc                                                                                                                            >  3:874672/1‑99 (MQ=255)
                                                                             acaGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGc                                                                         >  3:500327/1‑149 (MQ=255)
                                                                                 cGTGTGCAGGAAGACACCATGCGTTTTGCTACAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACAt                                                                     <  2:32484/149‑1 (MQ=255)
                                                                                      gCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGa                                                                                  >  4:303976/1‑131 (MQ=255)
                                                                                        aGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTAt                                                              >  3:755980/1‑149 (MQ=255)
                                                                                             gACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCt                                                         <  4:755980/149‑1 (MQ=255)
                                                                                              acacCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTg                                                        <  4:664163/149‑1 (MQ=255)
                                                                                              acacCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTg                                                        <  1:693916/149‑1 (MQ=255)
                                                                                              acacCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTg                                                        <  3:1392719/149‑1 (MQ=255)
                                                                                              acacCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTg                                                        <  4:82689/149‑1 (MQ=255)
                                                                                                acCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTcc                                                                                                        >  1:588710/1‑99 (MQ=255)
                                                                                                acCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCt                                                                                       <  4:132401/116‑1 (MQ=255)
                                                                                                acCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGt                                                      <  3:326730/149‑1 (MQ=255)
                                                                                                acCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGt                                                      <  2:635903/149‑1 (MQ=255)
                                                                                                 ccATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTg                                                     >  4:201914/1‑149 (MQ=255)
                                                                                                 ccATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCGCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTg                                                     >  4:1217843/1‑149 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGc                                                                                              >  1:646862/1‑106 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGc                                                                                              <  2:646862/106‑1 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAg                                                                                         <  3:1377125/111‑1 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGt                                                           >  2:783082/1‑141 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGt                                                           <  1:783082/141‑1 (MQ=255)
                                                                                                   aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                   <  1:989671/149‑1 (MQ=255)
                                                                                                    tGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAtt                                                  <  3:1150096/149‑1 (MQ=255)
                                                                                                       gTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGcc                                               >  3:831554/1‑149 (MQ=255)
                                                                                                        ttttGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCg                                              >  1:1105014/1‑149 (MQ=255)
                                                                                                            gCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAAt                                          >  4:477449/1‑149 (MQ=255)
                                                                                                             cTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATc                                         <  4:712698/149‑1 (MQ=255)
                                                                                                              ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                               <  4:1408560/126‑1 (MQ=255)
                                                                                                              ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                               >  3:1408560/1‑126 (MQ=255)
                                                                                                              ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATc                                         >  1:1234206/1‑148 (MQ=255)
                                                                                                              ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCg                                        >  1:1277490/1‑149 (MQ=255)
                                                                                                                        gagaATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGc                                                                                     <  3:942358/94‑1 (MQ=255)
                                                                                                                        gagaATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTgg                                                       <  3:593872/124‑1 (MQ=255)
                                                                                                                        gagaATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTgg                                                       >  4:593872/1‑124 (MQ=255)
                                                                                                                            atatGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATg                          <  2:74656/149‑1 (MQ=255)
                                                                                                                                gggCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGa                            <  3:1458879/143‑1 (MQ=255)
                                                                                                                                    gTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATg                          <  4:1408731/141‑1 (MQ=255)
                                                                                                                                       aGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGt                                                           >  2:989859/1‑105 (MQ=255)
                                                                                                                                         ttttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGt                                                      >  1:760028/1‑108 (MQ=255)
                                                                                                                                         ttttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGt                                                      <  2:760028/108‑1 (MQ=255)
                                                                                                                                         ttttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGc             >  1:30269/1‑149 (MQ=255)
                                                                                                                                         ttttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGc             <  3:201914/149‑1 (MQ=255)
                                                                                                                                         ttttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGc             >  1:532636/1‑149 (MQ=255)
                                                                                                                                          tttATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCt            <  4:102737/149‑1 (MQ=255)
                                                                                                                                          ccgattaaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCt            <  1:933940/143‑1 (MQ=255)
                                                                                                                                            tATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTgg          <  3:305340/149‑1 (MQ=255)
                                                                                                                                            gaTCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCaa                                           >  2:1267446/2‑116 (MQ=255)
                                                                                                                                            gaTCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCt                                     <  1:745419/121‑1 (MQ=255)
                                                                                                                                            gaTCAACGTCATAATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTgg          >  3:1993/2‑149 (MQ=255)
                                                                                                                                              tCAACGCCATCATGACGTTGATCGTCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATgggg                       <  4:914036/134‑1 (MQ=255)
                                                                                                                                               cAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCaa                                           <  2:392689/113‑1 (MQ=255)
                                                                                                                                                aaCGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGa                                                    <  2:1436221/103‑1 (MQ=255)
                                                                                                                                                aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                               <  3:678537/92‑1 (MQ=255)
                                                                                                                                                aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTg                                    >  4:332343/1‑119 (MQ=255)
                                                                                                                                                aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTg                                    <  3:332343/119‑1 (MQ=255)
                                                                                                                                                 aCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGTACATTCCGTAt                                                              <  4:380031/92‑1 (MQ=255)
                                                                                                                                                 aCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGTACATTCCGTAt                                                              >  3:380031/1‑92 (MQ=255)
                                                                                                                                                   gCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGc                                                <  2:1110678/104‑1 (MQ=255)
                                                                                                                                                   gCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTc                                      >  1:357994/1‑114 (MQ=255)
                                                                                                                                                   gCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGt   <  4:1092137/149‑1 (MQ=255)
                                                                                                                                                    ccATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCAtgt                                                                                                >  3:1078919/1‑55 (MQ=255)
                                                                                                                                                    ccATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTa  <  2:731138/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGTTCCGCTGGCGTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTA  >  NC_000913/397125‑397421

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