Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,716,508 A→T 5.4% intergenic (‑114/‑39) bisC ← / → yiaD biotin sulfoxide reductase/multicopy suppressor of bamB; outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,716,5080AT 5.4% 352.7 / 7.6 129intergenic (‑114/‑39)bisC/yiaDbiotin sulfoxide reductase/multicopy suppressor of bamB; outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (59/62);  new base T (5/2);  total (65/64)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 130 total aligned reads displayed.

TCAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTAAATCTGCTATC  >  NC_000913/3716361‑3716657
                                                                                                                                                   |                                                                                                                                                     
tCAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAAcc                                                                                                                                                     >  2:1682518/1‑150 (MQ=255)
 cAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAAccc                                                                                                                                                    >  1:1390220/1‑150 (MQ=255)
 cAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAAccc                                                                                                                                                    >  1:1173261/1‑150 (MQ=255)
 cAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGATCACGCGGCATAATCTCCCGCCACGGAAccc                                                                                                                                                    >  2:29885/1‑150 (MQ=255)
   gAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGt                                                                                                                                                  >  2:1367048/1‑150 (MQ=255)
     aCGGACTATCGTGTTGATGAGTTTGCCAAACCGAAGACTCATGAAAGTTTAGATGTATTGAGCATTTAAATATTACTTACACATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAAACCTTgg                                                                                                                                                <  1:1066450/150‑1 (MQ=255)
         aaTATCGTGAGGATGAGGTGGCCAATCCTAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTAAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGGTTCCGTGGcggg                                                                                                                                            <  1:64450/150‑4 (MQ=255)
           tatCGTGAGGATGAGTTGGCCAAACAGAAGACACCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGa                                                                                                                                          <  1:118765/150‑1 (MQ=255)
                    gATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAgg                                                                                                                                 >  2:1142963/1‑150 (MQ=255)
                     aTGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGATTAAAAAGGt                                                                                                                                <  1:562634/150‑1 (MQ=255)
                           tGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTaa                                                                                                                          <  1:336263/150‑1 (MQ=255)
                            ggCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAg                                                                                                                         >  2:1136075/1‑150 (MQ=255)
                             gCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAgg                                                                                                                        <  2:772390/150‑1 (MQ=255)
                              ccAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGa                                                                                                                       >  1:46898/1‑150 (MQ=255)
                               cAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGAt                                                                                                                      <  1:474927/150‑1 (MQ=255)
                               cAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGAt                                                                                                                      >  2:1298251/1‑150 (MQ=255)
                                 aaCCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTa                                                                                                                    <  2:1257226/150‑1 (MQ=255)
                                 aaCCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTa                                                                                                                    >  2:1224764/1‑150 (MQ=255)
                                   ccGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAAc                                                                                                                  >  2:388033/1‑150 (MQ=255)
                                   ccGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAAc                                                                                                                  >  1:686281/1‑150 (MQ=255)
                                     gAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACaa                                                                                                                <  2:327352/150‑1 (MQ=255)
                                      aaGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAAt                                                                                                               <  2:795034/150‑1 (MQ=255)
                                       aGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATg                                                                                                              <  2:1390220/150‑1 (MQ=255)
                                       aGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGGTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATg                                                                                                              <  2:1394071/150‑1 (MQ=255)
                                            ccTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATgaagaa                                                                                                         <  2:1334681/150‑1 (MQ=255)
                                            ccTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATgaagaa                                                                                                         >  2:1061780/1‑150 (MQ=255)
                                             cTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACAGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGaaa                                                                                                        >  2:1575567/1‑150 (MQ=255)
                                             cTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTAACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGaaa                                                                                                        >  1:162109/1‑150 (MQ=255)
                                                 aaGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGATCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACgtg                                                                                                    <  1:29885/150‑1 (MQ=255)
                                                  aGGTTAGATGTATTGAGCAGTTAAATATTTCTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACgtgt                                                                                                   >  2:514707/1‑150 (MQ=255)
                                                  aGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACgtgt                                                                                                   <  1:1061780/150‑1 (MQ=255)
                                                   ggTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGtt                                                                                                  >  2:842395/1‑150 (MQ=255)
                                                    gTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGttt                                                                                                 >  1:466149/1‑150 (MQ=255)
                                                     ttAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGTAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTa                                                                                                <  1:723538/150‑1 (MQ=255)
                                                        gATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGTCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAATGATTATCAATGAAGAAACGTGTTTATCt                                                                                             <  1:1099609/150‑1 (MQ=255)
                                                        gATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCt                                                                                             <  1:894836/150‑1 (MQ=255)
                                                            tATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTAtt                                                                                         <  2:899590/150‑1 (MQ=255)
                                                             aTTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTg                                                                                        <  2:1173261/150‑1 (MQ=255)
                                                             aTTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTg                                                                                        >  2:396576/1‑150 (MQ=255)
                                                             aTTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTg                                                                                        <  2:686281/150‑1 (MQ=255)
                                                                gAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTgccg                                                                                     <  1:618641/150‑1 (MQ=255)
                                                                gAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTgccg                                                                                     >  2:1595866/1‑150 (MQ=255)
                                                                gAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTgccg                                                                                     >  2:710379/1‑150 (MQ=255)
                                                                   cAGTTAAATATTACTTACAAATTTCCGAGTATTCCCTGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGATCCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGccgccg                                                                                  >  1:102798/1‑150 (MQ=255)
                                                                   cAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGccgccg                                                                                  >  2:974802/1‑150 (MQ=255)
                                                                      ttAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTag                                                                               >  2:388574/1‑150 (MQ=255)
                                                                        aaaTATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTagtg                                                                             <  2:1396985/150‑1 (MQ=255)
                                                                              tacttacAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGACCCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGgc                                                                       >  1:102348/1‑150 (MQ=255)
                                                                               acttacAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGgcg                                                                      >  1:41522/1‑150 (MQ=255)
                                                                                cttacAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGgcgc                                                                     <  2:932630/150‑1 (MQ=255)
                                                                                 ttacAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGct                                                                    >  2:565916/1‑150 (MQ=255)
                                                                                   acAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGctct                                                                  <  2:1201034/150‑1 (MQ=255)
                                                                                    cAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATATAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTg                                                                 >  1:1567599/1‑150 (MQ=255)
                                                                                       aTTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCg                                                              <  2:486939/150‑1 (MQ=255)
                                                                                          tCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTa                                                           >  2:1530944/1‑150 (MQ=255)
                                                                                            cGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATc                                                         >  2:557776/1‑150 (MQ=255)
                                                                                            cGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATc                                                         >  2:1189082/1‑150 (MQ=255)
                                                                                            cGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATc                                                         >  1:832382/1‑150 (MQ=255)
                                                                                            cGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATc                                                         >  1:879055/1‑150 (MQ=255)
                                                                                             gAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCt                                                        >  1:1486484/1‑150 (MQ=255)
                                                                                             gAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCt                                                        >  2:1705583/1‑150 (MQ=255)
                                                                                                tATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGc                                                     <  2:660885/150‑1 (MQ=255)
                                                                                                 aTTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCt                                                    <  2:162109/150‑1 (MQ=255)
                                                                                                 aTTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCt                                                    <  1:238564/150‑1 (MQ=255)
                                                                                                    cccAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGca                                                 >  2:1702717/1‑150 (MQ=255)
                                                                                                    cccAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGca                                                 <  1:1508538/150‑1 (MQ=255)
                                                                                                     ccAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGTATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGcac                                                >  1:491949/1‑150 (MQ=255)
                                                                                                     ccAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTAATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGACTGcac                                                >  1:112251/1‑150 (MQ=255)
                                                                                                      cAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGcaca                                               <  1:745994/150‑1 (MQ=255)
                                                                                                      cAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGcaca                                               <  2:466149/150‑1 (MQ=255)
                                                                                                      cAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGAAGTGAGTGGCGCTCTGGCGGTATCTGGCTGcaca                                               <  2:1680733/150‑1 (MQ=255)
                                                                                                      cAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCAGAGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATGGCCGCCGTAGTGAGTGGCGCTCTGGCGGTAGCTGGCTGcacc                                               >  1:331551/1‑149 (MQ=255)
                                                                                                       aGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACaa                                              >  1:985925/1‑150 (MQ=255)
                                                                                                        ggAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAAc                                             <  1:1367048/150‑1 (MQ=255)
                                                                                                        ggAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAAc                                             <  2:442762/150‑1 (MQ=255)
                                                                                                         gAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACt                                            <  1:1341545/150‑1 (MQ=255)
                                                                                                             aatCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAAcc                                        <  1:710379/150‑1 (MQ=255)
                                                                                                              atCTTCATGTTCACGCGGCATAATCTACCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAAccc                                       <  2:46898/150‑1 (MQ=255)
                                                                                                               tCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCt                                      >  2:1450009/1‑150 (MQ=255)
                                                                                                                 ttCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTa                                    <  1:183206/150‑1 (MQ=255)
                                                                                                                    aTGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTacac                                 <  1:1136075/150‑1 (MQ=255)
                                                                                                                       ttCACGCGGCATAATCTCGCGGCTCGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAATCGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATGTGGCTGCACAACTAACCCTTACACCgg                              <  1:1409548/150‑1 (MQ=255)
                                                                                                                        tCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGAAGAACAAGGAAGCAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCGGCCCAACTAGCCCTTACACCGGc                             >  1:344656/1‑150 (MQ=255)
                                                                                                                         cACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCg                            <  1:388574/150‑1 (MQ=255)
                                                                                                                            gcgGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAAc                         >  2:842743/1‑150 (MQ=255)
                                                                                                                             cgGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAAcg                        <  1:388033/150‑1 (MQ=255)
                                                                                                                             cgGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAAcg                        >  2:767164/1‑150 (MQ=255)
                                                                                                                                cATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGa                     >  1:593715/1‑150 (MQ=255)
                                                                                                                                 aTAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGaa                    >  2:1147341/1‑150 (MQ=255)
                                                                                                                                 aTAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGaa                    <  1:1595866/150‑1 (MQ=255)
                                                                                                                                   aaTCTCCCGCCACGGAACCCGTGGCAAAGAATAAAATGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGTGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGc                  <  1:901504/150‑1 (MQ=255)
                                                                                                                                     tctcCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAg                <  1:963497/150‑1 (MQ=255)
                                                                                                                                       ctttgCCACGGGTTCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAgg               >  2:309138/5‑149 (MQ=255)
                                                                                                                                       ctttgCCACGGGTTCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAgg               <  1:309138/145‑1 (MQ=255)
                                                                                                                                        tttgCCACGGGTTCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTa             >  2:759721/4‑150 (MQ=255)
                                                                                                                                        tttgCCACGGGTTCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTa             >  1:833035/4‑150 (MQ=255)
                                                                                                                                        cccGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTa             <  1:1189082/150‑1 (MQ=255)
                                                                                                                                          cGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTaaa           >  2:958251/1‑150 (MQ=255)
                                                                                                                                            ccACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTAAATc         <  1:1298251/150‑1 (MQ=255)
                                                                                                                                                   aCCCGTGGCAAAGAATAAAAAGGTTATTTAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTAAATCTGCTATc  >  2:152941/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGAACGGAATATCGTGAGGATGAGTTGGCCAAACCGAAGACTCCTGAAAGGTTAGATGTATTGAGCAGTTAAATATTACTTACAAATTTCCGAGTATTCCCAGGAATAATCTTCATGTTCACGCGGCATAATCTCCCGCCACGGAACCCGTGGCAAAGAATAAAAAGGTTATTAAGGATTAACAATGAAGAAACGTGTTTATCTTATTGCCGCCGTAGTGAGTGGCGCTCTGGCGGTATCTGGCTGCACAACTAACCCTTACACCGGCGAACGCGAAGCAGGTAAATCTGCTATC  >  NC_000913/3716361‑3716657

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