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→AAA) tolCouter 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

CATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGA  >  NC_000913/3179138‑3179626
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
cATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTTTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑c                                                                                                                                                                                                                                                <  1:343047/251‑1 (MQ=255)
     tAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTgg                                                                                                                                                                                                                                           >  2:115003/1‑251 (MQ=255)
       gCGTCGTGCAGACCGTGCGTTCCTCCTCCACCAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGt                                                                                                                                                                                                                                         <  1:254367/251‑1 (MQ=255)
              gCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTg                                                                                                                                                                                                                                   <  1:141083/250‑1 (MQ=255)
                  aCCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACg                                                                                                                                                                                                                               >  1:254000/1‑250 (MQ=255)
                   ccGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAg                                                                                                                                                                                                                             <  2:71691/251‑1 (MQ=255)
                    cGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATGAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAAGCAG‑CTCTGGGTACGTTGAACGAGc                                                                                                                                                                                                                            <  1:28817/251‑1 (MQ=255)
                              tcctTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATctgct                                                                                                                                                                                                                  >  1:209394/1‑251 (MQ=255)
                                        acaTTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAAc                                                                                                                                                                                                        <  1:366963/251‑1 (MQ=255)
                                               tGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGATTGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATgcgc                                                                                                                                                                                                 >  2:381542/1‑251 (MQ=255)
                                                 cATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCCTTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTg                                                                                                                                                                                               >  2:85187/1‑251 (MQ=255)
                                                 cATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCGGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATTCGCTg                                                                                                                                                                                               <  2:389409/251‑1 (MQ=255)
                                                      aTCAGTAGCAGTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGGACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGc                                                                                                                                                                                            <  1:80522/249‑1 (MQ=255)
                                                               aTTAACGCCTACAACCAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTc                                                                                                                                                                                 >  2:275879/1‑251 (MQ=255)
                                                                   aCGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACt                                                                                                                                                                             <  1:204357/251‑1 (MQ=255)
                                                                             acaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGaa                                                                                                                                                                    >  1:336692/1‑250 (MQ=255)
                                                                                         ttCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCg                                                                                                                                                        <  1:93265/250‑1 (MQ=255)
                                                                                             gCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGCGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                                                                                                                   >  2:131560/1‑251 (MQ=255)
                                                                                                       cATTAGACGCGGTGGAAGCGGGCTGCTCGGTCGGGACGCGTACCATTGGTGATGTGTTGGATGCGACCTCCACGTTGTACAACGCCAAGCAAGAGCTGTAGAATGCGCGTTATAACTACCTGAGTAAGCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTTAGCAATACGGTGTCTACTAATCCGGAAAACGTTGCACCGCCACCGCCGTAAc                                                                                                                                           <  2:372924/249‑1 (MQ=255)
                                                                                                       cATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACa                                                                                                                                          <  1:321165/250‑1 (MQ=255)
                                                                                                       cATTAGACGCGATGGAAGCGGGCTACGCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAg                                                                                                                                         <  1:275879/251‑1 (MQ=255)
                                                                                                           aGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAAt                                                                                                                                      <  1:150107/250‑1 (MQ=255)
                                                                                                                  aTGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGc                                                                                                                              <  1:115003/251‑1 (MQ=255)
                                                                                                                     gAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGa                                                                                                                           <  2:171501/251‑1 (MQ=255)
                                                                                                                           ggCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTa                                                                                                                     <  1:261998/251‑1 (MQ=255)
                                                                                                                                   ggtcggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTg                                                                                                             <  1:239357/251‑1 (MQ=255)
                                                                                                                                      cggtACGCGTACCATTGTTGATGTGTTGGCTGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGATTCCCCTACTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGAt                                                                                                           >  2:384170/1‑250 (MQ=255)
                                                                                                                                                   attgttactgtgTTGGATGCGACCACCACGTTGTAAACCGCCTAGCAAGAGCTGGCGGATGCGCGTTATAACTACCTGATTAATCAGCTGAATCTTAACTCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCCAACGCCGGAACCGAATGCTATTGCTGATAGTTATGCGCCTGATAGCCCGGCCCCAGt                                                                                             >  2:354586/9‑251 (MQ=255)
                                                                                                                                                      gttgATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGt                                                                                          <  1:245795/251‑1 (MQ=255)
                                                                                                                                                                     gCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCcg                                                                            <  2:155423/250‑1 (MQ=255)
                                                                                                                                                                           cccacGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCgcacgcac                                                                     <  1:177919/250‑1 (MQ=255)
                                                                                                                                                                                              agcaagAGCTGGCGAATGCGCGTTATAACTACCTGAATAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGTAACAGAATGCTATTGCTGATGGTTATGCGACTGATAGCCCGGCACCAGGCGTTCAGCAAAGATCCGCTCTCGCTACCGCCAGTAATGGTCAt                                                  <  2:254000/251‑1 (MQ=255)
                                                                                                                                                                                                 aagAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAAc                                               >  1:95826/1‑251 (MQ=255)
                                                                                                                                                                                                 aagAGCTGGCGAATGCGCGTTATAACGACCTGATTAATCAGCTGAAGATTAAAGCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCGGAGCAAACCGGTTTCCACTAATCCTGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTTTTGCTTATGGTTATTCGCCGGATAGCCCGTCACCAGTCGGTCAGCATACATCCTGACTCTCTACCTCGAATAACGGTCATAAc                                               <  2:391723/251‑1 (MQ=255)
                                                                                                                                                                                                      cTGGCGAGTGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCATGCTCTGCTGGCACTGACCAATGCGCTGAGCAAACCGGTTTCCACTACTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCATGCTAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCCCCAGTACCGGTCATAACCCttt                                          >  2:328669/1‑251 (MQ=255)
                                                                                                                                                                                                         gCGAATGCGCGTTATAACTACCTGATTAATCAGCGGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCg                                       <  2:226050/251‑1 (MQ=255)
                                                                                                                                                                                                           gAATGCGCGTTCTAGCTACCCGATTAATCAGCTGAATATTAACTAAG‑CTCTGGGTACGTTGACCGAGCAGGATCTGCTGGCACTGAAAATTGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGACGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCCGCAAACATCAGCACGCACTACCACCAGTACCGATCACAACCCTTTccat                                      >  2:343984/1‑248 (MQ=255)
                                                                                                                                                                                                              tGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGCACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCACACATCCGCACGCACTACCACCAGTAACGGACATAACCCTTTCCGTacct                                  >  2:179997/1‑248 (MQ=255)
                                                                                                                                                                                                                   gTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCAGCAGTAACGGTCATAACCCTTTCCGTAACtgatga                             <  2:268994/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTTATAACTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGCTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACtgatga                             <  1:292204/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaCTACCTTATGAATCATCTGAATATTAAATCAGTCACTGG‑TTCGTTGAACTAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCTCGAATCCGGCAAACGTTTCACCGCAAACCCCGGAACAGAAGGCTCTTGCTGATGGTTATGCGCCTGTTAGCCCGGCCCCAGTCGTTCAGTAAACATCCGCACGCACTACCACCAGTAACTGTCATAACACTTTCCGTAACTGATgacgatg                        <  2:15716/251‑3 (MQ=255)
                                                                                                                                                                                                                          cTACCTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATgacgacgac                      >  1:366980/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTGATTAATCAGATGAATCTTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCCCTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGCTGACGACGACggg                   >  2:19952/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgggg                  >  1:185559/1‑251 (MQ=255)
                                                                                                                                                                                                                              cTGATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAAAATCCGCACGCACTACAACCAGAAACGGACATAACCCTTTCCGTAACTGATGACGCCGACgggg                  >  1:353962/1‑251 (MQ=255)
                                                                                                                                                                                                                                gATTAATCAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGc                 <  1:249859/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGCTGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCGTCGGcc          <  2:274907/250‑1 (MQ=255)
                                                                                                                                                                                                                                          cgtaaTATTAAATCAG‑CTCTGGTTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATGGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACTATGGGGCTTCGTCCCCGt      <  2:343076/248‑1 (MQ=255)
                                                                                                                                                                                                                                           tGAATATTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTc     >  1:287751/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tatTAAATCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGa  >  1:150820/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CATCGTAGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAG‑CTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGA  >  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: ‑