Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 734,341 G→A 5.2% pseudogene (122/1434 nt) ybfO → conserved protein, rhs‑like YbfO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913734,3410GA 5.2% 226.5 / 2.1 77pseudogene (122/1434 nt)ybfOconserved protein, rhs‑like YbfO
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGAC  >  NC_000913/734192‑734486
                                                                                                                                                     |                                                                                                                                                 
ggAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcg                                                                                                                                                   <  1:306956/150‑1 (MQ=32)
ggAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcg                                                                                                                                                   <  2:1370501/150‑1 (MQ=32)
ggAGCTGCACCCGAACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCGGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAgcg                                                                                                                                                   <  1:1207088/150‑1 (MQ=17)
  aGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGca                                                                                                                                                 <  2:34464/150‑1 (MQ=33)
  aGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCCGAAGGGGTTATCCGCACGGATGATGAGCGca                                                                                                                                                 <  2:254048/150‑1 (MQ=18)
   gCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcac                                                                                                                                                <  2:602501/150‑1 (MQ=34)
       cACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccac                                                                                                                                            >  1:1200726/1‑150 (MQ=34)
        aCCCGGACAGCACCCTCAGCATGGGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGcacccacc                                                                                                                                           <  1:331408/150‑1 (MQ=21)
            ggACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCTGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGTATGATGAGCGCACCCACCGGTa                                                                                                                                       <  2:1388008/150‑1 (MQ=21)
              aCAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAcc                                                                                                                                     <  1:1356661/150‑1 (MQ=34)
              aCAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTAcc                                                                                                                                     <  1:1062497/150‑1 (MQ=34)
               cAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCa                                                                                                                                    <  2:418212/150‑1 (MQ=34)
                 gCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCAtt                                                                                                                                  >  2:468383/1‑150 (MQ=34)
                  cACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTa                                                                                                                                 >  1:128252/1‑150 (MQ=33)
                  cACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTa                                                                                                                                 <  1:228766/150‑1 (MQ=33)
                   aCCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTAc                                                                                                                                <  1:239062/150‑1 (MQ=32)
                    cccTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACg                                                                                                                               <  2:308193/150‑1 (MQ=32)
                           cATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGACGATGAGCGCACCCACCGGTACCATTACGACAGTCa                                                                                                                        <  1:879865/150‑1 (MQ=17)
                             tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGc                                                                                                                      <  1:1174163/150‑1 (MQ=32)
                             tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGc                                                                                                                      >  2:894434/1‑150 (MQ=32)
                             tgtgGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGc                                                                                                                      <  1:372190/150‑1 (MQ=32)
                                  ccGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCgg                                                                                                                 >  2:443043/1‑150 (MQ=32)
                                    ggATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCt                                                                                                               >  1:169751/1‑150 (MQ=32)
                                    ggATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCt                                                                                                               >  2:929388/1‑150 (MQ=32)
                                     gATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTg                                                                                                              <  2:994811/150‑1 (MQ=32)
                                      aTAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTgg                                                                                                             <  1:110522/150‑1 (MQ=32)
                                       tAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                            <  1:894434/150‑1 (MQ=32)
                                       tAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACACACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                            >  2:859845/1‑150 (MQ=17)
                                        aaCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                           <  1:929388/150‑1 (MQ=32)
                                         aCCGTATCTCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCACACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                          <  2:359139/150‑1 (MQ=17)
                                         aCCGTATCTCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCACACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                          >  1:359139/1‑150 (MQ=17)
                                         aCCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGc                                                                                                          >  1:609409/1‑150 (MQ=32)
                                          ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                           <  2:791639/148‑1 (MQ=32)
                                          ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                           >  1:791639/1‑148 (MQ=32)
                                          ccGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCa                                                                                                         <  2:1489493/150‑1 (MQ=32)
                                             tATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTa                                                                                                      <  1:887116/150‑1 (MQ=32)
                                               tCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTaca                                                                                                    >  1:1055593/1‑150 (MQ=32)
                                                 gCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcg                                                                                                  >  1:13380/1‑150 (MQ=32)
                                                  cccGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACAcgc                                                                                                 <  1:17523/150‑1 (MQ=32)
                                                        aCGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGacac                                                                                           <  2:463100/150‑1 (MQ=32)
                                                             cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAat                                                                                        >  2:342029/1‑148 (MQ=32)
                                                             cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAat                                                                                        <  1:342029/148‑1 (MQ=32)
                                                             cACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACCGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAatat                                                                                      >  2:601420/1‑150 (MQ=17)
                                                               cTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGc                                                                                    <  2:609409/150‑1 (MQ=32)
                                                                tATCTTTACCGGTATGACCGTCAAGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGATCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCa                                                                                   <  1:847798/150‑1 (MQ=17)
                                                                 aTCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCag                                                                                  >  2:1053640/1‑150 (MQ=32)
                                                                   cTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCagag                                                                                >  2:765889/1‑150 (MQ=32)
                                                                        ccGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCt                                                                           <  2:1351760/150‑1 (MQ=32)
                                                                            tATGACCGTCACGGCAGGCGGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTc                                                                       <  2:1355383/150‑1 (MQ=17)
                                                                                  cGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGc                                                                 >  1:311286/1‑150 (MQ=32)
                                                                                             gCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTa                                                      >  1:698568/1‑150 (MQ=34)
                                                                                              cggACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTAc                                                     >  2:285392/3‑150 (MQ=21)
                                                                                                  cGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACC                                                 >  1:1475560‑M1/1‑150 (MQ=255)
                                                                                                    gagaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGCCCCG                                               >  1:1405132‑M1/1‑150 (MQ=255)
                                                                                                    gagaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCG                                               >  1:692496‑M1/1‑150 (MQ=255)
                                                                                                        aaaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACGGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                           <  1:1033999/150‑1 (MQ=17)
                                                                                                        aaaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTgg                                           >  1:29262/1‑150 (MQ=32)
                                                                                                          aCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGttttgc                                         >  1:233173‑M1/1‑144 (MQ=255)
                                                                                                               ccTCATCCCGGAAGGGGTTAACCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAg                                    >  1:1427900/1‑150 (MQ=17)
                                                                                                                  cATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                 >  1:539465/1‑150 (MQ=32)
                                                                                                                   aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                                >  2:611236/1‑150 (MQ=32)
                                                                                                                   aTCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTg                                >  1:32957/1‑150 (MQ=32)
                                                                                                                        ggAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaa                           >  2:1116419/1‑150 (MQ=32)
                                                                                                                            ggggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGttttgccaccctgcggcccagcgg                       <  2:233173‑M1/150‑25 (MQ=255)
                                                                                                                             gggTTATCCGCACGGATGATGAGCGCACCCCCCGGTACCATTACGCCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg                      >  1:293044/1‑150 (MQ=17)
                                                                                                                               gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGgtg                    <  1:1053640/150‑1 (MQ=32)
                                                                                                                               gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGgtg                    >  2:1219498/1‑150 (MQ=32)
                                                                                                                               gTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCACAGGGTGGCAAAACGGgtg                    <  1:765889/150‑1 (MQ=17)
                                                                                                                                 tATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGtgtg                  <  1:443043/150‑1 (MQ=34)
                                                                                                                                  aTCCGCACGGATGATGAGCACACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTgg                 >  1:741010/1‑150 (MQ=21)
                                                                                                                                   tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGc                >  2:1139415/1‑150 (MQ=34)
                                                                                                                                    ccGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGttttgccaccctgcggcccagcgggtcgtaaa               <  1:255977‑M1/150‑33 (MQ=255)
                                                                                                                                     cGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                               <  2:136669/133‑1 (MQ=32)
                                                                                                                                     cGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                               >  1:136669/1‑133 (MQ=32)
                                                                                                                                      gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTTGAAAGCCGCTATCTTTACGACCCGttttgccaccctgcggcccagcgggtcgtaaaga             <  1:1147819‑M1/150‑35 (MQ=255)
                                                                                                                                      gCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGAc             >  1:171827/1‑150 (MQ=35)
                                                                                                                                       cACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACg            <  2:1358294/150‑1 (MQ=37)
                                                                                                                                       cACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACg            <  1:807515/150‑1 (MQ=37)
                                                                                                                                               atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACggg    <  1:323452/150‑1 (MQ=37)
                                                                                                                                                 gaGCACACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGAc  <  2:741010/150‑1 (MQ=37)
                                                                                                                                                     |                                                                                                                                                 
GGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACGGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGCCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGACGTGAACGGGAC  >  NC_000913/734192‑734486

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