Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,850,765 T→G 6.3% intergenic (‑39/+37) uhpA ← / ← ilvN DNA‑binding transcriptional activator UhpA/acetohydroxy acid synthase I subunit IlvN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,7650TG 6.3% 282.3 / 10.2 96intergenic (‑39/+37)uhpA/ilvNDNA‑binding transcriptional activator UhpA/acetohydroxy acid synthase I subunit IlvN
Reads supporting (aligned to +/- strand):  ref base T (42/48);  new base G (4/2);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCT  >  NC_000913/3850623‑3850914
                                                                                                                                              |                                                                                                                                                     
gAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTCTCCTGGCGGGTAAAAAAAAAACGCGCTTACCTTAACGaa                                                                                                                                                >  1:930056/1‑149 (MQ=255)
 aaCCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATa                                                                                                                                               >  1:1195947/1‑150 (MQ=255)
 aaCCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATa                                                                                                                                               <  1:1270740/150‑1 (MQ=255)
    cAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAgc                                                                                                                                            <  1:1520695/150‑1 (MQ=255)
     aaaCTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAgcg                                                                                                                                           <  1:411764/150‑1 (MQ=255)
      aaCTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAgcgc                                                                                                                                          <  1:1289146/150‑1 (MQ=255)
       aCTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTATCGTTAAGGTAAGcgcg                                                                                                                                         <  2:869988/150‑1 (MQ=255)
          cGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTATCGTTAAGGTAAGcgcgtat                                                                                                                                      <  1:808607/150‑4 (MQ=255)
            gCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGAtgtt                                                                                                                                    >  1:1142762/1‑150 (MQ=255)
               aCTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCATAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTGCTGTCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATgttgtt                                                                                                                                 <  2:1504438/150‑1 (MQ=255)
                cTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTc                                                                                                                                <  1:452637/150‑1 (MQ=255)
                 tACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGGCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCa                                                                                                                               <  1:1187488/150‑1 (MQ=255)
                  aCCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCaa                                                                                                                              <  1:204123/150‑1 (MQ=255)
                    cTGAAAATCAGGTTCCAGCCCCAGCAGCTGCGCAACGCCGGAGCGGACGATGAGGAGATCGTCTATAAGGGAAACGGTGATCATGGTCTTGTCCTGGCGGTAAATAAAAATAAGCCCTTACCTAAACGATCAGAGGGATGTTGTTAAAGc                                                                                                                            >  2:698443/1‑150 (MQ=255)
                      gCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAAAACGAGATTACCTTAACGATAAGCGCGCTGTTGTTCAAGCCt                                                                                                                          >  2:1575689/1‑150 (MQ=255)
                        aaaTCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTg                                                                                                                        >  2:539468/1‑150 (MQ=255)
                           tCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGc                                                                                                                     >  1:590015/1‑150 (MQ=255)
                           tCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGc                                                                                                                     <  1:129968/150‑1 (MQ=255)
                            cAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCg                                                                                                                    <  1:223290/150‑1 (MQ=255)
                               gTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGtt                                                                                                                 >  2:886238/1‑150 (MQ=255)
                                 tCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATGAGCGCGATGTTGTTCAAGCCTTGAGCGGTTAc                                                                                                               >  1:805193/1‑150 (MQ=255)
                                   cAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGa                                                                                                            >  1:1108278/1‑150 (MQ=255)
                                     gCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGaa                                                                                                           >  1:1356433/1‑150 (MQ=255)
                                       cccAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGaaaa                                                                                                         <  2:1500199/150‑1 (MQ=255)
                                        ccAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGaaaaa                                                                                                        >  1:1561550/1‑150 (MQ=255)
                                          agcagcTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAac                                                                                                      <  1:217156/150‑1 (MQ=255)
                                               cTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATCAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACcgc                                                                                                 <  2:805193/150‑1 (MQ=255)
                                                 gcgcAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGa                                                                                               >  2:370615/1‑150 (MQ=255)
                                                   gcAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATc                                                                                             >  1:1427865/1‑150 (MQ=255)
                                                     aaaGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCtt                                                                                           <  2:383950/150‑1 (MQ=255)
                                                        gCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGtt                                                                                        <  2:1195947/150‑1 (MQ=255)
                                                        gCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGtt                                                                                        <  2:10060/150‑1 (MQ=255)
                                                          cGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTaa                                                                                      >  1:1354871/1‑150 (MQ=255)
                                                                 gACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACAtcgt                                                                               >  1:805391/1‑150 (MQ=255)
                                                                        aGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGtcggatcg                                                                        <  2:1324794/150‑1 (MQ=255)
                                                                         cgTGATAGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGtcggatcgg                                                                       >  2:1245457/2‑150 (MQ=255)
                                                                           tGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGAc                                                                     <  2:1601257/150‑1 (MQ=255)
                                                                           tGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTAAATTAACGATAAGAGCGATGTTGTTAAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGacggatcggaa                                                                     >  1:663535/1‑149 (MQ=255)
                                                                              tCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGa                                                                  >  2:1427640/1‑150 (MQ=255)
                                                                               cGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGAt                                                                 <  1:57885/150‑1 (MQ=255)
                                                                                 tCTATAAGGTCAACGGTGATCATTGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTa                                                               <  1:698443/150‑1 (MQ=255)
                                                                                 tCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTa                                                               <  2:671952/150‑1 (MQ=255)
                                                                                 tCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTa                                                               >  1:942927/1‑150 (MQ=255)
                                                                                  cTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTAc                                                              <  1:1338883/150‑1 (MQ=255)
                                                                                  cTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTAc                                                              >  2:1062178/1‑150 (MQ=255)
                                                                                    ataAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGc                                                            <  2:619772/150‑1 (MQ=255)
                                                                                     taAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCt                                                           >  1:1228279/1‑150 (MQ=255)
                                                                                      aaGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTg                                                          >  1:443279/1‑150 (MQ=255)
                                                                                          gCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACt                                                      >  1:326014/1‑150 (MQ=255)
                                                                                          gCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACt                                                      >  2:128335/1‑150 (MQ=255)
                                                                                               ggTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACg                                                <  1:597766/150‑1 (MQ=255)
                                                                                                  gATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGAc                                              >  2:1288792/1‑150 (MQ=255)
                                                                                                   aTCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             <  1:601699/150‑1 (MQ=255)
                                                                                                   aTCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             >  2:972301/1‑150 (MQ=255)
                                                                                                   aTCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             <  2:942927/150‑1 (MQ=255)
                                                                                                   aTCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             <  2:831971/150‑1 (MQ=255)
                                                                                                    tCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACAt                                            >  1:117669/1‑150 (MQ=255)
                                                                                                     cATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATc                                           >  1:1236917/1‑150 (MQ=255)
                                                                                                      aTGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACCCCGCGATCTTGTTAAACATCGTCGGATCGGACTGCTTACGCTGCACTTTCACGACATCt                                          >  2:1126338/1‑150 (MQ=255)
                                                                                                      aTGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             <  2:1318527/147‑1 (MQ=255)
                                                                                                      aTGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACa                                             >  1:1318527/1‑147 (MQ=255)
                                                                                                       tGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGtcggatcg                                                                        >  1:324089/1‑119 (MQ=255)
                                                                                                       tGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGtcggatcg                                                                        <  2:324089/119‑1 (MQ=255)
                                                                                                          tCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCa                                      <  1:370615/150‑1 (MQ=255)
                                                                                                           cTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAg                                     >  2:793948/1‑150 (MQ=255)
                                                                                                            ttGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGc                                    <  2:1116091/150‑1 (MQ=255)
                                                                                                            ttGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGc                                    <  1:1427640/150‑1 (MQ=255)
                                                                                                             tGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCt                                   >  2:30402/1‑150 (MQ=255)
                                                                                                             tGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCt                                   <  2:296329/150‑1 (MQ=255)
                                                                                                             tGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCt                                   <  2:290487/150‑1 (MQ=255)
                                                                                                               tCCTGGCGGGTTAAAAAAATAGGCGCTTACCTTTACGATAAGCGGGATGTTGTTCAAGCCTTGAGCGGTTTCTGAATATACACCGCGATCTTGTTAAACATCGTCGGTTCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTa                                 <  2:930056/150‑1 (MQ=255)
                                                                                                                ccTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTAt                                >  2:295217/1‑150 (MQ=255)
                                                                                                                  tGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATCAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCg                              <  2:1051781/150‑1 (MQ=255)
                                                                                                                  tGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCg                              >  2:298101/1‑150 (MQ=255)
                                                                                                                  tGGCGGGTAAAAAAAATACGCGCTTACCTTAACCATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCg                              >  2:174800/1‑150 (MQ=255)
                                                                                                                     cGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGAtt                           >  1:957417/1‑150 (MQ=255)
                                                                                                                       ggTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTg                         <  2:1356433/150‑1 (MQ=255)
                                                                                                                         tAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGc                       >  2:542682/1‑150 (MQ=255)
                                                                                                                          aaaaaaaaTACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCt                      <  2:1142762/150‑1 (MQ=255)
                                                                                                                           aaaaaaaTACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCtt                     >  1:814745/1‑150 (MQ=255)
                                                                                                                            aaaaaaaaCGAGCTTACCTTAACGATAAGCGCCATGTTGATCAACCCTTGAGCGGTTACTGAAAAAACACAGCGATCTTGTTAAACATCGTCGGATCGGACTGATTAAGCTGCAATTTCACGACAACTTCCAGCTTATCGATTTGGCTTa                    >  2:26144/1‑150 (MQ=255)
                                                                                                                            aaaaaaTTCGCGCTTACCTTAACTATAGCCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTa                    <  1:1572005/150‑1 (MQ=255)
                                                                                                                             aaaaaTACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTat                   <  2:1228279/150‑1 (MQ=255)
                                                                                                                               aacatcgcgCTTATCGTTAAGGTAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTatca                 >  1:394089/6‑150 (MQ=255)
                                                                                                                               aaaTACGCGGTTACCTTAACGATAAGCGCGATGTTGTTCAAGCGTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTatca                 <  1:1390941/150‑1 (MQ=255)
                                                                                                                               aaaTACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTatca                 <  1:972301/150‑1 (MQ=255)
                                                                                                                                acatcgcgCTTATCGTTAAGGTAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTatcat                >  1:775441/5‑150 (MQ=255)
                                                                                                                                 catcgcgCTTATCGTTAAGGTAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTatcatc               >  1:332511/4‑150 (MQ=255)
                                                                                                                                  tACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCt              <  2:1108278/150‑1 (MQ=255)
                                                                                                                                  atcgcgCTTATCGTTAAGGTAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCt              >  2:844360/3‑150 (MQ=255)
                                                                                                                                    cgcgCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGc            <  2:590015/150‑1 (MQ=255)
                                                                                                                                    cgcgCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGc            <  1:155715/150‑1 (MQ=255)
                                                                                                                                        cTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGGCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCAGCTGCTCCa        <  2:1063946/150‑1 (MQ=255)
                                                                                                                                         ttACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAg       <  2:1427865/150‑1 (MQ=255)
                                                                                                                                           aCCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCTGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGAc     <  1:994131/150‑1 (MQ=255)
                                                                                                                                           aCCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGAc     >  2:998416/1‑150 (MQ=255)
                                                                                                                                              ttAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCt  <  1:533958/150‑1 (MQ=255)
                                                                                                                                              ttAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCt  <  2:1561550/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCCTGGCGGGTAAAAAAAATACGCGCTTACCTTAACGATAAGCGCGATGTTGTTCAAGCCTTGAGCGGTTACTGAAAAAACACCGCGATCTTGTTAAACATCGTCGGATCGGACTGATTACGCTGCACTTTCACGACATCTTCCAGCTTATCGATTTGGCTTATCATCTGCTCCAGACGCT  >  NC_000913/3850623‑3850914

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