Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,302,834 G→C 6.5% intergenic (+21/‑65) oppA → / → oppB oligopeptide ABC transporter periplasmic binding protein/oligopeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,302,8340GC 6.5% 186.5 / 4.7 62intergenic (+21/‑65)oppA/oppBoligopeptide ABC transporter periplasmic binding protein/oligopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base C (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGA  >  NC_000913/1302685‑1302977
                                                                                                                                                     |                                                                                                                                               
gATTCGGCCATTGTTGCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGTGGCAGGAg                                                                                                                                                 >  2:202193/1‑150 (MQ=255)
  ttCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAgtg                                                                                                                                               <  1:576804/150‑1 (MQ=255)
        cATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGGGGGGTGGTGGCTATACCGGCAAAGAGCCGCTGGATAATACCTATACCCGGAATATGTACATTGGGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGc                                                                                                                                         <  2:173542/150‑1 (MQ=255)
           tGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGcccc                                                                                                                                      <  1:1230782/150‑1 (MQ=255)
            gTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCa                                                                                                                                     <  2:175757/150‑1 (MQ=255)
            gTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGAGCAGGACACTCCTGCCCCa                                                                                                                                     <  2:151905/150‑1 (MQ=255)
               ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                  <  2:616772/150‑1 (MQ=255)
               ccTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACgg                                                                                                                                  <  1:947758/150‑1 (MQ=255)
               actctttATTACTACGTGAATGCGCGTCTTTTGAAACCGTTGGGTGGTGGCAATACCGGCATACATCCGCTGGATAATACCTATACCAGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGGCCTGCTCCACgg                                                                                                                                  <  2:25718/146‑1 (MQ=255)
                 tGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGgtg                                                                                                                                <  2:572462/150‑1 (MQ=255)
                 tGTTTATTACTACGGGAAGGCGCGTCTGGTGAAACCGTGGGTGGGTGGCTAGACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGGGAAGCACTAATGGCAATACGTGGGGCAGGAGTGGCCTGCTCCACGgtg                                                                                                                                <  1:540333/150‑1 (MQ=255)
                      attaCTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGa                                                                                                                           <  2:253/150‑1 (MQ=255)
                            acGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTa                                                                                                                     <  1:51814/150‑1 (MQ=255)
                             cGTGAATGCGCGTCTGGTGAAACCGTGGGTTTGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTAt                                                                                                                    <  1:303420/150‑1 (MQ=255)
                              gTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATc                                                                                                                   <  2:1073899/150‑1 (MQ=255)
                                     cgcgTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTAc                                                                                                            >  1:698654/1‑150 (MQ=255)
                                     cgcgTCTGGTGAAACCGTGGGTGGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTCAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCAGGGTGACTGATATATATCTCATTAc                                                                                                            >  2:632891/1‑150 (MQ=255)
                                                 aaCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAgg                                                                                                >  2:8734/1‑150 (MQ=255)
                                                 aaCCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCCGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAgg                                                                                                >  2:265629/1‑150 (MQ=255)
                                                   ccGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGcc                                                                                              >  2:1199433/1‑150 (MQ=255)
                                                        ggTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAg                                                                                         >  2:1181399/1‑150 (MQ=255)
                                                                gCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGTGAAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCa                                                                                 >  2:1196119/1‑150 (MQ=255)
                                                                 cTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGATTATGTACATTGTGAAGCACAAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCaa                                                                                <  1:428862/150‑1 (MQ=255)
                                                                   ataCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                              <  1:401425/150‑1 (MQ=255)
                                                                   ataCCGGCAAAGATCCGCGGGATAATACCTATACCCGGAATATGTACAGTGTGAAGCACTAATGGCAATACGTGGTGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATg                                                                              <  2:205379/150‑1 (MQ=255)
                                                                    taCCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGt                                                                             >  1:951319/1‑150 (MQ=255)
                                                                        ggCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaa                                                                         >  1:514179/1‑150 (MQ=255)
                                                                        ggCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTaaa                                                                         >  1:412683/1‑150 (MQ=255)
                                                                          cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAt                                                                       >  1:1138776/1‑150 (MQ=255)
                                                                          cAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCATGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAt                                                                       >  1:1229471/1‑150 (MQ=255)
                                                                           aaaGATCCGCTTGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACTGAACGCACAGGCCAGAAGGTAGGGCAATGTTAAAAtt                                                                      <  2:1036814/150‑1 (MQ=255)
                                                                           aaaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAtt                                                                      >  2:234498/1‑150 (MQ=255)
                                                                            aaGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAAttt                                                                     >  1:71262/1‑150 (MQ=255)
                                                                              gATCCGCTGGATAAAACCTATACCCGGAATATGTACATTGTGAAGCACAATTGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTAt                                                                   <  2:1247557/150‑1 (MQ=255)
                                                                                tCCGCTGGATAATACCTATACCCGGAATATGTACATTGGGAACCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTc                                                                 <  1:62203/150‑1 (MQ=255)
                                                                                 ccGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCt                                                                <  1:302822/150‑1 (MQ=255)
                                                                                  cGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTa                                                               >  2:104968/1‑150 (MQ=255)
                                                                                   ccTGGATAATACCTAAACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGATGTTCCGTCGCCACGGGGTCTTAATTTTAACGCATTACATAAGGCACATGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTAc                                                              <  2:709570/149‑1 (MQ=255)
                                                                                       gATAATACCTATACCCGTAATCTGAACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCGCGACGGTGTCTGATTTATATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCg                                                          <  1:1212951/150‑1 (MQ=255)
                                                                                        atactaCCTATACCCGGAATATGTACATTGTGAAGCACTAATGTCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCCCAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGc                                                         <  2:477658/150‑1 (MQ=255)
                                                                                         taataCCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCt                                                        >  1:39980/1‑150 (MQ=255)
                                                                                           ataCCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGt                                                      >  1:8156/1‑150 (MQ=255)
                                                                                              ccTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGTCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTg                                                   >  1:917455/1‑150 (MQ=255)
                                                                                                 ataCCCGGAATATGTACATTGTGAAGCACTAATGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAg                                               <  1:1037274/150‑1 (MQ=255)
                                                                                                          atatGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCg                                       <  2:490390/150‑1 (MQ=255)
                                                                                                           tatGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGa                                      >  2:1276226/1‑150 (MQ=255)
                                                                                                            atGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGAc                                     <  2:1221129/150‑1 (MQ=255)
                                                                                                             tGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACg                                    <  1:1199433/150‑1 (MQ=255)
                                                                                                             tGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACg                                    >  1:345405/1‑150 (MQ=255)
                                                                                                              gTACATTGTGAAGCACTACTGGCAATAAGTGGGGCTGGAGTGTCGTGCTCCAGGGTGGCTGATTTTTATCGCATTACAGAAGGCCCACGCCAGAAGGTAGGGCAATGTTAAATTTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGc                                   <  2:181875/150‑1 (MQ=255)
                                                                                                              gTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATATTTATCGCATTACAGAAGGCACAGGCAAGAAGGTAGGGAAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATGCCGACGc                                   >  1:979018/1‑150 (MQ=255)
                                                                                                                 cATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCGACGTCGCTGTCTGGAAGCGATTCCGACGCtat                                >  2:960523/1‑150 (MQ=255)
                                                                                                                   ttGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCtattt                              >  1:1006547/1‑150 (MQ=255)
                                                                                                                   ttGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCAACGGTGTCTGATTTTTATCGCATTACAGCAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATACCGACGCtattt                              >  2:1001721/1‑150 (MQ=255)
                                                                                                                         aGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCtt                        >  1:890714/1‑150 (MQ=255)
                                                                                                                               aaTGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACt                  <  1:687556/150‑1 (MQ=255)
                                                                                                                               aaTGGCAATACGTGGGGCAGGAGTGTCCTGCACCACGGTGTCTGATGTTGATCGCATTTAAGAAGGCCCAGGCCAGAAGGTAGGGCAATGATAAAATTTATTCTCCGTCGCTGTCTGGAAGCGAGTCCGAGGATATTTTTTGTTATTACt                  >  2:841364/1‑150 (MQ=255)
                                                                                                                                  ggCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAACGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTAtt               <  2:43994/150‑1 (MQ=255)
                                                                                                                                     aaTACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCg            <  1:118892/150‑1 (MQ=255)
                                                                                                                                        ccGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTc         >  2:162907/2‑150 (MQ=255)
                                                                                                                                        ccGTGGAGCAGGACACTCCTGCCCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTc         >  2:133071/2‑150 (MQ=255)
                                                                                                                                            ggggCAGGAGTGTCCTGCGCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACCGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTa     <  1:122836/150‑1 (MQ=255)
                                                                                                                                               gCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTatga  >  1:859203/1‑150 (MQ=255)
                                                                                                                                               gCAGGAGTGTCCTGCTCCACGGTGTCTGATATTTATCGCATTACAGAAGGCCCAGGACTGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAATCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTATTTttga  >  2:401438/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GATTCGGCCATTGTTCCTGTTTATTACTACGTGAATGCGCGTCTGGTGAAACCGTGGGTTGGTGGCTATACCGGCAAAGATCCGCTGGATAATACCTATACCCGGAATATGTACATTGTGAAGCACTAATGGCAATACGTGGGGCAGGAGTGTCCTGCTCCACGGTGTCTGATTTTTATCGCATTACAGAAGGCACAGGCCAGAAGGTAGGGCAATGTTAAAATTTATTCTACGTCGCTGTCTGGAAGCGATTCCGACGCTATTTATTCTTATTACTATTTCGTTCTTTATGA  >  NC_000913/1302685‑1302977

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