Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,179,383 G→A K423K (AAG→AAA tolC → outer membrane channel TolC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,179,3830GA96.0% 154.1 / NA 50K423K (AAG→AAAtolCouter membrane channel TolC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/30);  minor base C (2/0);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

CATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGA  >  NC_000913/3179138‑3179626
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTTTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGc                                                                                                                                                                                                                                                <  1:343047/251‑1 (MQ=255)
     tAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTgg                                                                                                                                                                                                                                           >  2:115003/1‑251 (MQ=255)
       gCGTCGTGCAGACCGTGCGTTCCTCCTCCACCAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGt                                                                                                                                                                                                                                         <  1:254367/251‑1 (MQ=255)
              gCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTg                                                                                                                                                                                                                                   <  1:141083/250‑1 (MQ=255)
                  aCCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACg                                                                                                                                                                                                                               >  1:254000/1‑250 (MQ=255)
                   ccGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAg                                                                                                                                                                                                                             <  2:71691/251‑1 (MQ=255)
                    cGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATGAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAAGCAGCTCTGGGTACGTTGAACGAGc                                                                                                                                                                                                                            <  1:28817/251‑1 (MQ=255)
                              tcctTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATctgct                                                                                                                                                                                                                  >  1:209394/1‑251 (MQ=255)
                                        acaTTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAAc                                                                                                                                                                                                        <  1:366963/251‑1 (MQ=255)
                                               tGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGATTGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATgcgc                                                                                                                                                                                                 >  2:381542/1‑251 (MQ=255)
                                                 cATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCCTTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTg                                                                                                                                                                                               >  2:85187/1‑251 (MQ=255)
                                                 cATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCGGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATTCGCTg                                                                                                                                                                                               <  2:389409/251‑1 (MQ=255)
                                                      aTCAGTAGCAGTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGGACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGc                                                                                                                                                                                            <  1:80522/249‑1 (MQ=255)
                                                               aTTAACGCCTACAACCAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTc                                                                                                                                                                                 >  2:275879/1‑251 (MQ=255)
                                                                   aCGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACt                                                                                                                                                                             <  1:204357/251‑1 (MQ=255)
                                                                             acaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGaa                                                                                                                                                                    >  1:336692/1‑250 (MQ=255)
                                                                                         ttCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCg                                                                                                                                                        <  1:93265/250‑1 (MQ=255)
                                                                                             gCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGCGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                                                                                                                   >  2:131560/1‑251 (MQ=255)
                                                                                                       cATTAGACGCGGTGGAAGCGGGCTGCTCGGTCGGGACGCGTACCATTGGTGATGTGTTGGATGCGACCTCCACGTTGTACAACGCCAAGCAAGAGCTGTAGAATGCGCGTTATAACTACCTGAGTAAGCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTTAGCAATACGGTGTCTACTAATCCGGAAAACGTTGCACCGCCACCGCCGTAAc                                                                                                                                           <  2:372924/249‑1 (MQ=255)
                                                                                                       cATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACa                                                                                                                                          <  1:321165/250‑1 (MQ=255)
                                                                                                       cATTAGACGCGATGGAAGCGGGCTACGCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAg                                                                                                                                         <  1:275879/251‑1 (MQ=255)
                                                                                                           aGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAAt                                                                                                                                      <  1:150107/250‑1 (MQ=255)
                                                                                                                  aTGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGc                                                                                                                              <  1:115003/251‑1 (MQ=255)
                                                                                                                     gAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGa                                                                                                                           <  2:171501/251‑1 (MQ=255)
                                                                                                                           ggCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTa                                                                                                                     <  1:261998/251‑1 (MQ=255)
                                                                                                                                   ggtcggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTg                                                                                                             <  1:239357/251‑1 (MQ=255)
                                                                                                                                      cggtACGCGTACCATTGTTGATGTGTTGGCTGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGATTCCCCTACTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGAt                                                                                                           >  2:384170/1‑250 (MQ=255)
                                                                                                                                                   attgttactgtgTTGGATGCGACCACCACGTTGTAAACCGCCTAGCAAGAGCTGGCGGATGCGCGTTATAACTACCTGATTAATCAGCTGAATCTTAACTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCCAACGCCGGAACCGAATGCTATTGCTGATAGTTATGCGCCTGATAGCCCGGCCCCAGt                                                                                             >  2:354586/9‑251 (MQ=255)
                                                                                                                                                      gttgATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGt                                                                                          <  1:245795/251‑1 (MQ=255)
                                                                                                                                                                     gCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCcg                                                                            <  2:155423/250‑1 (MQ=255)
                                                                                                                                                                           cccacGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCgcacgcac                                                                     <  1:177919/250‑1 (MQ=255)
                                                                                                                                                                                              agcaagAGCTGGCGAATGCGCGTTATAACTACCTGAATAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGTAACAGAATGCTATTGCTGATGGTTATGCGACTGATAGCCCGGCACCAGGCGTTCAGCAAAGATCCGCTCTCGCTACCGCCAGTAATGGTCAt                                                  <  2:254000/251‑1 (MQ=255)
                                                                                                                                                                                                 aagAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAAc                                               >  1:95826/1‑251 (MQ=255)
                                                                                                                                                                                                 aagAGCTGGCGAATGCGCGTTATAACGACCTGATTAATCAGCTGAAGATTAAAGCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCGGAGCAAACCGGTTTCCACTAATCCTGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTTTTGCTTATGGTTATTCGCCGGATAGCCCGTCACCAGTCGGTCAGCATACATCCTGACTCTCTACCTCGAATAACGGTCATAAc                                               <  2:391723/251‑1 (MQ=255)
                                                                                                                                                                                                      cTGGCGAGTGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCATGCTCTGCTGGCACTGACCAATGCGCTGAGCAAACCGGTTTCCACTACTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCATGCTAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCCCCAGTACCGGTCATAACCCttt                                          >  2:328669/1‑251 (MQ=255)
                                                                                                                                                                                                         gCGAATGCGCGTTATAACTACCTGATTAATCAGCGGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCg                                       <  2:226050/251‑1 (MQ=255)
                                                                                                                                                                                                           gAATGCGCGTTCTAGCTACCCGATTAATCAGCTGAATATTAACTAAGCTCTGGGTACGTTGACCGAGCAGGATCTGCTGGCACTGAAAATTGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGACGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCCGCAAACATCAGCACGCACTACCACCAGTACCGATCACAACCCTTTccat                                      >  2:343984/1‑248 (MQ=255)
                                                                                                                                                                                                              tGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGCACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCACACATCCGCACGCACTACCACCAGTAACGGACATAACCCTTTCCGTacct                                  >  2:179997/1‑248 (MQ=255)
                                                                                                                                                                                                                   gTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCAGCAGTAACGGTCATAACCCTTTCCGTAACtgatga                             <  2:268994/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGCTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACtgatga                             <  1:292204/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaCTACCTTATGAATCATCTGAATATTAAATCAGTCACTGGTTCGTTGAACTAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCTCGAATCCGGCAAACGTTTCACCGCAAACCCCGGAACAGAAGGCTCTTGCTGATGGTTATGCGCCTGTTAGCCCGGCCCCAGTCGTTCAGTAAACATCCGCACGCACTACCACCAGTAACTGTCATAACACTTTCCGTAACTGATgacgatg                        <  2:15716/251‑3 (MQ=255)
                                                                                                                                                                                                                          cTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATgacgacgac                      >  1:366980/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTGATTAATCAGATGAATCTTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCCCTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGCTGACGACGACggg                   >  2:19952/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgggg                  >  1:185559/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAAAATCCGCACGCACTACAACCAGAAACGGACATAACCCTTTCCGTAACTGATGACGCCGACgggg                  >  1:353962/1‑251 (MQ=255)
                                                                                                                                                                                                                                gATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGc                 <  1:249859/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCGTCGGcc          <  2:274907/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cgtaaTATTAAATCAGCTCTGGTTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATGGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACTATGGGGCTTCGTCCCCGt      <  2:343076/248‑1 (MQ=255)
                                                                                                                                                                                                                                           tGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTc     >  1:287751/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tatTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGa  >  1:150820/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGA  >  NC_000913/3179138‑3179626

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