Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,511,549 G→C 5.2% R266G (CGA→GGA)  yhfZ ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,511,5490GC 5.2% 317.9 / 5.8 91R266G (CGA→GGA) yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (43/43);  new base C (3/2);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTC  >  NC_000913/3511405‑3511697
                                                                                                                                                |                                                                                                                                                    
gTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTc                                                                                                                                                  >  3:1038440/1‑149 (MQ=255)
gTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTc                                                                                                                                                  >  4:597732/1‑149 (MQ=255)
   gAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATTGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAgg                                                                                                                                               >  4:694047/1‑149 (MQ=255)
    aGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGa                                                                                                                                              <  3:618515/149‑1 (MQ=255)
     gTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGAc                                                                                                                                             >  2:960911/1‑149 (MQ=255)
        tGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCAc                                                                                                                                          <  1:147383/149‑1 (MQ=255)
        tGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCAc                                                                                                                                          >  2:857861/1‑149 (MQ=255)
         ggTTTCCATATTAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACg                                                                                                                                         >  1:818887/1‑149 (MQ=255)
          gTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTTTCTACAACGGTACGCAGCAGTTGTTGCATGGGATAATCATCTACTCCCGTCAGGACCACgg                                                                                                                                        >  4:43539/1‑149 (MQ=255)
                atatCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCgg                                                                                                                                  >  4:430366/1‑149 (MQ=255)
                 tatCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGt                                                                                                                                 >  2:198432/1‑149 (MQ=255)
                 tatCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGt                                                                                                                                 <  1:163781/149‑1 (MQ=255)
                  atCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTg                                                                                                                                <  4:727196/149‑1 (MQ=255)
                    cAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGc                                                                                                                              <  3:730184/149‑1 (MQ=255)
                      gATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCt                                                                                                                            >  4:8311/1‑149 (MQ=255)
                           cTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTaaa                                                                                                                       >  3:919788/1‑149 (MQ=255)
                                 tttAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATcgc                                                                                                                 <  1:674914/149‑1 (MQ=255)
                                 tttAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATcgc                                                                                                                 <  3:449752/149‑1 (MQ=255)
                                  ttAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATcgcg                                                                                                                <  3:380707/149‑1 (MQ=255)
                                   taataaCTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGACGCGCGCGCAGCGAGTGTGTATCAACAACGGCCCGCAGCAGTTGTAGCTTTGGGTAGACATCGACTCGCGTCAGCACCACTGCTTAGGTAGCCGGTAAAAATCGCgg                                                                                                               >  2:1063714/1‑149 (MQ=255)
                                    aataaCTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGa                                                                                                              <  2:474558/149‑1 (MQ=255)
                                      taaCTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATc                                                                                                            >  3:392169/1‑149 (MQ=255)
                                         cTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTc                                                                                                         <  3:558347/149‑1 (MQ=255)
                                         cTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGACAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTc                                                                                                         <  3:43539/149‑1 (MQ=255)
                                                ccTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgc                                                                                                  <  3:290964/149‑1 (MQ=255)
                                                  tGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgcgg                                                                                                >  3:833342/1‑149 (MQ=255)
                                                  tGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgcgg                                                                                                >  1:927202/1‑149 (MQ=255)
                                                   gTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgcggc                                                                                               >  2:637744/1‑149 (MQ=255)
                                                    ttCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgcggcg                                                                                              >  3:779147/1‑149 (MQ=255)
                                                      ccccACTCAATACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCCGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACTACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTc                                                                                            <  2:636283/149‑1 (MQ=255)
                                                        ccACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGc                                                                                          >  3:234050/1‑149 (MQ=255)
                                                        ccACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGc                                                                                          <  3:512618/149‑1 (MQ=255)
                                                        ccACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGCGTAATCATCGACCCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGc                                                                                          >  2:401569/1‑149 (MQ=255)
                                                          actcactACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGACCGGCGTCGCct                                                                                        <  3:631561/149‑1 (MQ=255)
                                                           ctcactACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACTACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCctc                                                                                       <  3:687800/149‑1 (MQ=255)
                                                           ctcactACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCctc                                                                                       <  3:326168/149‑1 (MQ=255)
                                                              actACCCGTTGTTGATGCGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATc                                                                                                            <  1:1027188/125‑1 (MQ=255)
                                                              actACCCGTTGTTGATGCGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATc                                                                                                            >  2:1027188/1‑125 (MQ=255)
                                                               ctACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAAt                                                                                   <  1:411265/149‑1 (MQ=255)
                                                                tACCCGTTGTTGATTGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATc                                                                                  <  3:694047/149‑1 (MQ=255)
                                                                  ttcgaagttgATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATccc                                                                                >  2:807498/7‑149 (MQ=255)
                                                                  cccGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATccc                                                                                <  4:1038440/149‑1 (MQ=255)
                                                                   ccGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAAGCCCa                                                                               <  1:913168/149‑1 (MQ=255)
                                                                     gttgttGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGc                                                                             >  1:857680/1‑149 (MQ=255)
                                                                       tgttgATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCAt                                                                           <  1:280183/149‑1 (MQ=255)
                                                                       tgttgATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCAt                                                                           <  4:522526/149‑1 (MQ=255)
                                                                          tgATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGt                                                                        <  1:1057072/149‑1 (MQ=255)
                                                                           gATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTc                                                                       >  1:490413/1‑149 (MQ=255)
                                                                              gggCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGt                                                                    >  4:415847/1‑149 (MQ=255)
                                                                                  gAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGt                                                                >  3:1072344/1‑149 (MQ=255)
                                                                                     cAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAGTCATCGACTCGCGTCAGGACCACGGCTTCGGGGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGtttt                                                             >  2:208763/1‑149 (MQ=255)
                                                                                       gggCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTcc                                                           <  1:198432/149‑1 (MQ=255)
                                                                                           gtgtTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGAGCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGcca                                                       >  1:193107/1‑149 (MQ=255)
                                                                                              tttATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGccacca                                                    >  2:1048830/1‑149 (MQ=255)
                                                                                                 aTCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGt                                                                                                                                 <  1:1049846/69‑1 (MQ=255)
                                                                                                 aTCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGt                                                                                                                                 >  2:1049846/1‑69 (MQ=255)
                                                                                                 aTCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGt                                                 >  1:48009/1‑149 (MQ=255)
                                                                                                   cTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGcc                                                        >  4:79632/1‑140 (MQ=255)
                                                                                                   cTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGcc                                                        <  3:79632/140‑1 (MQ=255)
                                                                                                    tACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTcc                                              >  3:1042629/1‑149 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAAt                                                                                   >  2:876851/1‑111 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAAt                                                                                   <  1:876851/111‑1 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCa                                             <  3:584549/149‑1 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCa                                             >  1:173538/1‑149 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCa                                             <  1:927177/149‑1 (MQ=255)
                                                                                                     aCAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCa                                             >  3:915499/1‑149 (MQ=255)
                                                                                                       aaCGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGa                                           >  4:296492/1‑149 (MQ=255)
                                                                                                         cGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATc                                         <  3:8311/149‑1 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGa                                           <  4:439334/139‑1 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGa                                           >  3:439334/1‑139 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATc                                 <  4:500530/149‑1 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCCCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTaa                                                                                    >  4:362712/1‑98 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCCCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTaa                                                                                    <  3:362712/98‑1 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCCCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTaa                                                                                    <  3:641404/98‑1 (MQ=255)
                                                                                                                 agcagTTGTTGCATTGGGTAATCATCGACTCCCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTaa                                                                                    >  4:641404/1‑98 (MQ=255)
                                                                                                                       tgttgCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATcgcg                                                                                                                <  4:776602/64‑1 (MQ=255)
                                                                                                                       tgttgCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATcgcg                                                                                                                >  3:776602/1‑64 (MQ=255)
                                                                                                                            cATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCtt                      <  1:1096706/149‑1 (MQ=255)
                                                                                                                                 ggTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACa                 >  1:925644/1‑149 (MQ=255)
                                                                                                                                   tAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAAt               >  1:775253/1‑149 (MQ=255)
                                                                                                                                    aaTCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgc                                                                                                  <  2:238606/65‑1 (MQ=255)
                                                                                                                                    aaTCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAgc                                                                                                  >  1:238606/1‑65 (MQ=255)
                                                                                                                                     atcatcGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGc             <  1:1048830/149‑1 (MQ=255)
                                                                                                                                     atcatcGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGc             <  3:597732/149‑1 (MQ=255)
                                                                                                                                       catcGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGt           <  4:642619/149‑1 (MQ=255)
                                                                                                                                        atcGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGtt          >  4:1061411/1‑149 (MQ=255)
                                                                                                                                           gACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTa       >  2:419013/1‑149 (MQ=255)
                                                                                                                                             cTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATc                           >  4:120315/1‑127 (MQ=255)
                                                                                                                                             cTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATc                           <  3:120315/127‑1 (MQ=255)
                                                                                                                                             cTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTaaa     <  1:637744/149‑1 (MQ=255)
                                                                                                                                              tCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAAc    <  2:834473/149‑1 (MQ=255)
                                                                                                                                                gcgTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAActc  <  2:193107/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTTGAGTCTGGTTTCCATATCAGATGCCTTTCGTTTAATAACTCGGTTCCTGTTCCCCACTCACTACCCGTTGTTGATGGGCGAGCAGGGCGTGTTTATCTACAACGGCACGCAGCAGTTGTTGCATTGGGTAATCATCGACTCGCGTCAGGACCACGGCTTCGGTGGCCTGTAAAAATCGCGGATCGTCTGTGAGCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTC  >  NC_000913/3511405‑3511697

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: