Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,481,904 (A)8→7 coding (273/1539 nt) emrY ← putative multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,481,8970A.96.9% 292.5 / 3.7 66coding (280/1539 nt)emrYputative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (33/30);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGCCA  >  NC_000913_3_pae_tpiA/2481753‑2482042
                                                                                                                                                |                                                                                                                                                 
aTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaa                                                                                                                                              <  1:450255/150‑1 (MQ=255)
 tGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGACCGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaa                                                                                                                                             <  1:600092/150‑1 (MQ=255)
   ttCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAAAACTATAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaagt                                                                                                                                           <  2:268075/150‑3 (MQ=255)
        ttttCTGACGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATAACAGCACATCAAAATTCGTCGATAGGCTACACATTAATGAAGACAGGGAAAAAAAGTGACTg                                                                                                                                     >  2:970492/1‑150 (MQ=255)
           tCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGACTTACCGCCCCAGCCCTTAACCCATGACTGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGaaa                                                                                                                                  <  3:269939/150‑1 (MQ=255)
           tCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGaaa                                                                                                                                  <  1:261737/150‑1 (MQ=255)
           tCTGGCGGATAATTCCTTAATAATAAACTCAAAAAAAGTGAAAAAAACGGCCCCGACATTAACCCATGAACGACTCTAAAAAACATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATAAAAAAAAAAAAAAAAAGTGACTGaaa                                                                                                                                  <  2:216214/150‑1 (MQ=255)
               gCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTaa                                                                                                                              >  3:236727/1‑150 (MQ=255)
                cGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTaaa                                                                                                                             <  4:224062/150‑1 (MQ=255)
                      aTTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAaataat                                                                                                                       <  1:644960/150‑1 (MQ=255)
                           ttAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                                                                                                                  >  4:167820/1‑150 (MQ=255)
                           ttAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                                                                                                                  >  4:167627/1‑150 (MQ=255)
                            aaataataaACTCTGTGACAGCGGAATTAACGGCCCCGCCATTAACCCCTGAACGACCCTAAAAACTACCAGCACATCAAGATTGGTCGACAGGCTACACATTAATGAAGACAAAGAACAAAAGTGACTGAAAGTAAAAATAATCTTCAt                                                                                                                 >  2:400792/2‑150 (MQ=255)
                                aataaACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                                                                                                             >  3:323888/1‑150 (MQ=255)
                                aataaACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                                                                                                             >  3:484472/1‑150 (MQ=255)
                                 ataaACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGcc                                                                                                            >  2:47660/1‑150 (MQ=255)
                                         tgtgACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt                                                                                                    >  3:42475/1‑150 (MQ=255)
                                         tgtgACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt                                                                                                    <  3:112635/150‑1 (MQ=255)
                                          ataaCAGTAGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATAAAAACAAAAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtg                                                                                                   <  4:49417/147‑1 (MQ=255)
                                           tgACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgt                                                                                                  >  1:195841/1‑150 (MQ=255)
                                           tgACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgt                                                                                                  <  2:1391/150‑1 (MQ=255)
                                            gACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCCTAATTCGCCTATTCTTtgtg                                                                                                 >  2:640945/1‑150 (MQ=255)
                                                 tGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAc                                                                                            >  3:70008/1‑150 (MQ=255)
                                                       tAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCa                                                                                      >  3:345349/1‑150 (MQ=255)
                                                          cGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAACTAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAgca                                                                                   >  3:115186/1‑148 (MQ=255)
                                                             ccccGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACa                                                                                <  1:628757/150‑1 (MQ=255)
                                                              cccGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAg                                                                               <  4:323888/150‑1 (MQ=255)
                                                                       aaCCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCaa                                                                      >  2:979323/1‑150 (MQ=255)
                                                                          cccTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATgg                                                                   >  2:579220/1‑150 (MQ=255)
                                                                           ccTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGc                                                                  <  4:387650/150‑1 (MQ=255)
                                                                                 cGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGc                                                            <  4:236727/150‑1 (MQ=255)
                                                                                       tAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                                                                                                             >  2:756473/1‑95 (MQ=255)
                                                                                       tAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                                                                                                             <  1:756473/95‑1 (MQ=255)
                                                                                        aaaaaaTATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCaa                                                                      <  1:46385/133‑1 (MQ=255)
                                                                                        aaaaaaTATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCaa                                                                      >  2:46385/1‑133 (MQ=255)
                                                                                         aaaaaTATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCaa                                                    >  3:355058/1‑150 (MQ=255)
                                                                                                tCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGTCTATTCTTTGTGCCAACTTGCCAGTAACAGGGATCGCAATGGCATTTGCTACAACTAACGAggg                                             >  4:579263/1‑149 (MQ=255)
                                                                                                 cAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTg                                            <  3:167627/150‑1 (MQ=255)
                                                                                                 cAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTg                                            <  3:167820/150‑1 (MQ=255)
                                                                                                       aTCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGAt                                                                           >  1:809201/1‑113 (MQ=255)
                                                                                                       aTCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGAt                                                                           <  2:809201/113‑1 (MQ=255)
                                                                                                       aTCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAAcc                                      >  4:82189/1‑150 (MQ=255)
                                                                                                       aTCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAAcc                                      >  1:882792/1‑150 (MQ=255)
                                                                                                         cAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgt                                                                                                  <  3:136251/88‑1 (MQ=255)
                                                                                                         cAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTtgt                                                                                                  >  4:136251/1‑88 (MQ=255)
                                                                                                         cAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAgg                                                                              >  3:602272/1‑108 (MQ=255)
                                                                                                         cAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAgg                                                                              <  4:602272/108‑1 (MQ=255)
                                                                                                         cAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCa                                    <  3:401445/150‑1 (MQ=255)
                                                                                                          aaGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAg                                   <  1:853461/150‑1 (MQ=255)
                                                                                                           aGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAgg                                  >  4:35587/1‑150 (MQ=255)
                                                                                                                    gATAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGt                         >  1:812565/1‑150 (MQ=255)
                                                                                                                     aTAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAAt                                                                     <  3:392149/105‑1 (MQ=255)
                                                                                                                     aTAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAAt                                                                     >  4:392149/1‑105 (MQ=255)
                                                                                                                      tAGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCt                       <  3:326584/150‑1 (MQ=255)
                                                                                                                       aGGCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTg                      <  4:42475/150‑1 (MQ=255)
                                                                                                                         gCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATg                                                                    >  1:472573/1‑102 (MQ=255)
                                                                                                                         gCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATg                                                                    <  2:472573/102‑1 (MQ=255)
                                                                                                                         gCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                    <  2:812565/150‑1 (MQ=255)
                                                                                                                         gCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                    >  2:836102/1‑150 (MQ=255)
                                                                                                                         gCTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                    <  4:484472/150‑1 (MQ=255)
                                                                                                                          cTACACATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTg                   <  3:353164/150‑1 (MQ=255)
                                                                                                                            acacATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                    <  1:793910/147‑1 (MQ=255)
                                                                                                                            acacATTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGtt                    >  2:793910/1‑147 (MQ=255)
                                                                                                                              acaTTAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGATACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGc               >  1:258164/1‑150 (MQ=255)
                                                                                                                                  tAATGAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTccc           >  4:20436/1‑150 (MQ=255)
                                                                                                                                      gAAGACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGaa       <  4:378093/150‑1 (MQ=255)
                                                                                                                                         gACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                                                                                                               >  4:102528/1‑43 (MQ=255)
                                                                                                                                         gACAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                                                                                                               <  3:102528/43‑1 (MQ=255)
                                                                                                                                           cAGCGAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGCCa  <  1:836102/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCTGCCAGTAACAGGGATCGCAATGGCATTTGCTACACCAAACGAGGTGATAACCCAGGTGCCTTCGTCTGTTGATGCTCCCAGAAAGCCA  >  NC_000913_3_pae_tpiA/2481753‑2482042

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