Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,091,507 G→T 5.3% P137P (CCG→CCT hisD → bifunctional histidinal dehydrogenase/ histidinol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,091,5070GT 5.3% 262.0 / 3.4 76P137P (CCG→CCThisDbifunctional histidinal dehydrogenase/ histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

GTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAACGTCGGCG  >  NC_000913/2091364‑2091652
                                                                                                                                               |                                                                                                                                                 
gTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCtctt                                                                                                                                              <  8:113437/149‑1 (MQ=255)
 tGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCtcttc                                                                                                                                             <  1:195665/149‑1 (MQ=255)
      gTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTACGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAAc                                                                                                                                        >  1:124833/1‑149 (MQ=255)
       tAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACg                                                                                                                                       >  5:191689/1‑149 (MQ=255)
           gAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTAt                                                                                                                                   >  1:188928/1‑149 (MQ=255)
           gAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTAt                                                                                                                                   <  7:237959/149‑1 (MQ=255)
           gAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTGCAGTTGGGTGGTATATTCCTGGCGGCGCCGCCCCGCTCTTCTCAACGGTAt                                                                                                                                   <  8:126578/149‑1 (MQ=255)
              tatTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGTTCCGCCCCGCTCTTCTCAACGGTATTaa                                                                                                                                <  4:388692/149‑1 (MQ=255)
              tatTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTTCGCCCCGCtctt                                                                                                                                              >  1:197864/1‑135 (MQ=255)
              tatTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTTCGCCCCGCtctt                                                                                                                                              <  2:197864/135‑1 (MQ=255)
                ttGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATg                                                                                                                              <  2:188928/149‑1 (MQ=255)
                   aaaCCTTCCACACTTCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTg                                                                                                                           >  3:1543/1‑149 (MQ=255)
                       cTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGa                                                                                                                       <  8:194792/149‑1 (MQ=255)
                        ttCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGa                                                                                                                       >  2:405383/1‑148 (MQ=255)
                             cacTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCgg                                                                                                                 >  6:372143/1‑149 (MQ=255)
                                   gcAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGa                                                                                                                       >  4:184685/1‑137 (MQ=255)
                                   gcAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGa                                                                                                                       <  3:184685/137‑1 (MQ=255)
                                     aaaaaCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTAtt                                                                                                         >  4:4008/1‑149 (MQ=255)
                                         aCTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCgg                                                                                                     >  1:49052/1‑149 (MQ=255)
                                                ccgGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTaa                                                                                              <  8:12734/149‑1 (MQ=255)
                                                 cgGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTaa                                                                                              >  1:370383/1‑148 (MQ=255)
                                                   gTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGACTCCGCCCCTCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTaaaa                                                                                            >  3:173746/1‑148 (MQ=255)
                                                       aTGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAgtg                                                                                       <  6:191689/149‑1 (MQ=255)
                                                          tAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAgtggtg                                                                                    >  3:70669/1‑149 (MQ=255)
                                                            gAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGGGTATTGCGGGCTGTAAAAAAGTGgtgct                                                                                  >  5:168275/1‑149 (MQ=255)
                                                            gAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgct                                                                                  >  8:75131/1‑149 (MQ=255)
                                                            gAAACGCAGCCAGGCGTACGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACTGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgct                                                                                  <  2:124833/149‑1 (MQ=255)
                                                               aCGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgtg                                                                               <  8:179793/149‑1 (MQ=255)
                                                                 gCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgtgct                                                                             >  1:253576/1‑149 (MQ=255)
                                                                    gCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTATTGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCAc                                                                          <  8:199848/149‑1 (MQ=255)
                                                                     ccAGGCGTGCGTTGCCAGCAGGTCACTCGTCCGGTAGCTTCAGTTGGGTTGTATTTTGCTGGCGGCTCCGCCCCGCTCTTCTCTACTGTATTAATGCTGGGGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCAcc                                                                         <  1:296061/149‑1 (MQ=255)
                                                                      cAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCAccg                                                                        >  3:388856/1‑149 (MQ=255)
                                                                       aGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCAccgc                                                                       >  8:191825/1‑149 (MQ=255)
                                                                            tgcgtTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGTGTTGTATATTCCTGGCGGCTCCGCCCCGCACTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTAAccgcctcc                                                                   >  6:152526/1‑148 (MQ=255)
                                                                               gtTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTa                                                                                                           <  7:273743/105‑1 (MQ=255)
                                                                               gtTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTa                                                                                                           >  8:273743/1‑105 (MQ=255)
                                                                                   ccAGCAGGTCACGCGTCCGGTAGCTTCAGTAGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTTCTGGCTGTAAAAAACTGGTGCTGTGCTCACCGCCGCCGATTGcc                                                            <  5:152526/148‑1 (MQ=255)
                                                                                     agcagGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGAt                                                         <  7:191825/149‑1 (MQ=255)
                                                                                     agcagGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGAt                                                         <  2:253576/149‑1 (MQ=255)
                                                                                       cagGGCACGCGTCCGGTAGCTTCAGGTGGGTGGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATga                                                       <  4:70669/149‑1 (MQ=255)
                                                                                           tCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACAGCCGCCGATTGCCGATGAGATc                                                   <  5:282044/149‑1 (MQ=255)
                                                                                             aCGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCt                                                 >  3:215922/1‑149 (MQ=255)
                                                                                               gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgt                                                                                >  4:261351/1‑116 (MQ=255)
                                                                                               gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGgtgctgt                                                                                <  3:261351/116‑1 (MQ=255)
                                                                                               gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCttt                                               <  1:6937/149‑1 (MQ=255)
                                                                                               gcgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCttt                                               <  2:370383/149‑1 (MQ=255)
                                                                                                cgTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTa                                              <  5:372143/149‑1 (MQ=255)
                                                                                                     ggTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATgcgg                                         >  5:19537/1‑149 (MQ=255)
                                                                                                      gTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATgcggc                                        >  5:386847/1‑149 (MQ=255)
                                                                                                        aGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATgag                                                      >  2:71317/1‑133 (MQ=255)
                                                                                                        aGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATgag                                                      <  1:71317/133‑1 (MQ=255)
                                                                                                         gCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATgcggcg                                       >  5:194479/1‑147 (MQ=255)
                                                                                                         gCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGgcgc                                      <  3:4008/148‑1 (MQ=255)
                                                                                                          cTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCTCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAg                                    >  7:357952/1‑149 (MQ=255)
                                                                                                          cTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAg                                    <  6:194479/149‑1 (MQ=255)
                                                                                                          cTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAg                                    >  2:167357/1‑149 (MQ=255)
                                                                                                           ttCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGc                                   >  8:253609/1‑149 (MQ=255)
                                                                                                           ttCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGc                                   >  6:248136/1‑149 (MQ=255)
                                                                                                                ttGGGTGGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCGGTAAAAAAGTGGTGCGGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGc                              <  1:262236/149‑1 (MQ=255)
                                                                                                                    gTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTAAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGgtg                          <  6:19537/149‑1 (MQ=255)
                                                                                                                       gTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGtgtg                        >  3:244767/1‑148 (MQ=255)
                                                                                                                        tataTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTTAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGtgtg                        >  4:74613/1‑147 (MQ=255)
                                                                                                                        tataTTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGtgtg                        >  8:283744/1‑147 (MQ=255)
                                                                                                                          tatTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAg                     <  4:388856/148‑1 (MQ=255)
                                                                                                                           atTCCTGGCGGCTCCGCACCGCTCTTCTCAACGGTAGTAATGCTGGCGACTCCGGCGAGTATTGCGGGGTGTTAAAAAGTGGTGCTGTGCTCACCGCCGCCGAGTGCCGATGCGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGa                   >  8:399806/1‑149 (MQ=255)
                                                                                                                               cTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACgtg               >  8:178976/1‑149 (MQ=255)
                                                                                                                               cTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACgtg               >  3:239845/1‑149 (MQ=255)
                                                                                                                                tGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACgtgt              >  5:124838/1‑149 (MQ=255)
                                                                                                                                    gacTCCGCCCCTCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTaa          <  4:173746/147‑1 (MQ=255)
                                                                                                                                      cTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACgtg               <  7:332920/142‑1 (MQ=255)
                                                                                                                                      cTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACgtg               >  8:332920/1‑142 (MQ=255)
                                                                                                                                        ccGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAACGTc      <  2:49052/149‑1 (MQ=255)
                                                                                                                                         cGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTACGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAAc         <  1:286541/145‑1 (MQ=255)
                                                                                                                                         cGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTACGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAAc         >  2:286541/1‑145 (MQ=255)
                                                                                                                                           ccccGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCGGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTCAGGACGTGTTTAACGTcggcg  <  4:1543/149‑1 (MQ=255)
                                                                                                                                            cccTCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAACGTcggcg  <  8:357952/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTGGCAGTAAAGAATATTGAAACCTTCCACACTGCGCAAAAACTGCCGCCGGTAGATGTAGAAACGCAGCCAGGCGTGCGTTGCCAGCAGGTCACGCGTCCGGTAGCTTCAGTTGGGTTGTATATTCCTGGCGGCTCCGCCCCGCTCTTCTCAACGGTATTAATGCTGGCGACTCCGGCGAGTATTGCGGGCTGTAAAAAAGTGGTGCTGTGCTCACCGCCGCCGATTGCCGATGAGATCCTTTATGCGGCGCAGCTGTGCGGTGTGCAGGACGTGTTTAACGTCGGCG  >  NC_000913/2091364‑2091652

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: